w750636248的专栏

一个主要问题就是364875103=97*3761599;#include #include #include #define M (364875103LL)#define M1 (97LL)#define M2 (M/M1)#define LL long longusing namespace std;LL fac1[M1+9],fac2[M2+9];LL mod(LL a

一个主要问题，就是a和c都等于1。。。#include #include #include #define LL long longusing namespace std;#define N 100010LL el[N];struct NUM{ LL allp; LL prime[20]; NUM(){allp=0;}}num[N];voi

#include#includeusing namespace std;int pow(int n,int m){ int ans=1,a=10; while(n) { if(n&1) ans*=a; a*=a; ans%=m; a%=m; n>>=1; } retur

能兼容多个模数不互质的《中国剩余定理》

int china(int n,int *m,int *r){ int i,M=1,ans=0; for(i=0;i<n;i++) M*=m[i]; for(i=0;i<n;i++) { int w=M/m[i],x,y; gcd(w,m[i],x,y); ans=(ans+w*x*r[i])%M; }

#include#includeusing namespace std;int gcd(int a,int b,int &x,int &y){ if(b==0) { x=1;y=0;return a; } int g=gcd(b,a%b,y,x); y-=x*(a/b); return g;}int china(int

zoj3707 Calculate Prime S，不用大数搞出来

2014牡丹江 现场赛 F zoj 3824 Fiber-optic Network树形dp + 容斥原理