【CF819D】Mister B and Astronomers EXGCD
【CF819D】Mister B and Astronomers
题意:小鼠Jack想当太空人(哦不,太空鼠)!为此,它在夜晚带领一堆小朋友一起来到户外看星星。一共有 $n$ 只小鼠,这些小鼠围成一圈轮流观察夜空。具体地,第 $i$ 只小鼠会在第 $(i-1)\%n$ 只小鼠观察夜空之后的第 $a_i$ 秒,抬头观察 $1$ 秒钟的夜空。即: $1$ 号小鼠在第 $0$ 秒观察夜空,$2$ 号小鼠在第 $a_2$ 秒观察夜空,$3$ 号小鼠在第 $a_2+a_3$ 秒观察夜空 ... $1$ 号小鼠在第 $a_2+a_3+...+a_n+a_1$ 秒观察夜空。今晚小鼠们特别想观察到的是名为 $X$ 的星星,但是他们并不知道 $X$ 会在什么时候闪烁,只知道 $X$ 会每隔 $T$ 秒闪烁一次,每次闪烁一秒钟,且第一次闪烁的时间在 $0$ 到 $T-1$ 之间。如果一个小鼠在观察夜空时,$X$ 恰好在闪烁,则这个小鼠就会在这一秒发现 $X$ 闪烁。我们定义第 $i$ 只小鼠的幸运值为:有多少 $t\in [0,T),t\in Z$ ,满足 如果 $X$ 第一次闪烁的时间是第 $t$ 秒,那么 $i$ 将会是第一个发现 $X$ 闪烁的小鼠。现在Jack想知道每个小鼠的幸运值是多少。
题解:设$S=\sum a_i$,那么假如第i只小鼠在第k*S+x秒钟观察到了星星,而这个时刻的星星已经被观察过了,就是说明存在一只小鼠在第k'*S+y秒钟观察到了星星,其中k'<k或k'=k且y<x,对于后一种情况我们特判掉。对于前一种情况,我们得到同余式(k-k')S=y-x(mod T),我们想知道的是k-k'的最小值。
现在我们只需要解这个同余式即可,先将S和T都除以gcd(S,T),同时将所有x按%gcd分组,显然只有同一组内的才会产生贡献。然后我们用exgcd解出k*S=-x(mod T)的k,在set里找一下k的前驱就能得到最小的k-k'了。
#include <cstdio>
#include <cstring>
#include <iostream>
#include <set>
#include <algorithm>
#include <vector>
using namespace std;
const int maxn=200010;
typedef long long ll;
int n,m;
ll S,T,d;
int ban[maxn],p[maxn];
ll t[maxn],k[maxn],top[maxn],ans[maxn];
vector<int> v[maxn];
vector<int>::iterator vi;
set<ll> s;
set<ll>::iterator si;
ll exgcd(ll a,ll b,ll &x,ll &y)
{
if(!b)
{
x=1,y=0;
return a;
}
ll ret=exgcd(b,a%b,x,y),t=x;
x=y,y=t-a/b*x;
return ret;
}
ll gcd(ll a,ll b)
{
return !b?a:gcd(b,a%b);
}
inline bool cmp(const int &a,const int &b)
{
return (t[a]%d==t[b]%d)?(a<b):(t[a]%d<t[b]%d);
}
inline int rd()
{
int ret=0,f=1; char gc=getchar();
while(gc<'0'||gc>'9') {if(gc=='-') f=-f; gc=getchar();}
while(gc>='0'&&gc<='9') ret=ret*10+(gc^'0'),gc=getchar();
return ret*f;
}
int main()
{
T=rd(),n=rd();
int i;
S=rd();
for(i=2;i<=n;i++) t[i]=rd(),S+=t[i],t[i]+=t[i-1];
ll x,y,b;
d=exgcd(S,T,x,y);
b=T/d,x=(x%b+b)%b;
for(i=1;i<=n;i++)
{
if(s.find(t[i]%T)!=s.end()) ban[i]=1;
else s.insert(t[i]%T),p[++p[0]]=i;
}
if(S%T==0)
{
for(i=1;i<=n;i++) printf("%d ",1-ban[i]);
return 0;
}
sort(p+1,p+p[0]+1,cmp);
for(i=1;i<=p[0];i++)
{
if(i==1||t[p[i]]%d!=t[p[i-1]]%d) m++;
v[m].push_back(p[i]);
}
for(i=1;i<=m;i++)
{
s.clear();
for(vi=v[i].begin();vi!=v[i].end();vi++)
{
k[*vi]=x*(t[*vi]/d%b)%b;
if(s.find(k[*vi])!=s.end()) ban[*vi]=1;
else s.insert(k[*vi]);
}
for(vi=v[i].begin();vi!=v[i].end();vi++)
{
si=s.upper_bound(k[*vi]);
if(si==s.end()) si=s.begin(),ans[*vi]=(*si)+b-k[*vi];
else ans[*vi]=(*si)-k[*vi];
}
}
for(i=1;i<=n;i++) printf("%lld ",ans[i]);
return 0;
}
【CF819D】Mister B and Astronomers EXGCD的更多相关文章
- 【CF819C】Mister B and Beacons on Field 数学
[CF819C]Mister B and Beacons on Field 题意:外星人盯上了Farmer Jack的农场!我们假设FJ的农场是一个二维直角坐标系,FJ的家在原点.外星人向FJ的农场上 ...
- 【BZOJ3122】[Sdoi2013]随机数生成器 BSGS+exgcd+特判
[BZOJ3122][Sdoi2013]随机数生成器 Description Input 输入含有多组数据,第一行一个正整数T,表示这个测试点内的数据组数. 接下来T行,每行有五个整数p,a,b, ...
- 【BZOJ 1319】 Sgu261Discrete Rootsv (原根+BSGS+EXGCD)
1319: Sgu261Discrete Roots Time Limit: 1 Sec Memory Limit: 64 MBSubmit: 389 Solved: 172 Descriptio ...
- 【BZOJ5418】【NOI2018】屠龙勇士(数论,exgcd)
[NOI2018]屠龙勇士(数论,exgcd) 题面 洛谷 题解 考场上半个小时就会做了,一个小时就写完了.. 然后发现没过样例,结果大力调发现中间值爆\(longlong\)了,然后就没管了.. 然 ...
- 【bzoj2242】[SDOI2011]计算器 EXgcd+BSGS
题目描述 你被要求设计一个计算器完成以下三项任务: 1.给定y,z,p,计算Y^Z Mod P 的值: 2.给定y,z,p,计算满足xy≡ Z ( mod P )的最小非负整数: 3.给定y,z,p, ...
- 【Codeforces Round #421 (Div. 2) B】Mister B and Angle in Polygon
[题目链接]:http://codeforces.com/contest/820/problem/B [题意] 给你一个正n边形; 然后让你在这正n边行中选3个点,组成一个角; 找出角的大小和所给的角 ...
- 【Codeforces Round #421 (Div. 2) A】Mister B and Book Reading
[题目链接]:http://codeforces.com/contest/820/problem/A [题意] 每天看书能看v页; 且这个v每天能增加a; 但是v有上限v1; 然后每天还必须往回看t页 ...
- 【bzoj1004】[HNOI2008]Cards
1004: [HNOI2008]Cards Time Limit: 10 Sec Memory Limit: 162 MBSubmit: 2928 Solved: 1754[Submit][Sta ...
- 【BZOJ3456】【CDQ分治+FNT】城市规划
试题来源 2013中国国家集训队第二次作业 问题描述 刚刚解决完电力网络的问题, 阿狸又被领导的任务给难住了. 刚才说过, 阿狸的国家有n个城市, 现在国家需要在某些城市对之间建立一些贸易路线, 使得 ...
随机推荐
- mysql 找不到或无法加载已注册的 .Net Framework Data Provider和Unable to find the requested .Net Framework Data Provider. It may not be installed解决
需要安装 mysql-connector-net-6.7.4.msi 在C盘安装mysql的位置找到三个DLL,复制到Bin文件夹下 在Web.config文件中添加对应配置: <system. ...
- mysql 类型
1.bigint 范围(-2^63 (-9223372036854775808) 到 2^63-1 (9223372036854775807) ) 字节大小(8个字节) 2.int 范围(-2^31 ...
- How to make an HTTP request in Swift
from: http://stackoverflow.com/questions/24016142/how-to-make-an-http-request-in-swift You can use N ...
- MongoDB 之 幽灵操作避免
进行静态加载数据到集合的过程中可能会出现. 假设建立一个任务(Job):在MongoDB中进行千条更新操作,开始后迅速终止任务,终止所有更新操作,但依然发现新的更新任务在不断出现,即使任务已经停止. ...
- MongoDB的php可视化管理工具
使用MongoDB命令查看很不方便 于是想把爬来的数据导出来,看爬来的数据是否正确 打开cmd,执行 mongoexport -d test -c blogs --csv -f title,link, ...
- Docker技术原理
Docker技术原理之Linux Namespace(容器隔离) https://blog.csdn.net/songcf_faith/article/details/82748987 Docker技 ...
- ajax+json+Struts2实现list传递(转载)
一.首先需要下载JSON依赖的jar包.它主要是依赖如下: json-lib-2.2.2-jdk15 ezmorph-1.0.4 commons-logging-1.0.4 c ...
- Android学习——在Android中使用OpenCV的第一个程序
刚開始学习Android,因为之前比較熟悉OpenCV,于是就想先在Android上执行OpenCV试试 =============================================== ...
- Dubbo -- 系统学习 笔记 -- 成熟度
Dubbo -- 系统学习 笔记 -- 目录 成熟度 功能成熟度 策略成熟度 啦啦啦
- 【Python】Linux Acanoda PySpark Spark
1.安装 Acanoda 2.安装 Spark和Scala 3.安装 PySpark 4.将Spark的Python目录拷贝至 Acanoda目录下 5.安装py4j,切换anaconda中bin目 ...