_随心所欲_

摘要： 质因数分解；牛人推导公式(1^3+2^3+……+(1+a1)^3)*……*(1^3+2^3+……+(1+ai)^3)……链接http://poj.org/problem?id=3604 1 #include 2 #include 3 #include 4 #include 5 #include 6 #include 7 using namespace std; 8 int prime[3001],m; 9 bool f[3001];10 void init()11 {12 int i,j;13 m=0;14 memset(f,0,sizeof(f));15 fo... 阅读全文

摘要： 数论水题，利用筛选法…… 1 #include 2 #include 3 #include 4 #include 5 #include 6 #include 7 #include 8 using namespace std; 9 int an[1000001];10 int main()11 {12 int i,j,k,n;13 memset(an,0,sizeof(an));14 for(i=5;i1000000)19 break;20 if(an[j]==0&&an[i]==0)21 ... 阅读全文

摘要： 高精度+大数同余求模，采用千进制能够AC……链接http://poj.org/problem?id=2635 1 #include 2 #include 3 #include 4 #include 5 #include 6 #include 7 #include 8 using namespace std; 9 char a[10000];10 int b[10000],prime[1000101],m;11 bool f[1000105];12 void init()13 {14 __int64 i,j;15 m=0;16 memset(f,0,sizeof(f));... 阅读全文

摘要： 扩展的欧几里得算法……链接http://poj.org/problem?id=2115 1 #include 2 #include 3 #include 4 #include 5 #include 6 #include 7 #include 8 using namespace std; 9 __int64 extend_gcd(__int64 a,__int64 b,__int64 &x,__int64 &y)10 {11 __int64 d;12 if(b==0)13 {14 x=1;15 y=0;16 return a... 阅读全文

摘要： 将a^b进行素数分解为a1^p1*a2^p2*……ai^pi;则因子之和为：(1+a1+a1^2+……+a1^p1)*(1+a2+a2^2+……+a2^p2)……(1+ai+ai^2+……+ai^pi)这样就方便多了，使用二分快速幂可以快速求出结果……链接http://poj.org/problem?id=1845#include#include#include#include#include#include#includeusing namespace std;int prime[10001],k,f[10001],e[10001];bool a[10001];const int M=990 阅读全文

摘要： 这个有2种方法。一种是通过枚举p的值（p的范围是从1-32），这样不会超时，再就是注意下精度用1e-8就可以了，还有要注意负数的处理…… 1 #include 2 #include 3 #include 4 #include 5 #include 6 #include 7 using namespace std; 8 int fun(double n) 9 {10 int p=32,f=1;11 if(n=1;p-=f)18 {19 double a=pow(n,1.0/p);20 long t1=a;21 long t2... 阅读全文

摘要： 不互质的中国剩余定理……链接http://acm.hdu.edu.cn/showproblem.php?pid=3579#include#include#include#include#includeusing namespace std;int n,m;int extend_gcd(int a,int b,int &x,int &y){ int m; if(b==0) { x=1; y=0; return a; } else { m=extend_gcd(b,a%b,x,y); int t=x; x=y; y=t-a/b*y; } return m;}int main(){ 阅读全文

摘要： 主要考察卡特兰数，大数乘法，除法……链接http://acm.hdu.edu.cn/showproblem.php?pid=1134#include#include#include#include#includeusing namespace std;int an[102][102];void mul(int i){ int k,j,t=0,a=4*i-2,b=i+1,temp; k=100; while(an[i-1][k]9) k--; for(j=0;j9) { t=an[i][j]/10; an[i][j]%=10; } else t=0; } while(t!=0) { an[i][ 阅读全文

摘要： 中国剩余定理……、链接http://acm.hdu.edu.cn/showproblem.php?pid=1370 1 /************************************************************************* 2 > File Name: xh.cpp 3 > Author: XINHUA 4 > Mail: 525799145@qq.com 5 > Created Time: 2013/7/11 星期四 17:46:00 新华 6 *********************************... 阅读全文

摘要： 这题用扩展的欧几里得算法#include<iostream>using namespace std;void gcd(int a,int b,int &x,int &y){ if(b==0) { x=1;y=0; } else { gcd(b,a%b,x,y); int t=x; x=y; y=t-a/b*y; }}int main(){ int T,x,y; int n,m; cin>>T; while(T--) { scanf("%d%d",&n,&m); gcd(m,9973,x,y); x*=n; printf 阅读全文