poj 2917 Diophantus of Alexandria 因數分解解1/x+1/y=1/n
阿新 • • 發佈:2019-02-18
題意:
給定n,求丟番圖方程1/x+1/y=1/n(x<=y)的解數。
分析:
設x=n+a,y=n+b,化簡可得n^2=a*b.設n^2的因子個數為p,n^2的所有因子中除n外都是成對出現的,故方程解數為(p+1)/2。
程式碼:
//poj 2917 //sep9 #include <iostream> using namespace std; int num[64]; int main() { int cases,t=0; scanf("%d",&cases); while(cases--){ int n,p=0; memset(num,0,sizeof(num)); scanf("%d",&n); for(int i=2;i*i<=n;++i){ if(n%i==0){ while(n%i==0){ ++num[p]; n/=i; } ++p; } } if(n!=1){ ++num[p]; ++p; } for(int i=0;i<p;++i) num[i]*=2; int ans=1; for(int i=0;i<p;++i) ans*=(num[i]+1); printf("Scenario #%d:\n%d\n\n",++t,(ans+1)/2); } }