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

bzoj3643 Phi的反函数 数学 搜索

时间:2017-09-17 22:23:20      阅读:226      评论:0      收藏:0      [点我收藏+]

标签:ffffff   利用   sed   题意   color   ide   open   none   lld   

链接:http://www.lydsy.com/JudgeOnline/problem.php?id=3643

题意:已知$x$,求出$phi(n)==x$的最小$n$。

不得不说今天刷了一堆(相对而言)思博大爆搜之后来这么一道相对有思考价值的数学题真是强啊……

首先我们需要记住一个基本事实:在$int$范围内,每个数的同一个质因数不会成为这个数字的$31$次方,一个数存在的不同质因数不会超过$10$个。

理由很简单:首先$int_max==2147483647$,一定没有任何一个质数$31$次方小于这个数;至于第二个嘛……

技术分享

这个东西说明了一切……

了解到这个之后,我们再回顾一下欧拉函数的性质,最重要的一点就在于:它是一个积性函数。于是我们就可以对着一个合数进行三种分类讨论:

1)他自己就是一个质数的欧拉函数,那这个质数一定是他再加$1$;

2)他是几个质数乘积,就递归搞下去。

那么我们就爆搜就可以了。预处理出$sqrt(2147483647)$范围内的质数,然后利用这些扫来扫去,最后特判剩下的是不是质数即可。还有一个优化点:可以证明选择的质数一定是递增的,因此直接递增计算即可。

技术分享
 1 #include<iostream>
 2 #include<cstdio>
 3 #include<algorithm>
 4 #include<cstring>
 5 #include<cmath>
 6 using namespace std;
 7 const int maxn=50000;
 8 int prime[maxn],tot,n;bool notprime[maxn]={1,1};
 9 void shaisushu()
10 {
11     for(int i=2;i<=50000;i++)
12     {
13         if(!notprime[i])prime[++tot]=i;
14         for(int j=1;j<=tot&&i*prime[j]<=50000;j++)
15         {
16             notprime[i*prime[j]]=1;
17             if(!(i%prime[j]))break;
18         }
19     }
20 }
21 long long pro=1;long long ans=2147483648ll;int mudi;
22 bool isprime(int x)
23 {
24     int mubiao=(int)sqrt(x);
25     for(int i=1;prime[i]<=mubiao;i++)
26         if(!(x%prime[i]))return 0;
27     return 1;
28 }
29 void dfs(int now,int num)
30 {
31     int tpro=pro,tt=now;
32     if(pro>=ans)return;
33     if(now==1){ans=min(ans,pro);return;}
34     if(now>mudi&&isprime(now+1)==1)ans=min(ans,pro*1ll*(now+1));
35     if(prime[num]-1>mudi||now<prime[num]-1)return;
36     for(int i=num;prime[i]<=mudi+1;i++)
37     {
38         if(prime[i]-1>now)return;
39         if(!(now%(prime[i]-1)))
40         {
41             now/=(prime[i]-1);pro*=1ll*prime[i];
42             dfs(now,i+1);
43             while(!(now%prime[i]))
44             {
45                 now/=prime[i],pro*=1ll*prime[i];
46                 dfs(now,i+1);
47             }
48             pro=tpro,now=tt;
49         }
50     }
51 }
52 int haha()
53 {
54     scanf("%d",&n);mudi=(int)sqrt(n);shaisushu();
55     dfs(n,1);
56     if(ans<=0x7fffffff)printf("%lld\n",ans);else puts("-1");
57 }
58 int sb=haha();
59 int main(){;}
bzoj3643

 

bzoj3643 Phi的反函数 数学 搜索

标签:ffffff   利用   sed   题意   color   ide   open   none   lld   

原文地址:http://www.cnblogs.com/Loser-of-Life/p/7537748.html

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