摘要：
polya定理＋乘法逆元 1 #include 2 using namespace std; 3 #define maxn 65 4 #define maxs 25 5 int sr,sb,sg,n,m,p; 6 int a[maxn][maxn],f[maxs][maxs][maxs],vis... 阅读全文
摘要：
发现其实是一个下凸壳，所以先按斜率排序，然后判断当前直线与栈顶直线的交点是否更靠右 注意平行的情况 1 #include 2 using namespace std; 3 4 #define maxn 50005 5 #define esp 1e-8 6 struct node{ 7 ... 阅读全文
摘要：
[BZOJ1005]明明的烦恼 百度题解第一篇简直清晰得飞起，orz 懒得写高精度除法，因为结果一定是整数，所以可以分解质因数乱搞 1 #include 2 using namespace std; 3 #define mod 1000000 4 #define maxp 1000 5 int ... 阅读全文
摘要：
题意：求圆Ｏ: x^2+y^2=r^2(r>0)上坐标为整点的个数 移向项 y^2=r^2-x^2=(r-x)(r+x) 设d=gcd(r-x,r+x)得 y^2=(d^2)*(r-x)/d*(r+x)/d 设A=(r-x)/d,B=(r+x)/d，得A+B=2*r/d 因为A，B为整数，所... 阅读全文
摘要：
做tsinsen发现题解都舔不懂QAQ，只能学习一个bzoj1070修车 吃了hzwer的安利，是为下一道题做铺垫 n辆车，m个修车工，把修车工拆成n个，每辆车分别向n*m个修车工连边，第i辆车向第(j-1)*n+k个修车工连的边边权为t[i][j]*k，代表第j个修车工修的倒数第k辆车... 阅读全文
摘要：
只能过前60%，还在罗 1 #include 2 using namespace std; 3 4 #define maxn 550000 5 #define lson rt>1;39 build(lson),build(rson);40 tree[rt]=push_up(tree... 阅读全文
摘要：
题解已经讲的很清楚了，就不鬼扯了，直接上代码 1 #include 2 #define maxn 100005 3 double a[maxn],p[maxn]; 4 double qp(double base,int x){ 5 double ans=1; 6 while(x){ ... 阅读全文
摘要：
设si表示表示每个格子的状态(１表示腐地，０表示土地)，xi表示对每个格子的操作状态(１表示发射，０表示不发射)，可以列出异或方程si＾sigma(xj)=1然而介四Ｏ(n^6)，BOOM于是我们把格子分成横向(Yl)和纵向(Yr)的块，每一个块代表能对其中的格子产生影响的格子的集合，所以我们有，，... 阅读全文