1. 程式人生 > >Pairs Forming LCM (素數,唯一分解定理)

Pairs Forming LCM (素數,唯一分解定理)

題目來源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); } }