【P5253】丢番图【数论,数学】
题目大意:
题目链接:https://www.luogu.org/problem/P5253
给出nnn,求1x+1y=1n\frac{1}{x}+\frac{1}{y}=\frac{1}{n}x1+y1=n1本质不同的解。
思路:
和 樱花 那道题基本一样,所以很快就推出来了。
通分
xyx+y=n\frac{xy}{x+y}=nx+yxy=n
移项
−n(x+y)+xy=0-n(x+y)+xy=0−n(x+y)+xy=0
配方
n2−n(x+y)+xy=n2n^2-n(x+y)+xy=n^2n2−n(x+y)+xy=n2
(n−x)(n−y)=n2(n-x)(n-y)=n^2(n−x)(n−y)=n2
所以如果我们n2n^2n2分解质因数即可。然后两两配对。
若n=c1p1×c2p2×...×cnpnn=c_1^{p_1}\times c_2^{p_2}\times ...\times c_n^{p_n}n=c1p1×c2p2×...×cnpn,则n=c12p1×c22p2×...×cn2pnn=c_1^{2p_1}\times c_2^{2p_2}\times ...\times c_n^{2p_n}n=c12p1×c22p2×...×cn2pn
时间复杂度O(n)O(\sqrt{n})O(n)
代码:
#include <cstdio>
using namespace std;
typedef long long ll;const int N=1e7+10;
ll n,tot,ans,c[N];int main()
{scanf("%lld",&n);for (ll i=2;i*i<=n;i++)if (!(n%i)){tot++;for (;!(n%i);n/=i)c[tot]++;}if (n>1) c[++tot]++;ans=1;for (int i=1;i<=tot;i++)ans*=(c[i]*2+1);printf("%lld",(ans+1)/2LL);return 0;
}
【P5253】丢番图【数论,数学】相关推荐
- 丢番图生平编程java_丢番图与麦乐鸡购买问题
在MIT公开课<计算机科学与编程导论>的Assignment2中,提到了丢番图方程,并有趣地将丢番图方程和卖麦乐鸡的问题联系到了一起. 首先让我们来看看维基百科中介绍的丢番图方程. 一.丢 ...
- 丢番图与麦乐鸡购买问题
在MIT公开课<计算机科学与编程导论>的Assignment2中,提到了丢番图方程,并有趣地将丢番图方程和卖麦乐鸡的问题联系到了一起. 首先让我们来看看维基百科中介绍的丢番图方程. 一.丢 ...
- 丢番图生平编程java_丢番图
[拼音]:Diufantu [外文]:Diophantus (活动于250年前后) 希腊代数学家.对他的生平事迹人们知道得很少.但流传着一首短诗(或墓志铭),用谜语的形式叙述了他的生平:"丢 ...
- JZOJ 3158. 【JSOI2013】丢番图
Description 丢番图 是亚历山大时期埃及著名的数学家.他是最早研究整数系数不定方程的数学家之一. 为了纪念他,这些方程一般被称作丢番图方程.最著名的丢番图方程之一是 xn+yn=znx^n+ ...
- 丢番图(Diophantine)方程MATLAB求解
丢番图(Diophantine)方程MATLAB求解 丢番图 (Diophantine) 方程在多项式中的一般形式为 A(z−1)X(z−1)+B(z−1)Y(z−1)=C(z−1)A(z^{-1 ...
- JZOJ 3158 【JSOI2013】丢番图
丢番图 题目大意 1x 1 \over x+ 1y 1 \over y= 1n 1 \over n 求上述方程的解的个数,n为常数. 数据范围 对于 100 100%的数据, n n≤101410^{ ...
- 亿些模板【数论数学】
文章目录 前言 数论数学模板 GCD exgcd 快速幂 线性推逆元 线性推逆元(非连续) 逆元求组合数 矩阵乘法 线性筛素数-埃氏筛 线性筛素数-线性筛 线性筛欧拉-埃氏筛 线性求欧拉 龟速乘 FF ...
- P5253-丢番图【数论】
正题 题目链接:https://www.luogu.org/problem/P5253 题目大意 给一个nnn,求有多少对x,y(x≤y)x,y(x\leq y)x,y(x≤y)使得 1x+1y=1n ...
- 概率论与数理统计思维导图_数学思维到底有多重要?这个学科往往影响国家实力...
原标题:<关于加强数学科学研究工作方案>日前发布--数学思维今何在 密码学家王晓云日前获得了2019年未来科学大奖数学与计算机科学奖.她提出密码哈希函数的碰撞攻击理论,推动帮助新一代密码哈 ...
最新文章
- 图解ecshop之批量上传与批量处理
- KERNELRELEASE
- 141. Linked List Cycle 环形链表
- 从0实现三层神经网络
- 软硬件协同漏洞防护技术【会议】
- 搭建redis主从复制,遇到的问题总结
- 提交健康信息服务器拥挤,python hpilo 监控 hp ilo 服务器健康信息
- CF855G. Harry Vs Voldemort(边双,并查集,dp)
- AI 仿人类人工智能(超级智能)的本质
- oracle的路由器,如何充分使用Oracle?
- element ui tree树节点数据平面化
- unity wsad 鼠标_Unity键盘WASD实现物体移动|chu
- px(像素)与 dp, sp换算公式
- mockito的用法
- 2018 年,做平凡世界里的超人
- 提取EXCEL中事故里程
- 555定时器及其应用
- XMind思维导图安卓PC直装入门下载及常见问题解决教程
- 欧拉编程c语言作业数值分析,数值分析作业 欧拉 龙格库塔
- 计算机安装了打印驱动无法打印,打印机驱动安装后不能打印怎么办解决教程