LightOJ - 1236 Pairs Forming LCM（唯一分解定理）

最新推荐文章于 2017-08-25 11:29:47 发布

zhuxiyulu

最新推荐文章于 2017-08-25 11:29:47 发布

阅读量227

点赞数

分类专栏：数论基础

本文链接：https://blog.csdn.net/zhuxiyulu/article/details/75715925

版权

数论基础专栏收录该内容

50 篇文章 0 订阅

订阅专栏

#include<cstdio>
#include<cstring>
using namespace std;
typedef long long LL;
/*
唯一分解定理的运用
有多少组(i,j)(i<=j<=n)使得lcm(i,j)=n (1<=n<=e^14)

i=p1^a1 * p2^a2 *……* pk^ak  
j=p1^b1 * p2^b2 *……* pk^bk
gcd(i,j)=p1^min(a1,b1) * p2^min(a2,b2) *……* pk ^ min(ak,bk)
lcm(a,b)=p1^max(a1,b1) * p2^max(a2,b2) *……* pk ^ max(ak,bk)
n = p1^s1 * p2^s2 *……* pk^sk
if lcd(i,j)==n_ than max(a1,b1)==s1 …… max(ak,bk,)=sk;
if a1==s1 than (0<=b1<=s1) s1+1 种选择  if b1==s1 同理
则，总有2*(s1+1)种选择，但是，a1==b1==s1重复，那么有(2*s1+1)种选择
当 i<j ans=ans*(2*s1+1)*(2*s2+1)*……*(2*sk+1)/2 (i<j 只有一半)
当 i==j 只有一种情况 ans+=1 

*/
const int maxn=1e6+5;
const int maxm=1e7+5;
int prime[maxn];//素数
bool is_prime[maxm];//i是否是素数
int p;//素数个数
//素数筛
void sieve(int n)
{
    p=0;
    for(int i=0; i<n; i++) is_prime[i]=1;
    is_prime[0]=is_prime[1]=0;
    for(int i=2; i<n; i++)
    {
        if(is_prime[i])
        {
            prime[p++]=i;
            for(int j=2*i; j<n; j+=i)
            {
                is_prime[j]=0;
            }
        }
    }
}

LL solve(LL a)
{
    int cnt;
    LL ans=1;
    LL t=a;
    //运用唯一分解定理进行分解
    for(int i=0; i<p&&prime[i]*prime[i]<=t; i++)
    {
        if(t%prime[i]==0)
        {
            cnt=0;//计算质因数i的指数
            while(t%prime[i]==0)
            {
                t/=prime[i];
                cnt++;
            }
            ans*=(2*cnt+1);//计算i!=j的数量
        }
    }
    if(t>1)//如果a是质数
    {
        ans*=(2*1+1);
    }
    ans=(ans>>1)+1;// i<j -> ans>>1 i==j -> ans+1 
    return ans;
}

int main()
{
    sieve(maxm);
    int T;
    LL n;
    scanf("%d",&T);
    int kase=0;
    while(T--)
    {
        scanf("%lld",&n);
        printf("Case %d: %lld\n",++kase,solve(n));
    }
    return 0;
}