[poi2011]bzoj 2277 —— strongbox·[洛谷3518]
·问题描述·
有一个密码箱,0到n-1中的某些数是它的密码。且满足:如果a和b都是它的密码,那么(a+b)%n也是它的密码。某人试了k次密码,前k-1次都失败了,最后一次成功。
问:该密码箱最多有多少个密码?
·输入格式·
输入第一行两个整数分别表示n,k。
第二行为k个用空格隔开的非负整数,表示每次试的密码。(数据保证存在合法解)
·输出格式·
输出一行一个数,表示结果。
·输入样例·
42 5
28 31 10 38 24
·输出样例·
14
·数据范围·
对于100%的数据:1<=k<=250000,k<=n<=10^14。
Solution:
本题考察数学。由题意可知,若x为密码则(x+x)%n为密码,则p*x%n(0<p<n)也为密码。而对于p*x%n=q,等价于p*x-n*c=q。
由引理:a*x+b*y=c(均为整数),有整数解的充要条件是gcd(a,b)|c。证明很容易:不妨设a=p*gcd(a,b),b=q*gcd(a,b) --> a*x+b*y=(p+q)*gcd(a,b)=c,显然要有整数解,则gcd(a,b)|c。
回到本题的条件:p*x-n*c=q。有解则必定满足gcd(x,n)|q,所以必定有p*x-n*c=gcd(x,n)成立,等价于p*x%n=gcd(x,n),则gcd(x,n)一定为一个密码。类似的,对于不同的密码x和y,存在(p*x+q*y)%n为密码,由引理必定存在p*x+q*y=gcd(x,y),与单个x是密码同理gcd(x,y)一定是密码。
而要使得密码最多,由x是密码则p*x%n(0<p<n)为密码可知,当x为最小时,密码最多有n/x个。
具体实现时,我们先求出a[k]=gcd(a[k],n),再使a[i]=gcd(a[i],a[k]),然后从新的a[k]中删去所有是a[i]因子的因子,最后输出答案就是n除以没被删的最小的因子。
代码:
/*数学一本通上的例题——by 520*/
#include<bits/stdc++.h>
#define il inline
#define ll long long
using namespace std;
ll n,k,tot,a[],q[],p[],cnt=;
il ll gi(){
ll a=;char x=getchar();bool f=;
while((x<''||x>'')&&x!='-')x=getchar();
if(x=='-')x=getchar(),f=;
while(x>=''&&x<='')a=a*+x-,x=getchar();
return f?-a:a;
}
il ll gcd(ll a,ll b){return b?gcd(b,a%b):a;}
int main()
{
freopen("strongbox.in","r",stdin);
freopen("strongbox.out","w",stdout);
n=gi(),k=gi();
for(int i=;i<=k;i++)a[i]=gi();
a[k]=gcd(a[k],n);
for(int i=;i<k;i++)a[i]=gcd(a[k],a[i]);
for(ll i=;i*i<=a[k];i++)
if(a[k]%i==){
q[++tot]=i;
if(i*i!=a[k])q[++tot]=a[k]/i;
}
sort(q+,q+tot+);
for(int i=;i<k;i++)p[lower_bound(q+,q+tot+,a[i])-q]=;
for(int i=;i<=tot;i++)
if(p[i])
for(int j=;j<i;j++)
if(q[i]%q[j]==)p[j]=;
while(p[cnt])cnt++;
cout<<n/q[cnt];
return ;
}
[poi2011]bzoj 2277 —— strongbox·[洛谷3518]的更多相关文章
- [bzoj] 3263 陌上花开 洛谷 P3810 三维偏序|| CDQ分治 && CDQ分治讲解
原题 定义一个点比另一个点大为当且仅当这个点的三个值分别大于等于另一个点的三个值.每比一个点大就为加一等级,求每个等级的点的数量. 显然的三维偏序问题,CDQ的板子题. CDQ分治: CDQ分治是一种 ...
- BZOJ 2277 strongbox (gcd)
题意 有一个密码箱,0到n-1中的某些整数是它的密码. 且满足,如果a和b都是它的密码,那么(a+b)%n也是它的密码(a,b可以相等) 某人试了k次密码,前k-1次都失败了,最后一次成功了. 问:该 ...
- 「洛谷3338」「ZJOI2014」力【FFT】
题目链接 [BZOJ] [洛谷] 题解 首先我们需要对这个式子进行化简,否则对着这么大一坨东西只能暴力... \[F_i=\sum_{j<i} \frac{q_iq_j}{(i-j)^2}-\s ...
- BZOJ2120/洛谷P1903 [国家集训队] 数颜色 [带修改莫队]
BZOJ传送门:洛谷传送门 数颜色 题目描述 墨墨购买了一套N支彩色画笔(其中有些颜色可能相同),摆成一排,你需要回答墨墨的提问.墨墨会向你发布如下指令: 1. Q L R代表询问你从第L支画笔到第R ...
- BZOJ2212或洛谷3521 [POI2011]ROT-Tree Rotations
BZOJ原题链接 洛谷原题链接 线段树合并裸题. 因为交换子树只会对子树内部的逆序对产生影响,所以我们计算交换前的逆序对个数和交换后的个数,取\(\min\)即可. 对每个叶子节点建一棵动态开点线段树 ...
- 洛谷 P3307: bzoj 3202: [SDOI2013] 项链
题目传送门:洛谷P3307.这题在bzoj上是权限题. 题意简述: 这题分为两个部分: ① 有一些珠子,每个珠子可以看成一个无序三元组.三元组要满足三个数都在$1$到$m$之间,并且三个数互质,两个珠 ...
- bzoj 4816: 洛谷 P3704: [SDOI2017]数字表格
洛谷很早以前就写过了,今天交到bzoj发现TLE了. 检查了一下发现自己复杂度是错的. 题目传送门:洛谷P3704. 题意简述: 求 \(\prod_{i=1}^{N}\prod_{j=1}^{M}F ...
- bzoj 1014: 洛谷 P4036: [JSOI2008]火星人
题目传送门:洛谷P4036. 题意简述: 有一个字符串,支持插入字符,修改字符. 每次需要查询两个后缀的LCP长度. 最终字符串长度\(\le 100,\!000\),修改和询问的总个数\(\le 1 ...
- [BZOJ 3039&洛谷P4147]玉蟾宫 题解(单调栈)
[BZOJ 3039&洛谷P4147]玉蟾宫 Description 有一天,小猫rainbow和freda来到了湘西张家界的天门山玉蟾宫,玉蟾宫宫主蓝兔盛情地款待了它们,并赐予它们一片土地. ...
随机推荐
- Linux系统运维基础管理命令总结
1.查看系统负载命令:w.uptime [root@localhost ~]# w :: up days, :, user, load average: 0.00, 0.01, 0.05 USER T ...
- meta标签的常见用法
一.定义和用法 <meta> 标签始终位于 head 元素中.<meta> 元素可提供有关页面的元信息(meta-information),元数据不会显示在页面上,但是对于机器 ...
- python环境通过selenium实现自动化web登陆及终端邀请
自动化主要的就是识别对象,可以在网上搜到各种各样的方法,自行百度.下面仅附上一个简单的例子. 环境搭建参考如下链接: https://www.cnblogs.com/hepeilinnow/p/101 ...
- vue的ui库使用Element UI,纯html页面,不使用webpack那玩意
使用手册访问:https://cloud.tencent.com/developer/doc/1270 第一步:在head添加样式 <link rel="stylesheet" ...
- Spring单元测试集成H2数据库
项目源代码在:Spring-H2测试 H2简介 H2数据库是一种由Java编写的,极小,速度极快,可嵌入式的数据库.非常适合用在单元测试等数据不需要保存的场景下面. 以下时其官网的介绍: {% blo ...
- Batch Normalization: Accelerating Deep Network Training by Reducing Internal Covariate Shift
1. 摘要 训练深层的神经网络非常困难,因为在训练的过程中,随着前面层数参数的改变,每层输入的分布也会随之改变.这需要我们设置较小的学习率并且谨慎地对参数进行初始化,因此训练过程比较缓慢. 作者将这种 ...
- 使用 Mesos 管理虚拟机
摘要 为了满足渲染.基因测序等计算密集型服务的需求,UCloud 推出了“计算工厂”产品,让用户可以快速创建大量的计算资源(虚拟机).该产品的背后,是一套基于 Mesos 的计算资源管理系统.本文简要 ...
- php的大小写敏感问题整理
php的大小写敏感问题整理 今天在开发php的过程中,因为命名大小写的问题导致代码错误,所以从网上整理了php的大小写敏感的一些资料,需要的朋友可以参考下. PHP对大小写敏感问题的处理比较乱,写 ...
- SQLSERVER 根据身份证号码 往出生年月日 赋值
update CREW_SailorInfo set DT_DOB= ( case then , ) then , ) else null end) 注:此问题仅供参考 如有疑问 请加QQ群18153 ...
- Python:文件操作总结1——文件基本操作
一.文件的操作流程 1.打开文件,得到文件句柄并赋值给一个变量 2.通过句柄对文件进行操作 3.关闭文件 二.文件的打开与关闭 A.文件的打开——open函数 语法:open(file[,mode[, ...