当前位置:编程学习 > C/C++ >>

hdu 3187 (欧拉函数+dfs)

 


/*

给出一个数n,满足P(k)=n,其中k的素因子个数<=3;n<2^31;
    欧拉函数:
        phi(k)=k*(1-1/p1)(1-1/p2)(1-1/p3)
              =(p1-1)*p1^x * (p2-1)*p2^y * (p3-1)*p3^z;
    phi(1)=1;
    因为n=phi(k),枚举pi-1,if(n%(pi-1)==0&&is_prime(pi))则加入pi;
    然后用DFS,最大有三个数使phi(k)=n;
    注意phi(1)=1;
    当n==1时,有两种情况,k=1,2;


    注意这里枚举(pi-1),不是n的素因子;

 

 

 

2013/04/22-13:36

*/

 

 

[cpp]
#include"stdio.h"  
#include"string.h"  
#include"algorithm"  
using namespace std; 
int n,ans; 
int A[4010]; 
int cnt; 
int ret; 
int t[3]; 
int is_prime(int x) 

    int i; 
    for(i=2;i*i<=x;i++) 
        if(x%i==0)return 0; 
        return 1; 

//从第x个开始,已经找到y个,总共需要找z个。。。  
void dfs(int x,int y,int z) 

    int i;     
    if(y==z) 
    { 
        int tem=n; 
        int f=0; 
        for(i=0;i<z;i++) 
        { 
            if(tem%(t[i]-1)!=0) 
            { 
                f=1;break; 
            } 
            tem/=(t[i]-1); 
        } 
        if(f==0) 
        { 
            for(i=0;i<z;i++) 
            { 
                while(tem%t[i]==0) 
                    tem/=t[i]; 
            } 
            if(tem==1)ret++; 
        } 
    } 
    else 
    { 
        for(i=x;i<cnt;i++) 
        { 
            t[y]=A[i]; 
            dfs(i+1,y+1,z); 
        } 
    } 

 
int main() 

    while(scanf("%d",&n)!=-1) 
    { 
        if(n==1) 
        { 
            printf("2\n"); 
            continue; 
        } 
        else 
        { 
            int i,j; 
            cnt=0; 
            //将n的素因子存入A中。  
            for(i=1;i*i<=n;i++) 
            { 
                if(n%i==0&&is_prime(i+1)==1) 
                    A[cnt++]=i+1; 
                if(n%i==0&&is_prime(n/i+1)==1) 
                    A[cnt++]=n/i+1; 
            } 
            sort(A,A+cnt); 
            //可能有重复的。。。  
            for(i=1;i<cnt;i++) 
            { 
                if(A[i]==A[i-1]) 
                { 
                    for(j=i;j<cnt-1;j++) 
                        A[j]=A[j+1]; 
                    cnt--; 
                } 
            } 
            ans=0; 
            for(i=1;i<=3;i++) 
            { 
                ret=0; 
                dfs(0,0,i); 
                ans+=ret; 
            } 
            printf("%d\n",ans); 
        } 
    } 
    return 0; 

#include"stdio.h"
#include"string.h"
#include"algorithm"
using namespace std;
int n,ans;
int A[4010];
int cnt;
int ret;
int t[3];
int is_prime(int x)
{
    int i;
    f

补充:软件开发 , C++ ,
CopyRight © 2022 站长资源库 编程知识问答 zzzyk.com All Rights Reserved
部分文章来自网络,