欧拉公式,求解小于某数且与该数互质的数的个数

74 篇文章 1 订阅


 

欧拉函数
把n分解质因数=p1^(a1)*p2^(a2)……pn^(an)
则小于n且与n互质的数的个数
Φ=p1^(a1-1)*(p1-1)*p2^(a2-1)*(p2-1)……pn^(an-1)*(pn-1)
int eular(int n)
{
	int ret=1,i;
	for (i=2;i*i<=n;i++)
		if (n%i==0){
			n/=i,ret*=i-1;
			while (n%i==0)
				n/=i,ret*=i;
		}
		if (n>1)
			ret*=n-1;
		return ret;
}


 

#include <stdlib.h>
#include<stdio.h>
#define N 50000
int main()
{
    int *phi,i,j;
    char *prime;
    prime=(char*)malloc((N+1)*sizeof(char));
    prime[0]=prime[1]=0;
    for(i=2;i<=N;i++)
        prime[i]=1;
    for(i=2;i*i<=N;i++)
    {
        if(prime[i])
        {
            for(j=i*i;j<=N;j+=i)
                prime[j]=0;
        }
    } //这段求出了N内的所有素数 
    phi=(int*)malloc((N+1)*sizeof(int));
    for(i=1;i<=N;i++)
        phi[i]=i;
    for(i=2;i<=N;i++)
    {
        if(prime[i])
        {
            for(j=i;j<=N;j+=i)
                phi[j]=phi[j]/i*(i-1); //此处注意先/i再*(i-1),否则范围较大时会溢出
        }
    }
 int t;
 scanf("%d",&t);
 while(t--){
  int sum=0;
  int n;
  scanf("%d",&n);
  for(i=1;i<=n;i++){
   sum+=phi[i];
  }
   printf("%d\n",sum-1);

 }
    return 0;
}


 

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值