数论
文章平均质量分 53
不可不戒
这个作者很懒,什么都没留下…
展开
-
hdu1905 Pseudoprime numbers
#include int isPrime(int x) { int i; if(x==2) return 1; if(x%2==0) return 0; for(i=3;i*i<=x;i+=2) if(x%i==0) return 0; return 1; } long long func(long long a,int p) { // 计算n的p次方 lo原创 2013-07-24 10:34:19 · 867 阅读 · 0 评论 -
hdu2817 A sequence of numbers (a^i mod n)
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=2817 #include #define MOD 200907 _int64 pow_mod(_int64 a,_int64 k) { _int64 temp=1; while(k) { if(k&1) temp=(temp*a)%MOD; a=(a*a)%MOD;原创 2013-08-31 09:31:58 · 613 阅读 · 0 评论 -
欧拉函数
在数论,对正整数n,欧拉函数是少于或等于n的数中与n互质的数的数目。此函数以其首名研究者欧拉命名,它又称为Euler's totient function、φ函数、欧拉商数等。 例如φ(8)=4,因为1,3,5,7均和8互质。 从欧拉函数引伸出来在环论方面的事实和拉格朗日定理构成了欧拉定理的证明。 简介 φ函数的值 φ(1)=1(唯一和1互质的数就是1本身)。转载 2013-07-16 20:01:48 · 817 阅读 · 0 评论 -
hdu1420 Prepared for New Acmer (a^i mod n)
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=1420 #include int main() { _int64 a,b,c,temp; int test; scanf("%d",&test); while(test--) { scanf("%I64d %I64d %I64d",&a,&b,&c); a%=c; t原创 2013-09-02 22:09:00 · 907 阅读 · 0 评论 -
hdu4715 Difference Between Primes
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=4715 #include #include #define MAXN 1010000 #define SIZE 1010002 int prime[SIZE],valid[SIZE]; void getPrime()//素数筛选O(N) { int i,j,原创 2013-09-09 14:17:53 · 706 阅读 · 0 评论 -
hdu1286 找新朋友 (欧拉函数)
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=1286 题解:其中会长的号码是N号,凡是和会长是老朋友的,那么该会员的号码肯定和N有大于1的公约数,否则都是新朋友:求小于n的数中与n互质的数的个数,欧拉函数。 #include #include #define MAXN 32770 #define MAX 32768原创 2013-09-03 16:44:36 · 631 阅读 · 0 评论 -
hdu1098 Ignatius's puzzle
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=1098 题解:用数学归纳法,设f(x)能整除65,要证明f(x+1)能整除65的话就只需要证明f(x+1)-f(x)能整除65...... #include int main() { int k,a,i; while(scanf("%d",&k)!=EOF) { a原创 2013-10-01 16:49:00 · 685 阅读 · 0 评论 -
hdu4506 小明系列故事——师兄帮帮忙 (快速幂)
题目连接:http://acm.hdu.edu.cn/showproblem.php?pid=4506 #include #define MOD 1000000007 _int64 arr[10001]; _int64 pow_mod(_int64 x,int k) {//快速幂x^k%mod _int64 temp=1; while(k) { if(k&1)原创 2013-09-22 16:26:14 · 742 阅读 · 0 评论 -
hdu4548 美素数
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=4548 #include #include #define MAXN 1000000 #define SIZE 1000002 int prime[SIZE],valid[SIZE],ans[SIZE]; void getPrime()//素数筛选O(N) {原创 2013-09-22 19:11:30 · 701 阅读 · 0 评论 -
hdu2157 How many ways?? (矩阵快速幂)
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=2157 题解:见:http://www.matrix67.com/blog/archives/276 #include #include #define MAXN 21 #define MOD 1000 typedef struct Matrix {原创 2013-08-28 19:01:10 · 570 阅读 · 0 评论 -
hdu2161 Primes
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=2161 #include int isprime(int n) { int i; if(n==1||n==2) return 0; for (i=2;i*i<=n;i++) { if(n%i==0) r原创 2013-08-16 21:06:50 · 1032 阅读 · 0 评论 -
hdu1262 寻找素数对
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=1262 #include #include #define MAXN 10000 #define SIZE 10002 int prime[SIZE],valid[SIZE]; void getPrime()//素数筛选O(N) { int i,j,tot=0; memset(va原创 2013-08-16 21:00:34 · 628 阅读 · 0 评论 -
hdu3123 GCC
/*n>=m,m!+(m+1)!+...+n!这些项都是可以被m整除的, 只需要找比m小n的阶乘即可*/ #include #include int main() { int test,m,i,n,tpow; _int64 ans,tp; char nstr[105]; scanf("%d",&test); while(test--) { n=0; ans=tp=1;原创 2013-07-28 20:52:31 · 696 阅读 · 0 评论 -
hdu2601 An easy problem
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=2601 //n+1=(i+1)(j+1)=i*j+i+j #include #include int main() { int test,ans; _int64 n,i,temp; scanf("%d",&test); while(test--) { ans=0; s原创 2013-08-08 22:13:56 · 660 阅读 · 0 评论 -
hdu1108 最小公倍数
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=1108 #include #include int gcd(int a,int b) { int temp; while(b!=0) { temp=a%b; a=b; b=temp; } return原创 2013-08-10 12:24:41 · 629 阅读 · 0 评论 -
hdu2138 How many prime numbers (素数测试)
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=2138 米勒拉宾算法 #include _int64 pow_mod(_int64 a,_int64 b,_int64 mod) { _int64 ans=1; while(b) { if(b&1) ans=(ans*a)%mod; a=(a*a)%mod原创 2013-08-16 09:52:06 · 624 阅读 · 0 评论 -
hdu1575 Tr A (矩阵快速幂)
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=1575 #include #include #define MAXN 11 #define MOD 9973 typedef struct Matrix { int matr[MAXN][MAXN]; }Matrix; int n; Matrix init,unit,c; Matr原创 2013-08-28 12:41:01 · 532 阅读 · 0 评论 -
hdu2824 The Euler function (欧拉函数)
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=2824 打表 #include #include #define MAXN 3000002 #define MAX 3000000 int prime[MAXN],phi[MAXN],valid[MAXN],tot; void getPrime()//素数筛选O(N) {原创 2013-08-16 10:37:51 · 602 阅读 · 0 评论 -
hdu1215 七夕节 (因子分解)
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=1215 #include int main() { int t,sum,i,n; scanf("%d",&t); while(t--) { sum=0; scanf("%d",&n); for (i=1;i*i原创 2013-08-16 17:40:40 · 577 阅读 · 0 评论 -
hdu1128 Self Numbers
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=1128 #include #include #define MAXN 1000000 #define SIZE 1000002 bool arr[SIZE]; int main() { int i,j,pre; memset(arr,false,sizeof(false)); f原创 2013-08-16 18:21:12 · 538 阅读 · 0 评论 -
hdu1395 2^x mod n = 1
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=1395 #include int main() { int n,i,x; while(scanf("%d",&n)!=EOF) { if((n&1)==0) { printf("2^? mod %d = 1\n",n); continue; } x=1原创 2013-09-25 19:06:09 · 595 阅读 · 0 评论