欧拉函数的应用,以后看到互质的数第一个就要想到欧拉函数。今天又学到了好多家伙。
欧拉定理:
欧拉定理表明,若n,a为正整数,且n,a互质,(a,n) = 1,则a^φ(n) ≡ 1 (mod n)
费马小定理:
且(a,p)=1,那么 a^(p-1) ≡1(mod p) 假如p是质数,且a,p互质,那么 a的(p-1)次方除以p的余数恒等于1 。
筛选法求欧拉函数,时间复杂度O(nloglogn), CODE:
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <math.h>
using namespace std;
const int SIZE = 1001;
int phi[SIZE];
void init()
{
int i, j;
memset(phi, 0, sizeof(phi));
phi[1] = 1;
for(int i = 2; i < SIZE; i++) if(!phi[i])
{
for(j = i; j < SIZE; j+=i)
{
if(!phi[j]) phi[j] = j;
phi[j] = phi[j] / i * (i-1);
}
}
return ;
}
int main()
{
init();
int n;
while(~scanf("%d", &n))
{
for(int i = 1; i <= n; i++)
{
printf(i != n?"%d ":"%d\n", phi[i]);
}
}
return 0;
#include <stdlib.h>
#include <string.h>
#include <math.h>
using namespace std;
const int SIZE = 1001;
int phi[SIZE];
void init()
{
int i, j;
memset(phi, 0, sizeof(phi));
phi[1] = 1;
for(int i = 2; i < SIZE; i++) if(!phi[i])
{
for(j = i; j < SIZE; j+=i)
{
if(!phi[j]) phi[j] = j;
phi[j] = phi[j] / i * (i-1);
}
}
return ;
}
int main()
{
init();
int n;
while(~scanf("%d", &n))
{
for(int i = 1; i <= n; i++)
{
printf(i != n?"%d ":"%d\n", phi[i]);
}
}
return 0;
}
求一个整数的欧拉函数值,CODE:
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <math.h>
using namespace std;
int euler_phi(int n)
{
int m = floor(sqrt(n+0.5));
int ans = n;
for(int i = 2; i <= m; i++) if(n%i == 0)
{
ans = ans / i * (i-1);
while(n%i == 0)
{
n /= i;
}
}
if(n > 1) ans = ans / n *(n-1);
return ans;
}
int main()
{
int n;
while(~scanf("%d", &n), n)
{
printf("%d\n", euler_phi(n));
}
return 0;
}
#include <stdlib.h>
#include <string.h>
#include <math.h>
using namespace std;
int euler_phi(int n)
{
int m = floor(sqrt(n+0.5));
int ans = n;
for(int i = 2; i <= m; i++) if(n%i == 0)
{
ans = ans / i * (i-1);
while(n%i == 0)
{
n /= i;
}
}
if(n > 1) ans = ans / n *(n-1);
return ans;
}
int main()
{
int n;
while(~scanf("%d", &n), n)
{
printf("%d\n", euler_phi(n));
}
return 0;
}
题目大意:求小于n并且与n不互质的数的和。
思路:由于数据范围比较大,显然需要用到__int64,暴力解决是不可能的。由于在poj做了一道求小于n且与n互质的数的和,所以可以直接总和减一下就行了。
设小于N且与N互质的正整数之和, 设为ans.
不妨设这些数为a[0],a[1], a[2], ..., a[ phi(N)-1 ],
由gcd(N, a[i]) =1,那么gcd(N,N - a[i]) =1
这样, N - a[0], N - a[1], ..., N - a[ phi(N)-1]与原数列相同, 从而:
ans = a[0] + a[1] + ... + a[ phi(N) -1]
ans= (N - a[0]) + (N - a[1]) + ... + (N - a[ phi[N]-1 ]);两式相加得 ans=N*phi[N]/2; //那么结果就显然了res=(N-1)*N/2-ans;
CODE:
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <math.h>
using namespace std;
#define MOD 1000000007;
__int64 euler_phi(__int64 n)
{
int m = floor(sqrt(n+0.5));
__int64 ans = n;
for(int i = 2; i <= m; i++) if(n%i == 0)
{
ans = ans / i * (i-1);
while(n%i == 0)
{
n /= i;
}
}
if(n > 1) ans = ans / n * (n-1);
return ans;
}
int main()
{
__int64 n, dif;
while(~scanf("%I64d", &n), n)
{
dif = ((n*(n-1)-n*euler_phi(n))/2)%MOD;
while(dif < 0) dif += MOD;
printf("%I64d\n", dif);
}
return 0;
#include <stdlib.h>
#include <string.h>
#include <math.h>
using namespace std;
#define MOD 1000000007;
__int64 euler_phi(__int64 n)
{
int m = floor(sqrt(n+0.5));
__int64 ans = n;
for(int i = 2; i <= m; i++) if(n%i == 0)
{
ans = ans / i * (i-1);
while(n%i == 0)
{
n /= i;
}
}
if(n > 1) ans = ans / n * (n-1);
return ans;
}
int main()
{
__int64 n, dif;
while(~scanf("%I64d", &n), n)
{
dif = ((n*(n-1)-n*euler_phi(n))/2)%MOD;
while(dif < 0) dif += MOD;
printf("%I64d\n", dif);
}
return 0;
}