hdu5514

题意

\(m\) 个石子绕成一圈,编号\([0, m - 1]\)。有 \(n\) 个青蛙从 \(0\) 号石子出发,给出每个青蛙的步长,青蛙无限跑圈。问哪些石子至少被一个青蛙经过,求这些石子的编号之和。

分析

假设某个青蛙的步长 \(x\),则一共会经过 $ \frac{m - 1}{x} + 1$ 个石子(包括 \(0\) )。可以用等差数列求和公式计算贡献。

然后找出 \(m\) 的所有因子并考虑哪些因子可能要计算贡献。

例如步长为 \(2\) 和 \(3\) 的计算后,步长为 \(6\) 的情况被计算了两次(或者说 \(6\) 的倍数的编号被计算了两次),应该减去一次,用一个数组记录某个步长被计算了几次贡献,在统计答案的时候减去即可。

code

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int MAXN = 1e4 + 10;
int a[MAXN];
int fac[MAXN];
int vis[MAXN]; // vis[i]表示fac[i]是否可能要计算贡献,如果为 0 一定不计算贡献
int num[MAXN]; // num[i]表示fac[i]被计算了几次贡献
int main() {
int T, kase = 1;
scanf("%d", &T);
while(T--) {
memset(vis, 0, sizeof vis);
memset(num, 0, sizeof num);
int n, m;
scanf("%d%d", &n, &m);
int fc = 0;
for(int i = 1; i * i <= m; i++) {
if(m % i == 0) {
fac[fc++] = i;
if(i * i != m) fac[fc++] = m / i;
}
}
sort(fac, fac + fc);
for(int i = 0; i < n; i++) {
scanf("%d", &a[i]);
a[i] = __gcd(a[i], m);
for(int j = 0; j < fc; j++) {
if(fac[j] % a[i] == 0) {
vis[j] = 1;
}
}
}
ll ans = 0;
for(int i = 0; i < fc - 1; i++) {
ll k = (m - 1) / fac[i];
ans += (k + 1) * k * fac[i] / 2 * (vis[i] - num[i]);
for(int j = i + 1; j < fc - 1; j++) {
if(fac[j] % fac[i] == 0) {
num[j] += vis[i] - num[i];
}
}
}
printf("Case #%d: %lld\n", kase++, ans);
}
return 0;
}

最新文章

  1. Centos7上安装dnf-plugins-core
  2. rocketmq简单搭建
  3. java并发ThreadLocal
  4. java下properties属性文件操作
  5. L9-3.安装PHP软件包
  6. hibernate事务并发问题(脏读,不可重复读,幻读)
  7. [转载] KAFKA分布式消息系统
  8. Openstack_O版(otaka)部署_镜像服务glance部署
  9. 什么是TensorFlow?
  10. elk 中kafka启动脚本和配置文件
  11. php发送短信验证码
  12. spring boot 入门及示例
  13. Hibernate添加日志--log4j
  14. Bootstrap3基础 list-inline 无序列表横向显示
  15. HTML5 实现获取 gzip 压缩数据,并进行解压,同时解决汉字乱码,相关 pako.js
  16. [z]Libevent使用例子,从简单到复杂
  17. BZOJ4818 序列计数
  18. Neo4j资料 Neo4j教程 Neo4j视频教程 Neo4j 图数据库视频教程
  19. 来自MSDN的RibbonGadgets练习
  20. GitHub和git和repo的使用

热门文章

  1. hdu6103 Kirinriki(trick+字符串)
  2. Xinetd服务的安装与配置
  3. hdu4035 Maze 【期望dp + 数学】
  4. 【BZOJ 2241 打地鼠】
  5. Consumer [分组背包]
  6. [poj 2796]单调栈
  7. HDU 多校对抗赛第二场 1004 Game
  8. idea 导入spring 源码注意的问题
  9. Oracle SQL 疑难解析读书笔记(二、汇总和聚合数据)
  10. 移动端list布局,左边固定,右边自适应