Pairs Forming LCM

http://acm.hust.edu.cn/vjudge/contest/view.action?cid=109329#problem/B   

题意：在a,b中(a,b<=n)(1 ≤ n ≤ 1014),有多少组(a,b)  (a<b)满足lcm(a,b)==n;

 

先来看个知识点：

素因子分解：n = p1 ^ e1 * p2 ^ e2 *..........*pn ^ en

for i in range(1,n):

        ei 从0取到ei的所有组合

必能包含所有n的因子。

现在取n的两个因子a,b

a=p1 ^ a1 * p2 ^ a2 *..........*pn ^ an

b=p1 ^ b1 * p2 ^ b2 *..........*pn ^ bn

gcd(a,b)=p1 ^ min(a1,b1) * p2 ^ min(a2,b2) *..........*pn ^ min(an,bn)

lcm(a,b)=p1 ^ max(a1,b1) * p2 ^ max(a2,b2) *..........*pn ^ max(an,bn)

哈哈，又多了种求gcd,lcm的方法。

 

题解：

先对n素因子分解，n = p1 ^ e1 * p2 ^ e2 *..........*pk ^ ek，

lcm(a,b)=p1 ^ max(a1,b1) * p2 ^ max(a2,b2) *..........*pk ^ max(ak,bk)

所以，当lcm(a,b)==n时，max(a1,b1)==e1,max(a2,b2)==e2,…max(ak,bk)==ek

当ai == ei时，bi可取 [0, ei] 中的所有数  有 ei+1 种情况，bi==ei时同理。

那么就有2(ei+1)种取法,但是当ai = bi = ei 时有重复，所以取法数为2(ei+1)-1=2*ei+1。
除了 (n, n) 所有的情况都出现了两次  那么满足a<=b的有 (2*ei + 1)) / 2 + 1 个

#include <iostream>
#include<stdio.h>
#include <bits/stdc++.h>
#define ll  long long
using namespace std;
ll t,ans,n,top,prime[1000006];
bool isprime[10000006];
void check()
{
    for(ll i=2; i<10000006; i++)
    {
        if(isprime[i]==0)
            prime[ans++]=i;
        for(ll j=0; j<ans&&i*prime[j]<10000006; j++)
        {
            isprime[i*prime[j]]=1;
            if(i%prime[j]==0)
                break;
        }
    }
}
int main()
{
    check();
    scanf("%lld",&t);
    while(t--)
    {
        ll sum=1;
        scanf("%lld",&n);
        for(int i=0; i<ans&&prime[i]*prime[i]<=n; i++)
        {
            ll s=0;
            while(n%prime[i]==0)
            {
                s++;
                n/=prime[i];
            }
            sum*=(2*(s+1)-1);
        }
        if(n>1)
            sum*=3;
        printf("Case %lld: %lld\n",++top,(sum+1)/2);
    }
    return 0;
}

卡了几发re，想了好久才发现开了一个long long的数组1e7的，就是这里的问题，开不了这么大的，开bool的可以，一个布尔变量就一个字节，long long应该是四个，一个数组的内存有限。