数论
关注
分享
扫一扫
huanzhizun
这个作者很懒，什么都没留下…
展开
-
中国剩余定理模板
#include#include#includeusing 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 · 729 阅读 · 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 · 515 阅读 · 0 评论 -
莫比乌斯反演
#include#include#include#include#define max 100000using 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 · 587 阅读 · 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 · 574 阅读 · 0 评论