Pairs Forming LCM (素數,唯一分解定理)
阿新 • • 發佈:2019-01-06
題目來源:https://vjudge.net/problem/LightOJ-1236
【題意】
求a,b的最小公倍數是n的個數。
【思路】
想到了素數,想到了唯一分解定理,但是唯獨沒有想到最小公倍數的另外一種定義(我也沒見過不是):假設有數a,b,他們的最小公倍數的值就等與兩著的素因子指數最大的相乘。例如,12與8,12可以拆成2 2 3,而8則是2 2 2,所以裡面一共包含了兩種素因子,其中2的最高次冪是3,3的最高次冪是1,乘起來就是2 2 2 3,得24。
先把n分解成指數積的形式:n=x( ei^y),那麼對於a,b來說,必須有一個數的ei的次冪等於y,另外一個小於等於y,先不考慮a,b的大小,對於每一個素因子:
在a中,若ay==y,那麼就有[0,y],共有y+1中選擇。
若ay
#include<map>
#include<stack>
#include<queue>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<iostream>
#include<string>
#define mem(a,b) memset(a,b,sizeof(a))
using namespace std;
const int INF=1e9;
typedef long long LL;
const int maxn=1e7+9;
bool vis[maxn];
int prime[maxn/10];
int main()
{
int l=0;
mem(vis,0);
for(int i=2; i<=maxn; i++)
{
if(!vis[i])
{
prime[l++]=i;
for(int j=i*2; j<=maxn; j+=i)
{
vis[j]=1;
}
}
}
int T,cases=1;
scanf("%d",&T);
while(T--)
{
LL n,ans=1;
scanf("%lld",&n);
for(int i=0; i<l; i++)
{
if(LL(prime[i]*prime[i])>n) break;
int num=0;
while(n%prime[i]==0)
{
n/=prime[i];
num++;
}
ans*=num*2+1;
}
if(n>1) ans*=1*2+1;
printf("Case %d: %lld\n",cases++,ans/2+1);
}
}