UVA-10837 一个研究课题

最新推荐文章于 2024-09-04 17:18:01 发布

bluebayou

最新推荐文章于 2024-09-04 17:18:01 发布

阅读量165

点赞数 1

分类专栏：刘汝佳紫书题目文章标签： c++ dfs

本文链接：https://blog.csdn.net/bluebayou/article/details/103773758

版权

刘汝佳紫书题目专栏收录该内容

7 篇文章 0 订阅

订阅专栏

题目大意：

给出phi(n)(phi(n)≤10⁸)，输出所有n中最小的一个。输入保证n< $2*10^8$ 。

phi(n)的定义：比n小、与n互素的数的个数。

phi(n)的公式：phi(n)= $n\displaystyle\prod_{p|n}(1-\frac{1}{p})$

唯一分解n和phi(n)，得到的结果约一下，然后设m=phi(n)，则可以求出 $\frac{n}{p~1~p~2~...p~k~}=\frac{m}{(p~1~-1)(p~2~-1)...(p~k~-1)}$ 。
通过这个规律，我们可以先将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&&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;
}