1. 程式人生 > >poj 2917 Diophantus of Alexandria 因數分解解1/x+1/y=1/n

poj 2917 Diophantus of Alexandria 因數分解解1/x+1/y=1/n

題意:

給定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);
	}
}