1/X+1/Y=1/N！的两种O(n)做法

题目大意：

求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：

待补。