hdu 1299 Diophantus of Alexandria

 1/x + 1/y = 1/n  1<=n<=10^9
给你 n 求符合要求的x,y有多少对 x<=y
// 首先 x>n 那么设 x=n+m 那么 1/y= 1/n - 1/(n+m)
//  1/y = m/(n*n+n*m) 所以满足 n*n|m 的m都是可以的 
// 问题转化成求n*n 的约数个数 因数分解 然后用求约数个数公式
// 最后答案记得 除 2 因为重复算了
#include <iostream>
#include <algorithm>
#include <queue>
#include <vector>
#include <math.h>
#include <stdio.h>
#include <string.h>
using namespace std;
#define maxm 100010
// #define maxn 47010
#define LL __int64
int prim[],p;
void getPrime(){
    int maxn=sqrt(1000000001.0);
    int i,j;
    for(i=;i<=maxn;i+=)
        prim[i]=;
    for(i=;i*i<=maxn;i+=)
        if(!prim[i])
         for(j=i*i;j<=maxn;j+=i)
           prim[j]=;
   for(i=;i<=maxn;i++)
    if(!prim[i]) prim[p++]=i;
}
int main()
{
    getPrime();
    int n,T;
    LL ans;
    int i,ct,Case=;
    scanf("%d",&T);
    while(T--){
      scanf("%d",&n);
       i=;
       ans=;
       while(i<p){
         if(n%prim[i]==){
            ct=;
            while(n%prim[i]==){
                ct++;
                n=n/prim[i];
            }
            ans=ans*(*ct+);
         }
         if(n==) break;
         i++;
       }
       if(n!=) ans=ans*;
       printf("Scenario #%d:\n",Case++);
       printf("%I64d\n\n",(ans+)/);
    }
    return ;
}