题目求\(\frac{1}{x}+\frac{1}{y}=\frac{1}{n!}\)已知n, x和y的正整数解的个数

设z=\(n!\)

\(\frac{1}{x}+\frac{1}{y}=\frac{1}{z}\)

\(x=\frac{yz}{y-z}\)

令\(t=y-z\)则\(x=\frac{z(t+z)}{t}=\frac{zt+z^2}{t}=z+\frac{z^2}{t}\)

当t为整数时 x,y有整数解

故求出\(z^2\)的因数即可

运用约数和公式和线性筛即可

#include<bits/stdc++.h>
#define ll long long
using namespace std;
const ll mo=1e9+7;
ll t,ans,tot,x,n;
ll prime[500005],f[1000005];
void init(ll x)
{
t=0;
for(ll i=2;i<=x;i++)
{
if(!f[i])prime[++t]=i;
for(ll j=1;j<=t;j++)
{
if(prime[j]*i>x)break;
f[prime[j]*i]=1;
if(i%prime[j]==0) break;
}
}
}
int main()
{
scanf("%lld",&n);
init(n);ans=1;
for(ll i=1;i<=t;i++)
{
tot=0;x=n;while(x){tot=tot+x/prime[i];x/=prime[i];}
ans=ans*(tot*2+1)%mo;
}
printf("%lld",ans);
return 0;
}

最新文章

  1. hibernate学习笔记之二 基本环境搭建
  2. jqGrid插件getCol方法的一个改进
  3. ./configure,make,make install的作用(转)
  4. 设置 Firewalld 防火墙控制对系统的访问
  5. mem 族函数的实现
  6. HOOK API(二)—— HOOK自己程序的 MessageBox
  7. CAS单点登录服务器搭建
  8. 用eclipce编写 MR程序 MapReduce
  9. Windows下静态库与动态库的创建与使用
  10. Python学习笔记第二十七周(Bootstrap)
  11. Web常见安全漏洞-SQL注入
  12. leetcode — letter-combinations-of-a-phone-number
  13. Ubuntu16.04下安装破解secureCRT和secureFX的操作记录
  14. taro 知识点
  15. [SoapUI] 通过JSONAssert比较两个环境的JSON Response,定制化错误信息到Excel
  16. Spring之c3p0连接池xml配置和使用举例
  17. js全局属性 全局变量
  18. 【BZOJ4036】按位或(Min-Max容斥,FWT)
  19. K-means + PCA + T-SNE 实现高维数据的聚类与可视化
  20. MVC Partial页面的使用

热门文章

  1. 关于js执行机制的理解
  2. yarn-RM-check-by-curl
  3. javascript的数组之sort()
  4. ES6 数组
  5. 002-红黑树【B-树】、二叉查找树
  6. block,inline和inline-block概念和区别(转载)
  7. linux-grep-tail-find
  8. 160A
  9. ionic3 在ios9.0 系统下 会出现ReferenceError:Can&#39;t find variable:Intl 错误提示
  10. Django进阶之中间件