链接:http://acm.hdu.edu.cn/showproblem.php?pid=5514

题意:给定n,m,有n只青蛙,m块编号0~m-1的围成环石头。初始时所有青蛙都在第0块石头上,然后给定n个整数表示n只青蛙每次跳的固定长度,青蛙一直绕环跳,求最后被跳过的石头的编号之和,同一块石头不多次计算。

分析:一看到这个题的时候就想到容斥,以为很容易看到会有多只青蛙跳到同一块石头上。首先我们能想到每只青蛙跳的长度a与总长度m的关系可以确定青蛙在环上跳的“最小长度”为gcd(a,m)。然后有多个最小公约数我们要怎么计算和去重呢?如果只是按普通的容斥去加加减减的话复杂度是O(2^n),并不可取,然后我们观察可以发现每一个gcd都是m的因子,并且!gcd与gcd之间的最小公倍数(即容斥的地方)要么<=m也会是m的一个因子或者>m就不用容斥啦。那么题目就变得简单啦,我们只要计算每一个m的因子所需要做的贡献,多做了就减少做了就加咯,时间复杂度的话看到别人分析说m的因子个数不会超过logm那么时间复杂度就是O(logm*logm+n*logm)了。

代码:

#include<map>
#include<set>
#include<cmath>
#include<queue>
#include<math.h>
#include<cstdio>
#include<vector>
#include<string>
#include<cstring>
#include<iostream>
#include<algorithm>
#pragma comment(linker, "/STACK:102400000,102400000")
using namespace std;
const int N=100010;
const int MAX=151;
const int MOD=1000000007;
const int MOD1=100000007;
const int MOD2=100000009;
const int INF=1000000000;
const double EPS=0.00000001;
typedef long long ll;
typedef unsigned long long uI64;
int read()
{int 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 gcd(int a,int b) {if (b) return gcd(b,a%b);return a;
}
int b[N],q[N],num[N];
int main()
{int i,j,k,n,m,t,x;ll ans;scanf("%d", &t);for (int ca=1;ca<=t;ca++) {scanf("%d%d", &n, &m);k=0;for (i=1;i*i<=m;i++)if (m%i==0) {b[++k]=i;if (i*i!=m) { b[++k]=m/i; }}sort(b+1,b+k+1);k--;memset(q,0,sizeof(q));for (i=1;i<=n;i++) {scanf("%d", &x);x=gcd(m,x);for (j=1;j<=k;j++)if (b[j]%x==0) q[j]=1;}ans=0;memset(num,0,sizeof(num));for (i=1;i<=k;i++) {if (q[i]!=num[i]) {int w=(m-1)/b[i];ans-=(ll)w*(w+1)/2*b[i]*(q[i]-num[i]);w=q[i]-num[i];for (j=i;j<=k;j++)if (b[j]%b[i]==0) num[j]+=w;}}printf("Case #%d: %lld\n", ca, ans);}return 0;
}/*
10
2 12
9 10
3 60
22 33 66
9 96
81 40 48 32 64 16 96 42 72
10 1000000189
1 2 3 5 7 13 17 19 23 29*/

hdu5514Frogs相关推荐

  1. HDU-5514-Frogs

    ACM模版 描述 题解 可以用容斥解,也可以用欧拉函数解,推导过程真心强大. Ps. 图片来源 ITAK's blog.大佬的博客中也有这个题的容斥解法,很强很强很强. 代码 #include < ...

  2. hdu5514Frogs 容斥

    //n只青蛙,每只青蛙每一步跳的步数为a[i] //从0点开始,在一个坐标为0...m-1的环中跳 //青蛙跳的次数没有限制 //问所有的青蛙中任意一只青蛙能够达到的所有的 //坐标之和 #inclu ...

  3. 信息学竞赛中的数学知识 --- 容斥原理

    C++基础数论-----容斥原理 C++基础数论-----容斥原理_C2020lax的博客-CSDN博客_容斥原理c++ C++数论容斥原理----无关的元素 C++数论容斥原理----无关的元素 - ...

最新文章

  1. Hive时间是String格式截取字串和转换数据类型小贴士
  2. AI 经典书单 | 人工智能学习该读哪些书
  3. 关于域帐户将计算机加入域登陆上限问题
  4. Bootstrap警告框
  5. linux强制停止压缩任务,6.25任务(linux压缩、打包命令:zip、tar)
  6. [Vue.js]实战 -- 电商项目(八)
  7. 时富金融:八年左右后内地房价会下降
  8. 云图说|SAP技术画册“一点通”
  9. 华为手机投屏电脑_华为手机如何实现无线投屏?
  10. 330UF16V 10*7.7片式铝电解电容封装
  11. Vue引入百度地图API,添加点击地图拾取坐标并且标记和地址搜索功能
  12. 连接MyBatis内部SqlSession与业务接口的代理类MapperProxy
  13. PS 导入笔刷和导入字体和导入滤镜
  14. 第 4 代 4 端口 HDMI/DVI 分配器龙迅LT86104SXE
  15. CSS选择器和层叠性
  16. 颓废了很久了,需要振作了
  17. gif、bmp、jpg、jpeg图片有什么区别
  18. IMPL10. formality 常用变量浅析
  19. 计算机桌面图标不可以移动怎么调,电脑桌面图标不能随意拖动,怎么随意拖动桌面图标...
  20. 阻止事件冒泡,阻止默认事件,event.stopPropagation()和event.preventDefault(),return false的区别

热门文章

  1. Redis Guava 布隆过滤器实现和准确率测试
  2. 无线连接服务器 很慢,手机无线WiFi很卡很慢怎么办?wifi信号满格但网速慢解决方法...
  3. 宽依赖和窄依赖_Spark宽依赖和窄依赖深度剖析
  4. C++11之防止类型收窄(列表初始化)
  5. 深度学习笔记(6)BatchNorm批量标准化
  6. 计算机配置主要看哪些东西,买电脑主要看哪些配置
  7. Android FaceBook登录 分享获取HashKey(密钥散列)的简单方法
  8. python求股票收益率_python – 有没有办法自动获取许多股票的一般信息,如市盈率,收益率等等?...
  9. ThreadPoolExecutor的四种抛弃策略
  10. 基于Android平台的虚拟试衣关键技术研究