Qianyri的博客

欧拉函数线性筛杜教筛莫比乌斯反演容斥原理HDU2588 GCD欧拉函数给定N,M,求1<=X<=N 且gcd(X,N)>=M条件下X的个数#include<bits/stdc++.h>using namespace std;int euler(int n){ int ans=n; for(int i=2;i*i&...

 单点更新 区间查询void add(int x,int v){ while(x<=n) c[x]+=v,x+=x&(-x);}int get(int r){ int sum=0; while(r) sum+=c[r],r-=r&(-r); return sum;}int get(int l, int r){ r...

HDU6354 Everything Has Changed求圆弧的周长#include<bits/stdc++.h>using namespace std;const int MAX=1e2+5;const double PI=3.14159265358979323846;const double eps=1e-8;struct point{ doubl...

转载一.  巴什博奕（Bash Game）：  A和B一块报数，每人每次报最少1个，最多报4个，看谁先报到30。这应该是最古老的关于巴什博奕的游戏了吧。其实如果知道原理，这游戏一点运气成分都没有，只和先手后手有关，比如第一次报数，A报k个数，那么B报5-k个数，那么B报数之后问题就变为，A和B一块报数，看谁先报到25了，进而变为20,15,10,5，当到5的时候，不管A怎么报数，最后一...

洛谷P1129匈牙利算法 求二分图最大匹配#include <bits/stdc++.h>using namespace std;const int MAX1=5*1e3+1;const int MAX2=5*1e4+1;struct Edge{ int to,next;}edge[MAX2];int head[MAX1],tot;int linker...

洛谷P2742求凸包周长#include<bits/stdc++.h>using namespace std;const int MAX=1e4+1;const double eps=1e-8;struct point{ double x,y; point(){} point(double a,double b) { ...

洛谷P2801 教主的魔法区间修改\区间查询大于C的个数#include<bits/stdc++.h>using namespace std;const int MAX=2*1e6;const int eps=1e-8;int a[MAX],aa[MAX],b[MAX],c[MAX],cc[MAX][2],block,cnt;inline void read(int...