题目来源: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);
}
}