20160315 LightOJ 1236 Pairs Forming（素因子分解）

题意：

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;
}

求上述代码的返回值，即满足i<=j&&lcm(i,j)==n的数对的个数。

思路：

可以将n进行素因子分解，得到π(ni^ei)的形式，对于i,j可进行如下分析：

      对于每个素因子ni：

         1.当i所含ni 因子的个数为ei 时，j所含ni个数有ei+1种方案。

         2.当i所含ni 因子的个数小于ei 时，j所含ni 因子个数必须为ei。

因此所有lcm(i,j)==n的数对数为π(2*ei+1)，因为有i<=j的条件，所以要去掉重复的部分。

刚开始对于每一个数进行单独分解结果超时了，后来用打素数表的方法进行分解。

#include <iostream>
#include <stdio.h>
#include <math.h>
#include <string.h>
using namespace std;
typedef long long LL;
int maxn=1000;
bool valid[10000005];
int p[1000000];
void solve(int n,int &tot)
{
    memset(valid,true,sizeof(valid));
    for(int i=2;i<=n;i++)
    {
        if(valid[i])
        {
            tot++;
            p[tot]=i;
        }
        for(int j=1; ((j<=tot)&&(i*p[j]<=n)); j++)
        {
            valid[i*p[j]]=false;
            if(i%p[j]==0) break;
        }
    }
}

int main()
{
    int T,c=1,cnt=0;
    LL n;
    solve(10000005,cnt);
    scanf("%d",&T);
    while(T--)
    {
        scanf("%lld",&n);
        LL ans=1;
        //factor(n,f,e,cnt);
        for(int i=1;i<=cnt;i++)
        {
            if((long long)p[i]*p[i]>n) break;
            if(n%p[i]==0)
            {
                int e=0;
                while(n%p[i]==0)
                {
                    e++;
                    n/=p[i];
                }
                ans*=2*e+1;
            }
        }
        if(n!=1)
        {
            ans*=3;
        }
        printf("Case %d: %lld\n",c++,ans/2+1);
    }
    return 0;
}