2013年11月_cugbbaoge

11月 10月 09月 08月 07月 05月 04月 03月

原创 Hdu 2089 数位dp

#include #include #include using namespace std; int m,n; int p[15]; int f[15][3]; int w; int dfs(int pos,int s,int e) { if(pos==0) { if(s==0||s==1)return 1; else return 0;

2013-11-26 13:04:07 557

原创中国剩余定理模版

#include #include #include #include using namespace std; long long exgcd(long long a,long long b,long long &x,long long &y) { if(b==0) { x=1; y=0; return a; } l

2013-11-23 13:51:39 737

原创扩展欧几里德求逆元模版

a*x=1（mod b），r为a，b的最大公约数long long exgcd(long long a,long long b,long long &x,long long &y) { if(b==0) { x=1; y=0; return a; } long long r=exgcd(b,a%b,x,y);

2013-11-23 12:30:33 645

原创欧拉函数模版

long long euler(long long n) { long long ans=1,i; for(i=2;i*i<=n;i++) { if(n%i==0) { ans*=(i-1); n/=i; while(n%i==0) {

2013-11-23 12:26:01 590

原创 uvalive 10692 欧拉定理

a^x=a^(x%phi(c)+phi(c)) (mod c) a^phi(c)=1 (mod c) #include #include #include using namespace std; //m^n%k long long eular(long long n) { long long ans=1,i; for(i=2; i*i<=n; i++) {

2013-11-21 16:13:04 654

原创 hdu 3037 lucas定理

思路：求c（m+n，n）%p lucas： a=ak*p^k+(ak-1)*p^(k-1)+(ak-2)*p^(k-2)+…+a0 b=bk*p^k+(bk-1)*p^(k-1)+(bk-2)*p^(k-2)+…+b0 C(a,b)=C(ak,bk)*C((ak-1),(bk-1))*…C(a0,b0)(modp) #include #include #include using n

2013-11-19 21:27:19 777

原创 hdu 1724 辛普森积分

思路：求积分。辛普森公式 #include #include #include #include #define eps 1e-8 using namespace std; double a,b; double F(double x) { return 2.0*b*sqrt(1.0-x*x/(a*a)); } double simpson(double l,double r

2013-11-19 17:17:42 1074

空空如也

TA创建的收藏夹 TA关注的收藏夹

TA关注的人