[noip模拟赛]跑跑步
https://www.zybuluo.com/ysner/note/1298652
题面
小胡同学是个热爱运动的好孩子。
每天晚上,小胡都会去操场上跑步,学校的操场可以看成一个由\(n\)个格子排成的一个环形,格子按照顺时针顺序从\(0\)到\(n−1\)标号。
小胡观察到有\(m\)个同学在跑步,最开始每个同学都在起点(即\(0\)号格子),每个同学都有个步长\(a_i\),每跑一步,每个同学都会往顺时针方向前进\(a_i\)个格子。
由于跑道是环形的,如果一个同学站在\(n−1\)这个格子上,如果他前进一个格子,他就会来到\(0\)。
他们就这样在跑道上上不知疲倦地跑呀跑呀。
小胡同学惊奇地发现,似乎有些格子永远不会被同学跑到,他想知道这些永远不会被任何一个同学跑到的格子的数目,你能帮帮他吗?(我们假定所有同学都跑到过\(0\)号格子)。
- \(n\leq10^9,m\leq50\)
解析
因没打表猜结论被踩系列
首先要知道一个结论,即第\(i\)个同学能走到的格子一定是\(k*gcd(a_i,n)\)。(也许需要打表)
然后容斥一下就行了。
但是复杂度不太对劲啊。好像只能过\(60pts\)。
考虑优化。
于是我把\(50\)个数中的倍数去掉,从大到小排序让它尽早乘爆,加了诸多此类优化。似乎续到了\(90pts\)。
然后看看这玩意儿
il void dfs(re int x,re ll tot,re int num)
{
if(tot>n) return;
if(x>m)
{
if(num&1) ans-=(n/tot);else ans+=(n/tot);
return;
}
fp(i,0,1)
if(i) dfs(x+1,lcm(a[x],tot),num+1);else dfs(x+1,tot,num);
}
如果\(lcm(a[x],tot)==tot\),那左右两边的结果不是能一正一负恰好抵消吗?
加上这个优化,就快得不行了。
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
#define ll long long
#define re register
#define il inline
#define fp(i,a,b) for(re int i=a;i<=b;i++)
#define fq(i,a,b) for(re int i=a;i>=b;i--)
using namespace std;
const int N=100;
int n,m;
ll ans,a[N],tot;
bool vis[N];
il ll gi()
{
re ll x=0,t=1;
re char ch=getchar();
while(ch!='-'&&(ch<'0'||ch>'9')) ch=getchar();
if(ch=='-') t=-1,ch=getchar();
while(ch>='0'&&ch<='9') x=x*10+ch-48,ch=getchar();
return x*t;
}
il ll lcm(re ll x,re ll y){return x*y/__gcd(x,y);}
il bool cmp(re ll x,re ll y){return x>y;}
il void dfs(re int x,re ll tot,re int num)
{
if(tot>n) return;
if(x>m)
{
if(num&1) ans-=(n/tot);else ans+=(n/tot);
return;
}
if(tot%a[x]==0) return;
fp(i,0,1)
if(i) dfs(x+1,lcm(a[x],tot),num+1);else dfs(x+1,tot,num);
}
int main()
{
freopen("running.in","r",stdin);
freopen("running.out","w",stdout);
n=gi();m=gi();
fp(i,1,m) a[i]=__gcd(gi(),1ll*n);
sort(a+1,a+1+m,cmp);//然后被xzy造了个第i个数为2^i的数据卡了。。。
dfs(1,1,0);
printf("%lld\n",ans);
fclose(stdin);
fclose(stdout);
return 0;
}
再看看标算。
只要枚举\(n\)的约数,令当前枚举到的数为\(d\),若存在一个\(i\),使得\(gcd(a_i,n)|d\)。说明所有\(gcd(i,n)=d\)的格子都能被到达,答案加上\(\phi(\frac{n}{d})\)即可。
证明:
\]
#include<cstdio>
#include<algorithm>
using namespace std;
int gi(){
int x=0,w=1;char ch=getchar();
while ((ch<'0'||ch>'9')&&ch!='-') ch=getchar();
if (ch=='-') w=0,ch=getchar();
while (ch>='0'&&ch<='9') x=(x<<3)+(x<<1)+ch-'0',ch=getchar();
return w?x:-x;
}
int n,m,a[55],ans;
int phi(int x){
int res=x;
for (int i=2;i*i<=x;++i)
if (x%i==0){
res/=i;res*=i-1;
while (x%i==0) x/=i;
}
if (x>1) res/=x,res*=x-1;
return res;
}
void cal(int x){
for (int i=1;i<=m;++i) if (x%a[i]==0) return;
ans+=phi(n/x);
}
int main(){
freopen("running.in","r",stdin);
freopen("running.out","w",stdout);
n=gi();m=gi();
for (int i=1;i<=m;++i) a[i]=__gcd(n,gi());
for (int i=1;i*i<=n;++i)
if (n%i==0){
cal(i);if (i*i<n) cal(n/i);
}
printf("%d\n",ans);return 0;
}
[noip模拟赛]跑跑步的更多相关文章
- NOIP模拟赛20161022
NOIP模拟赛2016-10-22 题目名 东风谷早苗 西行寺幽幽子 琪露诺 上白泽慧音 源文件 robot.cpp/c/pas spring.cpp/c/pas iceroad.cpp/c/pas ...
- NOIP模拟赛 by hzwer
2015年10月04日NOIP模拟赛 by hzwer (这是小奇=> 小奇挖矿2(mining) [题目背景] 小奇飞船的钻头开启了无限耐久+精准采集模式!这次它要将原矿运到泛光之源的矿 ...
- 大家AK杯 灰天飞雁NOIP模拟赛题解/数据/标程
数据 http://files.cnblogs.com/htfy/data.zip 简要题解 桌球碰撞 纯模拟,注意一开始就在袋口和v=0的情况.v和坐标可以是小数.为保险起见最好用extended/ ...
- NOIP模拟赛-2018.11.5
NOIP模拟赛 好像最近每天都会有模拟赛了.今天从高二逃考试跑到高一机房,然而高一也要考试,这回好像没有拒绝的理由了. 今天的模拟赛好像很有技术含量的感觉. T1:xgy断句. 好诡异的题目,首先给出 ...
- Nescafe #29 NOIP模拟赛
Nescafe #29 NOIP模拟赛 不知道这种题发出来算不算侵权...毕竟有的题在$bz$上是权限题,但是在$vijos$似乎又有原题...如果这算是侵权的话请联系我,我会尽快删除,谢谢~ 今天开 ...
- NOIP模拟赛 6.29
2017-6-29 NOIP模拟赛 Problem 1 机器人(robot.cpp/c/pas) [题目描述] 早苗入手了最新的Gundam模型.最新款自然有着与以往不同的功能,那就是它能够自动行走, ...
- contesthunter暑假NOIP模拟赛第一场题解
contesthunter暑假NOIP模拟赛#1题解: 第一题:杯具大派送 水题.枚举A,B的公约数即可. #include <algorithm> #include <cmath& ...
- 队爷的讲学计划 CH Round #59 - OrzCC杯NOIP模拟赛day1
题目:http://ch.ezoj.tk/contest/CH%20Round%20%2359%20-%20OrzCC杯NOIP模拟赛day1/队爷的讲学计划 题解:刚开始理解题意理解了好半天,然后发 ...
- 队爷的Au Plan CH Round #59 - OrzCC杯NOIP模拟赛day1
题目:http://ch.ezoj.tk/contest/CH%20Round%20%2359%20-%20OrzCC杯NOIP模拟赛day1/队爷的Au%20Plan 题解:看了题之后觉得肯定是DP ...
随机推荐
- linux命令 iperf-网络性能测试工具
博主推荐:更多网络测试相关命令关注 网络测试 收藏linux命令大全 iperf命令是一个网络性能测试工具.iperf可以测试TCP和UDP带宽质量.iperf可以测量最大TCP带宽,具有多种参数和 ...
- 等待某(N)个线程执行完再执行某个线程的几种方法(Thread.join(),CountDownLatch,CyclicBarrier,Semaphore)
1.main线程中先调用threadA.join() ,再调用threadB.join()实现A->B->main线程的执行顺序 调用threadA.join()时,main线程会挂起,等 ...
- codeforces 361A
//这题看着吓人,为何这么水 #include<stdio.h> int main() { int n,m,i,j; while(scanf("%d%d",& ...
- Linux下汇编语言学习笔记13 ---
这是17年暑假学习Linux汇编语言的笔记记录,参考书目为清华大学出版社 Jeff Duntemann著 梁晓辉译<汇编语言基于Linux环境>的书,喜欢看原版书的同学可以看<Ass ...
- iOS点击cell时,控件背景色消失的解决方法
同时调用一下两个方法: - (void)setSelected:(BOOL)selected animated:(BOOL)animated { [super setSelected:selected ...
- ***CodeIgnite/CI 去掉 index.php的 配置
CI有效删除URL中的index.php 参考: http://codeigniter.org.cn/forums/thread-15444-1-1.html 读CI的使用手册的话,关于如何有效删除U ...
- MVC view页面需要多个model,复杂网页的处理
需求描述 一个比较复杂的页面,界面中包含的元素数据来自于许多个有关联或者无关联的表,然后我们要做的就是将数据呈现在界面上. 10年前大概都是这么干的 直接写一个复杂的SQL语句,返回一个包含所需数据的 ...
- Batch update returned unexpected row count from update [0]; actual row count: 0; expected: 1
在操作hibernate数据库时,调用saveOrUpdate方法进行更新保存对象时, (1)ID为null时执行SAVE,但是前端jsp通过<input type="hidden&q ...
- GNS3模拟的硬件
Hardware emulated by GNS3 Cisco 1700 Series 1700s have one or more interfaces on the motherboard, 2 ...
- 015 WAN
Router#config t Enter configuration commands, one per line. End with CNTL/Z. Router(config)#int s0/ ...