Ngshily

摘要： 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)的块，每一个块代表能对其中的格子产生影响的格子的集合，所以我们有，，... 阅读全文

摘要： 把题目所给的三元组看成三维空间里的点坐标，记录下来每个点在每个坐标平面上的映射，对于每一个询问，查询每个面上分别有多少个点在该点的映射与原点所组成的矩形内部 原先想到二维树状数组，然而老爷吃着鸡块说了一句空间炸了，，，于是炸了 然后发现其实还是在树状数组上lg，但是要把插入和询问和（huo）... 阅读全文

摘要： 点分治+trie树+迷之复杂度分析 调代码时候的鬼畜错误也是不想说撒了 1 #include 2 using namespace std; 3 typedef pairpii; 4 #define fir first 5 #define sec seco... 阅读全文