hdu 1299

主题思想： 核心技术点是，求一个数的素因子个数。
x 、y、n都是正整数，并且 显然，x >= n ， y >= n ，现在假设 y = n +k (k为正整数) ，那么带入公式，可以得出 x = (n*(n+k))/k = n*n/k + n; 由于x 是正整数，现在的关键问题就是要求出 n*n/ k 有多少组正整数的可能，显然，所要求的就是 n*n 因子的个数// 问题已经非常接近答案了，但是最后还有一个问题，n<= 10^9 ， 那么n*n <= 10^18 ,对于一个这么大的数字怎样才能求出它因子的个数呢？

命题1： 一个正整数 n 可以用素因子唯一表示为 p1^r1 * p2^r2 * … pk^rk (其中 pi 为素数) , 那么这个数的因子的个数就是，(r1+1)(r2+1)…*(rk+1).

如果一个数字 n = p1^r1 * p2^r2 * … pk^rk ，那么 n*n = p1^r1 * p2^r2 * … pk^rk * p1^r1 * p2^r2 * … pk^rk ，它的因子的个数就是 (2*r1+1)(2*r2+1)…*(2*rk+1).

这个问题就转化成了求 n <= 10^9 的素因子的问题，只需要先筛选出 sqrt(10^9) 内的素数，然后用n去试除 sqrt(n) 中的每一个即可，当然，n可能会有大于sqrt(n) 的因子。

最后注意当n为质数时根据公式要乘3，（因为到最后n>1时还要乘以最后一项，2*1+1=3），考虑去掉重复的最后答案除2加1

参考博客：
http://fudq.blog.163.com/blog/static/19135023820126245711462/
http://www.aichengxu.com/other/646339.htm

#include <iostream>
#include<cstdio>
#include<cstring>

using namespace std;

const int INF=40000;

const int maxn=40000;
int prim[maxn];

bool isPrime(int n){

    if(n==2||n==3||n==5||n==7) return true;

    if(n%2==0||n%3==0||n%5==0||n%7==0) return false;


    for(int i=2;i*i<=n;i++){

        if(n%i==0) return false;
    }

    return true;
}

int Index=0;
void init(){

    for(int i=2;i<=INF;i++){

        if(isPrime(i)) prim[Index++]=i;
    }

}

int main()
{
   init();

   int T;
   scanf("%d",&T);
   int t=T;
   int n;
   while(T--){


    scanf("%d",&n);

    int ans=1;

    for(int i=0;i<Index&&n>1;i++){

        int tmp=0;
        while(n%prim[i]==0){
            tmp++;
            n/=prim[i];
        }
        ans*=(1+2*tmp);
    }

    if(n>1) ans*=3;

    printf("Scenario #%d:\n",t-T);
    printf("%d\n\n",(ans+1)/2);
   }
    return 0;
}