[容斥原理] hdu 2204 Eddy's爱好

标签：io   os   for   sp   数据   代码   amp   ef   as   

题意：

中文题目！

思路：

首先 M^k可以分解成 (M^(k*p)) p是素数

这么我们只要枚举素因子就好了

由于数据 所以只要枚举60以内的素数就够了

然后因为2*3*5*7就超过60了 做容斥原理就最多就只有三次

代码：

#include"cstdlib"
#include"cstdio"
#include"cstring"
#include"cmath"
#include"queue"
#include"algorithm"
#include"iostream"
using namespace std;
#define ll __int64
int ss[]={2,3,5,7,11,13,17,19,23,29,31,37,41,43,47,53,59,999};
int b[12];
ll n;
ll dfs(int k,int x,int lit)
{
    ll ans=0;
    if(x==lit)
    {
        int tep=1;
        for(int i=0;i<lit;i++) tep*=b[i];
        if(tep>60) return 0;
        //printf("%d %I64d\n",tep,(ll)(pow(n+0.0,1.0/tep)+0.00001)-1);
        return (ll)(pow(n*1.0,1.0/tep)+0.00001)-1;
    }
    for(int i=k+1;ss[i]<=59;i++)
    {
        b[x]=ss[i];
        ans+=dfs(i,x+1,lit);
    }
    return ans;
}
int main()
{
    while(scanf("%I64d",&n)!=-1)
    {
        ll ans=0;
        for(int i=1;i<=3;i++)
        {
            if(i%2!=0) ans+=dfs(-1,0,i);
            else ans-=dfs(-1,0,i);
        }
        printf("%I64d\n",ans+1);
    }
    return 0;
}

[容斥原理] hdu 2204 Eddy's爱好

标签：io   os   for   sp   数据   代码   amp   ef   as   

原文地址：http://blog.csdn.net/wdcjdtc/article/details/40840311