数论
huanzhizun
这个作者很懒,什么都没留下…
展开
-
中国剩余定理模板
#include #include #include using namespace std; int x, y,q; int exgcd(int a,int b) { int p; if(b==0) { x=1; y=0; return a; } q=exgcd(b,a%b); p=x; x=原创 2014-05-20 12:35:39 · 704 阅读 · 0 评论 -
扩展欧几里得
int x,y; int exgcd(int a,int b) { if(b==0) { x=1; y=0; return a; } int p,t; p=exgcd(b,a%b); t=x; x=y; y=t-a/b*y; return p; }原创 2014-04-29 09:28:01 · 496 阅读 · 0 评论 -
莫比乌斯反演
#include #include #include #include #define max 100000 using namespace std; int u[max+20]; long long f[max+20]; long long g[max+20]; bool vis[max+20]; int mmax(int a,int b) { int p; p=a; i原创 2014-04-29 08:42:33 · 565 阅读 · 0 评论 -
Acdream 1148 GCD SUM 莫比乌斯
公式很好推,s1=sigma(u[i]*(m/i)*(n/i),s2=sigma(u[i]*n/i*(m/i+1)*m/i*i/2),s3=sigma(u[i]*原创 2014-08-04 23:50:58 · 552 阅读 · 0 评论