Lightoj 1236 LCM性质 + 唯一分解定理 prime数组一定要用long long，大小不超过1e6

题意：

long long pairsFormLCM( int n ) {
    long long res = 0;
    for( int i = 1; i <= n; i++ )
        for( int j = i; j <= n; j++ )
           if( lcm(i, j) == n ) res++; // lcm means least common multiple
    return res;
}

输出这个函数的结果，1 <= n <= 10^14。

题解：

1.看到GCD，LCM，就要联想到唯一分解定理。

2.唯一分解定理：一个大于1的数n可以分解为：n=p1^a1 * p2^a2 * p3^a3 ...的形式，其中pi是质数，ai是幂次。

3.LCM(a,b)=c，假如a = p1^a1 * p2^a2 , b = p1^b1 * p2^b2  ，那么c = p1^max(a1,b1) * p2^max(a2,b2)...，这个想一想就知道了。

4.设c = p1^c1 * p2^c2，假如a1=c1，则b1可取[0,c1]，共c1+1种可能，假如b1=c1，则a1可取[0,c1]，共c1+1种可能。a1=c1,b1=c1被计算了两次，所以c的幂次为ci的质因子pi有2*ci+1种可能。假设c有x个质因子，那么共有sum=(2*c1+1)*(2*c2+1)*...*(2*cx + 1)对a和b，因为题目要求a <= b，除了a=c且b=c这种情况外，都有a<b和b<a的情况，因此最终答案是(sum+1)/2。

5.注意：prime数组一定要用long long，大小不超过1e6！！！！！！！！！！！！！！！（数组开大了RE到哭。

#include<math.h>
#include<stdio.h>
#include<algorithm>
#include<string.h>
#include<vector>
#define N 10000005
using namespace std ;
int cnt = 0 ;
bool vis[N] ;
int num = 0 ;
long long prime[1000000] ;
void init()
{
	int i , j ;
    memset(vis , 0 , sizeof(vis)) ;
    for(i = 2 ; i < N  ; i ++)
    {
    	if(vis[i])
    	   continue ;
    	prime[num ++] = i ;
    	for(j = 1 ; j * i < N  ; j ++)
    	    vis[j * i] = 1 ;
	}
}
void solve(long long x)
{
	int i , j ;
	long long temp ;
	long long ans = 1 ;
	for(i = 0 ; i < num && prime[i] * prime[i] <= x ; i ++)
	{
		temp = 0 ;
		while(x % prime[i] == 0)
		{
			x /= prime[i] ;
			temp ++ ;
		}
		ans *= (2 * temp + 1) ;
	}
	if(x > 1) //剩一个无法分解的大于sqrt(x)的质因子 
	   ans *= (2 * 1 + 1) ;
	ans = (ans + 1) / 2 ;
	printf("Case %d: %lld\n" , ++cnt , ans) ;
} 
int main()
{
	long long n ;
	int t ; 
    init() ;
	scanf("%d" , &t) ;
	while(t --)
	{
		scanf("%lld" , &n) ;
		solve(n) ;
	}
}