码迷,mamicode.com
首页 > 其他好文 > 详细

hunnu--11550--欧拉函数

时间:2015-06-08 23:24:59      阅读:218      评论:0      收藏:0      [点我收藏+]

标签:

欧拉函数
Time Limit: 3000ms, Special Time Limit:6000ms, Memory Limit:65536KB
Total submit users: 73, Accepted users: 59
Problem 11550 : No special judgement
Problem description
  一个数x的欧拉函数Φ(x)定义为所有小于x的正整数中与x互质的数的数目,如小于5且和5互质的数有1、2、3、4,一共4个,故Φ(5)=4。

对于任意正整数x,我们定义两种操作: 
1、f(x) = x + Φ(x);
2、g(x) = x * Φ(x);

现在,给定一个数a,问从1开始,需要多少步操作能得到a。
(如,当a = 2时,f(1)即为所求,故答案为1,而当a = 3时,f(f(1))即为所求,故答案为2) 

Input
  每行输入一个整数a(0<a<=100000)。 
Output
  输出需要的步数,如果无法得到,输出-1;
Sample Input
2
3
Sample Output
1
2
Problem Source
  HUNNU Contest 

如果你不知道求欧拉函数的话那么这题还是不要看了,就像我一样~~

打欧拉表的是网上找的,我自己写的只有下面的dp

#include <cstdio>
#include <cstring>
#include <iostream>
#include <malloc.h>
#define LL long long
using namespace std;
const LL N = 100001;
LL n;
LL *phi,i,j;
    char *prime;
LL flog;
LL visit[N];
LL dp[N],num[N];
int main()
{
    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=(LL*)malloc((N+1)*sizeof(LL));
    for(i=1; i<N; i++)
    {
        dp[i]=N;//dp初始化
        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),否则范围较大时会溢出
            }
        }
    }//以上都是打表部分,表示没细心去看,导致连注释都懒得改了,原创大哥看到了别喷我~
    dp[1]=0;
    //dp的递推思路就是用小点推出大点,因为上面两个公式是递增的,所以从小到大来遍历的话,每次遍历到的点都是已经求出了最优方案的
    for(i=1;i<N;i++)
    {
        if(i+phi[i]<N&&dp[i+phi[i]]>dp[i]+1)//用公式1如果没有过界,同时当前方案用的步数更少,那么替换
        dp[i+phi[i]]=dp[i]+1;
        if(i*phi[i]<N&&dp[i*phi[i]]>dp[i]+1)//用公式2如果没有过界,同时当前方案用的步数更少,那么替换
        dp[i*phi[i]]=dp[i]+1;
    }

    while(scanf("%lld",&n)!=-1)
    {
        if(dp[n]<N)//因为dp初始化就是N,所以这一步不能少。
            printf("%lld\n",dp[n]);
        else
            printf("-1\n");
    }
    return 0;
}

 

hunnu--11550--欧拉函数

标签:

原文地址:http://blog.csdn.net/jingdianitnan/article/details/46418461

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!