[Violet]樱花

BZOJ
luogu
一道有意思的推式子题
\[\frac{1}{x}+\frac{1}{y}=\frac{1}{N!}\]
\[\frac{x+y}{xy}=\frac{1}{N!}\]
\[xy-(x+y)N!=0\]
\[(N!)^2+xy-(x+y)N!=(N!)^2\]
\[(x-N!)(y-N!)=(N!)^2\]
由于x,y为任意正整数,所以x-N!,y-N!也为任意正整数
于是有\[xy=(N!)^2\]
考虑把N!分解质因数,\(N!=p_1^{k_1}+p_2^{k_2}+p_3^{k_3}+...\)
那么\((N!)^2=p_1^{2k_1}+p_2^{2k_2}+p_3^{2k_3}+...\)
显然如果确定了x,那么y是唯一的
x有多少种可能?\((2k_1+1)(2k_2+1)(2k_3+1)...\)
然后线筛一下,\(O(nlogn)\)分解2~n即可

#include<bits/stdc++.h>
using namespace std;
const int N=1e6+5,_=1e5,p=1e9+7;
int n,cnt,ans=1,pri[_],mn[N],t[_];
bool isp[N];
void fact(int x){if(x==1)return;t[mn[x]]++;fact(x/pri[mn[x]]);
}
int main(){cin>>n;memset(isp,1,sizeof(isp));for(int i=2;i<=n;i++){if(isp[i]){pri[++cnt]=i;mn[i]=cnt;}for(int j=1;j<=cnt&&pri[j]*i<=n;j++){isp[pri[j]*i]=0;mn[pri[j]*i]=j;if(i%pri[j]==0)break;}}for(int i=2;i<=n;i++)fact(i);for(int i=1;i<=cnt;i++)ans=1ll*ans*(2*t[i]+1)%p;cout<<ans<<endl;return 0;
}

转载于:https://www.cnblogs.com/sdzwyq/p/9840346.html

[Violet]樱花相关推荐

  1. Luogu P1445[Violet]樱花/P4167 [Violet]樱花

    Luogu P1445[Violet]樱花/P4167 [Violet]樱花 真·双倍经验 化简原式: \[\frac{1}{x}+\frac{1}{y}=\frac{1}{n!}\] \[\frac ...

  2. P4167 [Violet]樱花

    题目背景 又到了一年樱花盛开的时节.Vani 和妹子一起去看樱花的时候,找到了一棵大大的樱 花树,上面开满了粉红色的樱花.Vani 粗略估计了一下,一共有足足  n! 片花瓣. Vani 轻柔地对她说 ...

  3. 【[Violet]樱花】

    就是化柿子 我们求 \[\frac{1}{x}+\frac{1}{y}=\frac{1}{n!}\] 的正整数解的个数 喜闻乐见的化柿子了 \[\frac{x+y}{xy}=\frac{1}{n!}\ ...

  4. 触觉互联网_触觉服是VR的未来吗?

    触觉互联网 Haptic suits represent the next step towards true immersion in virtual reality gaming. Virtual ...

  5. 好消息轮询_在大流行中寻找数据科学工作而不是好消息

    好消息轮询 There's good news if you've been thinking of getting into the field of data science or have be ...

  6. 【BZOJ 2721】 2721: [Violet 5]樱花 (筛)

    2721: [Violet 5]樱花 Time Limit: 5 Sec  Memory Limit: 128 MB Submit: 599  Solved: 354 Description Inpu ...

  7. bzoj2721 [Violet 5]樱花

    分析:这道题对于我这种蒟蒻来说还是很有难度啊. 思路非常巧妙,既然不定方程要有有限个数解,那么这个肯定会对解有所限制,也就是本题中的正整数.这个时候我们要表示出方程中的一个根x,设z = n!,那么x ...

  8. 1624 Violet 5 樱花(Bzoj2721 LOJ10202 LUOGU1445 提高+/省选-) 交叉相乘 阶乘质因数分解技巧 约数个数计算

    总目录 在线测评地址(ybt) 在线测评地址(LOJ) 在线测评地址(LUOGU) 1.交叉相乘 阶乘质因数分解技巧 约数个数计算 ybt 通过 测试点 结果 内存 时间 测试点1 答案正确 604K ...

  9. ikbc机械键盘打字出现重复_抖音开啤酒小姐姐同款面具,ikbc白无垢樱花机械键盘开箱体验...

    现如今机械键盘同质化越来越严重,就跟智能手机领域一样,不看看背面的LOGO的话,还真没法认出是哪家的.目前,对于绝大多数人来说,还是属于颜值党的,总之,帅就完事了.外设圈很小,愿意掏几百块买一把键盘的 ...

最新文章

  1. slave-pxc后GTID不一致
  2. 进程管理supervisor的简单说明
  3. webpack.config.js====插件purifycss-webpack,提炼css文件
  4. 数据结构:基数排序(Radix sort)
  5. 2021-01-26 粒子群优化算法(PSO)以及Matlab实现
  6. ITK:提取具有多个分量的图像通道
  7. oneinstack_OneinStack lnmp网站环境部署Namecheap SSL安全证书方法
  8. c#web服务器 虚拟目录,C#建立自己的Web服务器
  9. Oracle使用数据泵 (expdp/impdp)实施迁移
  10. 第一次总结:整体框架
  11. GDAL工具使用示例(一)
  12. 蚂蚁金服自主研发的三地五中心异地多活解决方案获金融科技创新大奖
  13. centos7 jenkins 安装
  14. 高并发、高访问量如何计算【测试评估】
  15. 毕业设计python数据分析_毕业设计:爬虫及数据分析
  16. 顶点计划:关于大学生自习的调研
  17. word2vec python实现
  18. Unity3D 动态加载资源方式
  19. 基于OCR识别引擎的识别表格文字并将结果以Excel电子表格的形式原样导出的Android客户端代码
  20. 习题6-5 巡逻机器人

热门文章

  1. 几行代码实现c++/qt程序进程单例(文件锁)
  2. HTML小项目之双色球
  3. 序列模型第二周 自然语言处理与词嵌入
  4. 鲸鱼优化算法(WOA)及其优秀变体(含MATLAB代码)
  5. GC 日志解读与分析(番外篇可视化工具)
  6. SharePoint Fundation 2013中SecurityTokenServiceApplication错误
  7. python 实现青蛙跳台阶问题
  8. 【转载】任正非致员工信公布:全面提升软件工程能力和实践
  9. 使用webuploader上传大文件
  10. 7月火影忍者服务器维护,7月2日停机更新公告