Divisors 欧拉函数的应用

136 篇文章 0 订阅
48 篇文章 0 订阅
#include <stdio.h>
#include <cstring>
#define mm 500
bool fuck[mm];
int wokao[mm][100];
long long ca[mm][mm];
int p[100];
int main()
{
    int t=0,x,n;
    memset(fuck,0,sizeof(fuck));
    fuck[0]=fuck[1]=0;
    for(int i=2; i<=mm; i++)
    {
        if(!fuck[i]) p[t++]=i;
        for(int j=0; j<t&&i*p[j]<mm; j++)
        {
            fuck[i*p[j]]=true;
            if(!(i%p[j])) break;
        }
    }
    memset(wokao,0,sizeof(wokao));
    for(int i=0; i<t; i++)
        for(int j=2; j<mm; j++)
            wokao[j][i]=j/p[i]+wokao[j/p[i]][i];
    for(int i=2; i<mm; i++)
        for(int j=1; j<i; j++)
        {
            ca[i][j]=1;
            for(int k=0; k<t&&wokao[i][k]; k++)
            {
                int a=wokao[i][k]-wokao[j][k]-wokao[i-j][k];
                if(a)
                {
                    ca[i][j]*=(a+1);
                }
            }
        }
    while(scanf("%d%d",&n,&x)==2)
    {
        if(x==0||x==n)
        {
            printf("1\n");
            continue;
        }
        printf("%lld\n",ca[n][x]);
    }
    return 0;
}

  • 1
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值