题意:
求使不定方程1 / x + 1 / y = 1 / n成立的正整数解有多少个。
题解:
先说一下唯一分解定理的大概概念,即任何一个大于1的数字可以被n个素数相乘来表示。其主要有两条公式供大家使用:
(1)一个整数n可以表示为若干素数乘积: n = p1^a1 * p2^a2*…*pk^ak;(其中p1,p2…pk都是n的素数,且p1
#include<stdio.h>
#include<string.h>
#include<math.h>
#include<algorithm>
using namespace std;
#define LL long long int
const int MAXN=40000+7; //40000*40000>1e9了,所以足够了
bool prime[MAXN];
int a[MAXN],m=0,n;
void Prime()
{
prime[0]=prime[1]=false,prime[2]=true;
for(int i=3;i<MAXN;i++)
{
if(i%2) prime[i]=true;
else prime[i]=false;
}
for(int i=3;i<=sqrt(MAXN);i++)
{
if(prime[i])
for(int j=i+i;j<MAXN;j+=i)
prime[j]=0;
}
for(int i=2;i<MAXN;i++)
if(prime[i])
a[m++]=i;
}
int main()
{
Prime();
int k=1,t;
scanf("%d",&t);
while(t--)
{
scanf("%d",&n);
LL ans=1;
for(int i=0;i<m&&a[i]<=n;i++)
{
int num=0;
while(n%a[i]==0)
{
num++;
n/=a[i];
}
ans*=(2*num+1);
}
if(n>1)//一定要加上这个条件,防止漏掉一个素数的情况,为什么是素数因为一个合数可以由多个素数组成,这个合数除去多个素数之后还大于1的话那肯定还是素数了。
ans*=3;
printf("Scenario #%d:\n%lld\n\n",k++,(ans+1)/2);
}
}