五校联考 running (欧拉函数)
题面
$ solution: $
讲真吧,这道题真的出得,嗯,太恐怖了。考场上这道题真的把我看懵了,这道题以前是见过的,但欧拉函数?我学过吗?一道容斥都要超时的题目,我都要为我自己点根香了,拿着gcd一顿乱搞,果然搞出了个0分。不得不承认博主的数学真的太渣了,但这道题的解题思路真的太妙了(因为渣所以必须学习!)。
首先,一个必须要知道的东西,操场是环形的(即 \((mod~n)\) 意义下的)。若第$ k $个格子可以被第 $ i $ 个同学踩到,那么必定存在一个 $ x $ 使 $ x\times a_i $ 在 \((mod~n)\)之后等于 $ k $ !而这也就相当于:存在一对 $ x $ 和 $ y $ 使得 $ x\times a_i-y*n=k $
而这也就牵扯到我们整除的一些性质了,要满足上面这个式子有整数解,充要条件是 $ (a_i,n)|k $ 也就是说所有能被 $ (a_i,n) $ 整除的格子都可以被第 $ i $ 位同学踩到!!!可是这样的话,有可能一个格子会被多个同学踩到,我们怎么去重呢?
我们发现产生重复的原因是我们的充要条件即 $ (a_i,n)|k (i\in m)$ ,有可能存在多个 $ i $ 能使条件成立,并存在多个 $ i $ 使它们的 $ (a_i,n) $ 存在公因数,然后我们考虑如何让重复的情况只出现一次:就是枚举\(n\)的约数(把 $ (a_i,n) $ 所有可能的值枚举一遍),枚举并判断它是否有贡献之后我们不算存在多少 $ k $ 可以整除它,我们只算有多少 $ (k,n) $ 能等于它。(因为我们枚举的 $ n $ 的约数会包括它的倍数,当下我们只算好它的贡献!)
我们把枚举的约数定为 $ d $ ,那我们就只要求有多少个 $ k $ ,使 $ (k,n)=d $ ,也就是求有多少个 $ k $ ,使 $ (\frac{k}{d} ,\frac{n}{d})=1 $ 也就是求欧拉函数了。
$ code: $
#include<iostream>
#include<cstdio>
#include<iomanip>
#include<algorithm>
#include<cstring>
#include<cstdlib>
#include<ctime>
#include<cmath>
#include<vector>
#include<queue>
#include<map>
#include<set>
#define ll long long
#define db double
#define inf 0x7fffffff
#define rg register int
using namespace std;
int n,m,top,ans;
int a[55];
int b[100005];
inline int qr(){
char ch;
while((ch=getchar())<'0'||ch>'9');
int res=ch^48;
while((ch=getchar())>='0'&&ch<='9')
res=res*10+(ch^48);
return res;
}
inline int gcd(int x,int y){
if(!y)return x;
return gcd(y,x%y);
}
inline int phi(int x){//单独求欧拉函数
int y=x;
for(rg i=2;i<=top&&x>1;++i){
if(!(x%b[i]))y=y/b[i]*(b[i]-1);
while(!(x%b[i]))x/=b[i];
}return y;
}
int main(){
//freopen("running.in","r",stdin);
//freopen("running.out","w",stdout);
ans=n=qr(),m=qr();
for(rg i=1;i<=m;++i)a[i]=qr();
for(rg i=1;i<=sqrt(n);++i)
if(n%i==0)b[++top]=i,b[++top]=n/i;//枚举所有约数
sort(b+1,b+top+1);
for(rg i,k=1;k<=top;++k){
i=b[k];
for(rg j=1;j<=m;++j)
if(i%gcd(a[j],n)==0)//如果这个约数有贡献(即存在一个人使得(ai,n)|k)
{ans-=phi(n/i);break;}
} printf("%d\n",ans);
return 0;
}
五校联考 running (欧拉函数)的更多相关文章
- 【五校联考1day2】JZOJ2020年8月12日提高组T2 我想大声告诉你
[五校联考1day2]JZOJ2020年8月12日提高组T2 我想大声告诉你 题目 Description 因为小Y 是知名的白富美,所以自然也有很多的追求者,这一天这些追求者打算进行一次游戏来踢出一 ...
- 【五校联考1day2】JZOJ2020年8月12日提高组T1 对你的爱深不见底
[五校联考1day2]JZOJ2020年8月12日提高组T1 对你的爱深不见底 题目 Description 出乎意料的是,幸运E 的小R 居然赢了那个游戏.现在欣喜万分的小R 想要写一张明信片给小Y ...
- 【JZOJ4813】【NOIP2016提高A组五校联考2】running
题目描述 小胡同学是个热爱运动的好孩子. 每天晚上,小胡都会去操场上跑步,学校的操场可以看成一个由n 个格子排成的一个环形,格子按照顺时针顺序从0 到n-1 标号. 小胡观察到有m 个同学在跑步,最开 ...
- 【NOIP2016提高A组五校联考2】running
题目 小胡同学是个热爱运动的好孩子. 每天晚上,小胡都会去操场上跑步,学校的操场可以看成一个由n个格子排成的一个环形,格子按照顺时针顺序从0 到n- 1 标号. 小胡观察到有m 个同学在跑步,最开始每 ...
- 【2020五校联考NOIP #4】今天的你依旧闪耀
题面传送门 题意: 对于一个长度为 \(n\)(\(n\) 为偶数)的排列 \(p\),定义一次"变换"后得到的排列 \(p'\) 为: \(p'_i=\begin{cases}p ...
- 【五校联考3day2】C
題意: 現有一平面直角坐標系,有n個點,每一個點必須向某一個方向發射射線,且任意一條射線必須與某一條坐標軸平行.定義一種發射射線的方案是合法的,則方案必須滿足: 1.沒有一條射線交叉 2.沒有一條射線 ...
- 五校联考R1 Day1T3 平面图planar(递推 矩阵快速幂)
题目链接 我们可以把棱柱拆成有\(n\)条高的矩形,尝试递推. 在计算的过程中,第\(i\)列(\(i\neq n\))只与\(i-1\)列有关,称\(i-1\)列的上面/下面为左上/左下,第\(i\ ...
- 五校联考R1 Day2T2 矩阵matrix(容斥)
题目链接 容易想到容斥,但是很恶心,因为要对行和列都容斥,然后行+列又要容斥.. 于是得到\(O(nm\log)\)的做法. 就有70分了: #include <cstdio> #incl ...
- 五校联考模拟赛Day2T2矩阵(容斥原理)
题意 $n * m$的网格,对其进行黑白染色,问每一行每一列至少有一个黑格子的方案数. Sol 考场上只会$n^3$的dp,还和指数级枚举一个分qwq 设$f[i][j]$表示到了第$i$行,已经有$ ...
随机推荐
- lambda函数对象与作用域 (重点)
Python函数在定义的时候,默认参数L的值就被计算出来了,即[],因为默认参数L也是一个变量,它指向对象[],每次调用该函数,如果改变了L的内容,则下次调用时,默认参数的内容就变了,不再是函数定义时 ...
- BZOJ1041 HAOI2008圆上的整点(数论)
求x2+y2=r2的整数解个数,显然要化化式子.考虑求正整数解. y2=r2-x2→y2=(r-x)(r+x)→(r-x)(r+x)为完全平方数→(r-x)(r+x)/d2为完全平方数,d=gcd(r ...
- Debug时含有的子元素,在代码里获取不到的问题
比如,Debug时如下图展示: 我想要获取的是:ansList.get(i).getComponent().getConnectorId() debug时明明有这个元素,但是当我写出来的时候却发现:a ...
- MT【240】6*6放黑白子
$6*6$的方格中放三个完全相同的黑子和三个完全相同的白子,要求每行每列都有一个棋子,且每一格只有一个棋子.问有多少不同放法? 解:$\dfrac{36*25*16*9*4*1}{3!*3!}=144 ...
- MT【15】证明无理数(1)
证明:$tan3^0$是无理数. 分析:证明无理数的题目一般用反证法,最经典的就是$\sqrt{2}$是无理数的证明. 这里假设$tan3^0$是有理数,利用二倍角公式容易得到$tan6^0,tan1 ...
- Leetcode 66.加一 By Python
思路 如果单独操作最后一个元素,令其加一,满十进一,会挺麻烦的,要分情况. 所以我的思路是将list还原到字符串,再变成数值,直接+1,然后再还原到list.详见代码 代码 class Solutio ...
- if语句实例
if 单分支#!/bin/bashRATE=`df -hT | grep "/boot" | awk '{print $6}' | cut -d "%" -f1 ...
- 网络调试 adb connect
$su #stop adbd #setprop service.adb.tcp.port 5555 #start adbd 然后在电脑端打开cmd,输入adb命令: #adb connect < ...
- P1274 魔术数字游戏 naive搜索+剪枝
真的naive...... 我把所有能剪的枝都剪了才过的.否则就是TTT 还有个很神奇的事:数组作为参数传进递归函数时会造成上一层函数里的数组的改变.这个我TM调了一天. 下面奉上代码 #includ ...
- A1059. Prime Factors
Given any positive integer N, you are supposed to find all of its prime factors, and write them in t ...