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=50100;
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+=1;long long y;
    for(;;i++){
        y=q[i];
        if(x%y==0){
            long long ans=1,k=1;
            while(x%y==0){
                x/=y,k*=y,ans+=k;
            }
            if(x==1)return ans;
            return ans*dfs(x,i);
        }
    }
}
bool doit(long long x){
    if(a==1&&x==1)return 1;
    tot1=0;long long ma=(long long)sqrt(double(x));
    for(int i=1;i<=ma;i++){
        if(x%i==0){
            q[++tot1]=i;q[++tot1]=x/i;
        }
    }sort(q+1,q+tot1+1);
    if(dfs(x,1)==a)return 1;
    return 0;
}
int main(){
    while(~scanf("%I64d%I64d%I64d",&a,&b1,&b2)){
        if(a==0&&b1==0&&b2==0){break;}
        long long w=a*b2;
        if(w%b1==0){
            if(doit(w/b1))printf("%d %I64d\n",1,w/b1);
            else printf("%d\n",0);
        }
        else printf("%d\n",0);
    }
    return 0;
}