题目大意:
给出phi(n)(phi(n)≤108),输出所有n中最小的一个。输入保证n< 2 ∗ 1 0 8 2*10^8 2∗108。
phi(n)的定义:比n小、与n互素的数的个数。
phi(n)的公式:phi(n)= n ∏ p ∣ n ( 1 − 1 p ) n\displaystyle\prod_{p|n}(1-\frac{1}{p}) np∣n∏(1−p1)
唯一分解n和phi(n),得到的结果约一下,然后设m=phi(n),则可以求出
n
p
1
p
2
.
.
.
p
k
=
m
(
p
1
−
1
)
(
p
2
−
1
)
.
.
.
(
p
k
−
1
)
\frac{n}{p~1~p~2~...p~k~}=\frac{m}{(p~1~-1)(p~2~-1)...(p~k~-1)}
p 1 p 2 ...p k n=(p 1 −1)(p 2 −1)...(p k −1)m。
通过这个规律,我们可以先将m分解,把分解出来的素数每个加上1打到一个vector里。既然是素数,那么最好加一个筛素数。
然后,我们可以递归这个vector里的元素在n的唯一分解式中的a。既然是递归,那么一定要有回溯。回溯条件:除剩下的数不是这个因数的倍数。
注意:最后可能也会除剩下,这时用判素数的方法判它就可以了。
警告:数组狂慎用此代码,里面用了很多vector和set。
#include <iostream>
#include <cmath>
#include <cstring>
#include <vector>
#include <set>
using namespace std;
int vis[14150],phi,ans;
vector<int> prime,e;
set<int>usee;
void dfs(int cur,int n,int rem){
if(cur==e.size()){
if(rem==1){ans=min(ans,n);return;}
int ok=1,tmp=rem+1;
for(int i=0;i<prime.size();i++){int k=prime[i];
if(k*k>tmp)break;
if(tmp%k==0){ok=0;break;}
}
if(ok&&!usee.count(tmp))ans=min(ans,n*tmp);
return;
}
dfs(cur+1,n,rem);
int p=e[cur];
if(rem%(p-1))return;
usee.insert(p);n*=p;rem/=p-1;
while(1){
dfs(cur+1,n,rem);
if(rem%p)break;
n*=p;rem/=p;
}
usee.erase(p);
}
int main(){
for(int i=2;i<=120;i++)if(!vis[i])
for(int j=i*i;j<=14143;j+=i)vis[j]=1;
for(int i=2;i<=14143;i++)
if(!vis[i])prime.push_back(i);
for(int o_o=1;cin>>phi&φo_o++){
e.clear();usee.clear();
for(int i=0;i<prime.size();i++){int k=prime[i];
if(k>phi)break;
if(phi%(k-1)==0)e.push_back(k);
}
ans=200000010;dfs(0,1,phi);
cout<<"Case "<<o_o<<": "<<phi<<' '<<ans<<endl;
}
return 0;
}