JZYZOJ1379天才的约数和 数论 约数和

http://172.20.6.3/Problem_Show.asp?id=1379

 

易得n=a*b2/b1；

需要注意算出n之后要判断n的约数和是否等于a，这里需要用约数和定理递归，递归前求一下约数，避免数字太大浪费时间，这种数字即使是n的复杂度也承受不起吧括弧笑，相比之下sqrt(n)就好很多了。

 

代码

 #include<iostream>
 #include<cstdio>
 #include<cstring>
 #include<algorithm>
 #include<cmath>
 using namespace std;
 const int maxn=;
 long long a,b1,b2,tot,tot1;
 long long q[maxn]={};
 long long su[maxn]={};
 long long dfs(long long x,long long i){
     i+=;long long y;
     for(;;i++){
         y=q[i];
         if(x%y==){
             long long ans=,k=;
             while(x%y==){
                 x/=y,k*=y,ans+=k;
             }
             if(x==)return ans;
             return ans*dfs(x,i);
         }
     }
 }
 bool doit(long long x){
     if(a==&&x==)return ;
     tot1=;long long ma=(long long)sqrt(double(x));
     for(int i=;i<=ma;i++){
         if(x%i==){
             q[++tot1]=i;q[++tot1]=x/i;
         }
     }sort(q+,q+tot1+);
     if(dfs(x,)==a)return ;
     return ;
 }
 int main(){
     while(~scanf("%I64d%I64d%I64d",&a,&b1,&b2)){
         if(a==&&b1==&&b2==){break;}
         long long w=a*b2;
         if(w%b1==){
             if(doit(w/b1))printf("%d %I64d\n",,w/b1);
             else printf("%d\n",);
         }
         else printf("%d\n",);
     }
     return ;
 }