题目大意:
求1/X+1/Y=1/N!的答案对数。
解题思路1:
设 m=n! ,由等式知x,y必定大于n!,所以再设 x=n!+k=m+k 带入 1/m=1/x+1/y 中化简得到y=m*m/k+m,因y为整数,所以要求k整除m*m,即k为m*m的因子,问题便转化为求n!*n!的因子个数, 设n!=p1^e1 * p2^e2 * p3^e3 *...*pk^ek,则 n!*n!= p1^(2*e1) * p2^(2*e2) *...*pk^(2*ek) 。 则因子个数sum=(2*e1+1)*(2*e2+1)*...*(2*ek+1)。
对于每个质数p,它的出现次数等于[N/p]+[N/p^2]+...,这样对于每个p是O(log_p(N))的,而1~N之间的质数有O(N/lnN)个,均摊一下就是O(n)的。
至于怎么算每个质数p出现的次数,p的次数等于[n/p]+[n/p^2]+[n/p^3]+....这个可以看成是把1~n里p的倍数的数字都除掉一个p,这样的数有[n/p]个,除完之后这[n/p]个数会变成1~[n/p],然后再除掉一个p,就是[[n/p]/p]个,然后变成1~[n/p^2],一直做下去。
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<queue>
#include<map>
#include<vector>
#include<set>
#include<ctime>
#define LL long long
#define db double
#define EPS 1e-15
#define inf 1e16
#define pa pair<int,int>
using namespace std;
LL read()
{
LL x=0,f=1;char ch=getchar();
while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
return x*f;
}
int n,cnt=0;
bool f[1000010];
int pri[1000010],num[1000010];
void pre(){
memset(f,1,sizeof(f));
int j;
for (int i=2;i<=1000000;i++) {
if (f[i]) {
j=i;
cnt++; pri[cnt]=i;
while (j+i<=1000010) {
j+=i;
f[j]=0;
}
}
}
}
int main(){
pre();
while (1){
scanf("%d",&n);
if (n==0) break;
for (int i=0;i<=1000000;i++) num[i]=0;
for (int i=1;i<=cnt && pri[i]<=n;i++){
int tmp=n;
while (true){
if(tmp==0 || pri[i]>n) break;
num[i]+=tmp/pri[i];
tmp/=pri[i];
}
}
LL ans=1;
for (int i=1;i<=cnt && pri[i]<=n;i++)
if (num[i]>=1) ans*=(2*num[i]+1);
ans=(ans-1)/2+1;
printf("%lld\n",ans);
}
}
解题思路2:
待补。