题目链接

bzoj1220: [HNOI2002]跳蚤

题解

根据裴蜀定理,不定方程的解为未知数的gcd,所以选取的n个数的gcd为1
那么n - 1个数保证没有公约数为m的约数,枚举质因数容斥
质因数的个数上届是log的啊,我真傻,还想了半天QAq
那啥,bzoj高精,你们去做吧Qwq

代码

#include<cstdio>
#include<algorithm>
#define LL long long
inline LL read() { LL x = 0,f = 1; char c = getchar(); while(c < '0' || c > '9')c = getchar(); while(c <= '9' && c >= '0')x = x * 10 + c - '0',c = getchar(); return x * f;
}
const int maxn = 100007;
LL ans = 0,n,m;
LL pow(LL x,LL res) { LL ret =  1; for(;res;res >>= 1,x *= x) if(res & 1) ret *= x ;return ret;
}
int num = 0,rime[maxn];
void divede(int x) { for(int i = 2;i * i <= x;++ i)  if(x % i == 0) { rime[++ num] = i; while(x % i == 0) x /= i; } if(x != 1) rime[++ num]  = x;
}
void solve(int cnt = 1,LL x = 1,int tp = 1) { if(cnt == num + 1) { ans += (tp & 1) ? pow(m / x,n) : -pow(m / x,n); return ;} solve(cnt + 1,x,tp); solve(cnt + 1,x * rime[cnt],tp ^ 1);
}
int main() { n = read(),m = read(); //ans = 2 * pow(m,n); divede(m); solve(); printf("%lld\n",ans); return 0;
}

转载于:https://www.cnblogs.com/sssy/p/9279144.html

luoguP2231 [HNOI2002]跳蚤相关推荐

  1. P2231 [HNOI2002]跳蚤(裴蜀定理/莫比乌斯反演)

    P2231 [HNOI2002]跳蚤 给定一个长度为n+1的一列数,第n+1位为m,前n位小于m 求解使得他n+1个数的加减可以凑出1的方案数 首先可以凑出1,这显然是裴蜀定理,推一推就发现他要求所有 ...

  2. 洛谷P2231 [HNOI2002]跳蚤 [数论,容斥原理]

    题目传送门 跳蚤 题目描述 Z城市居住着很多只跳蚤.在Z城市周六生活频道有一个娱乐节目.一只跳蚤将被请上一个高空钢丝的正中央.钢丝很长,可以看作是无限长.节目主持人会给该跳蚤发一张卡片.卡片上写有N+ ...

  3. [HNOI2002]跳蚤

    题目描述 Z城市居住着很多只跳蚤.在Z城市周六生活频道有一个娱乐节目.一只跳蚤将被请上一个高空钢丝的正中央.钢丝很长,可以看作是无限长.节目主持人会给该跳蚤发一张卡片.卡片上写有N+1个自然数.其中最 ...

  4. P2231 [HNOI2002]跳蚤

    文章目录 R e s u l t Result Result H y p e r l i n k Hyperlink Hyperlink D e s c r i p t i o n Descripti ...

  5. bzoj1588[HNOI2002]营业额统计

    bzoj1588[HNOI2002]营业额统计 题意: n天,每天得到一个值,要求输出每一天和这天得到的值相差最小的之前天得到的值与这个值的差的和.n不知道,不过O(nlog2n)可写. 题解: 说是 ...

  6. 【链表】BZOJ1588: [HNOI2002]营业额统计

    1588: [HNOI2002]营业额统计 Time Limit: 5 Sec  Memory Limit: 162 MB Submit: 17555  Solved: 7179 [Submit][S ...

  7. POJ 1091 跳蚤

    跳蚤 Time Limit: 1000ms Memory Limit: 10000KB This problem will be judged on PKU. Original ID: 1091 64 ...

  8. [LOJ 6042]「雅礼集训 2017 Day7」跳蚤王国的宰相(树的重心+贪心)

    [LOJ 6042]「雅礼集训 2017 Day7」跳蚤王国的宰相 description solution 一个到所有节点距离和最小的节点 ⇔\Leftrightarrow⇔ 树的重心(满足最重的儿 ...

  9. 毕业设计上线啦!----跳蚤部落与基于Comet的WebIM系统开发

    我不清楚把我的毕业设计的东西放上来之后,毕业论文答辩的时候会不会说我是在网上抄袭的,不过我还是果断的发上来与大家分享了!!呵呵,请大家支持!高手就绕道吧! 现在已经放到公网上,并且开始使用,兼容IE6 ...

最新文章

  1. 深入实践Spring Boot3.1 模型设计
  2. 《深入理解ES6》笔记——块级作用域绑定(1)
  3. 进程间通讯,临界区,互斥
  4. php cgi漏洞,CGI漏洞攻击合集上
  5. php psot传值_三种方法教你如何用PHP模拟post提交数据
  6. 顶尖学府 加州伯克利大学开发高效机器人操纵框架
  7. cad vba 打开文件对话框_Excel VBA和文件夹-1.7通过对话框灵活选定文件夹并打开对应文件...
  8. lisp实心圆点怎么画_实心圆点_实心小圆圈符号怎么打出来?
  9. 巧记单词:odd和even
  10. unity粒子系统笔记
  11. 在持续集成 (CI) 环境中使用 Android 模拟器 | AndroidDevSummit 中文字幕视频
  12. 英语四级XYZ字母开头单词+托福高频单词
  13. 关于openfire
  14. 用php把图片合成视频,图片音乐合成视频 多张图片合成视频|图片合成视频软件...
  15. 简练软考知识点整理-智慧城市参考模型
  16. OSM数据如何下载使用(地图数据篇.11)
  17. 技术博文: 头条,知乎、CSDN 三大平台对比
  18. 太阳软件分享丨分布式ip代理池,可用率超高,赶快试试吧!
  19. SOLIDWORKS中的两个基本概念
  20. jQuery的基本使用方法

热门文章

  1. 关于淫秽×××网站泛滥的思考
  2. 硬件测试工程师系列0--写在最开始
  3. 计算机配色在纺织中的应用,浅析纺织染色电脑配色程序
  4. 批量删除的三种实现方式
  5. 排列组合-9.1 排列组合系列问题
  6. kafka性能测试之Comsumer
  7. 深入Preact源码分析(4.20更新)
  8. Naptime 【POJ - 2228】
  9. 移动考勤满足企业实时管理
  10. matlab仿真分形树,使用MATLAB实现随机分形树模拟