博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
[POI2011]SEJ-Strongbox
阅读量:6343 次
发布时间:2019-06-22

本文共 1955 字,大约阅读时间需要 6 分钟。

题目大意:

一个有密码箱,数字是0~n-1,其中有若干个密码,密码的特点:若x是密码,y是密码,(x可以等于y)则(x+y)%n也是密码。

给一个n(<=10^14),一个k(k<=min(250000,n)),给k个数(a[k]<n),前k-1个数不是密码,第k个数是密码。

求在0~n-1中,最多有多少个数字是密码?

 

题解:

看起来和数论有一些关系。

而且一定是一个性质题。

结论1:若x是密码,则gcd(n,x)是密码

发现,x是密码,则k*x%n都是密码。

所以,一定存在一个t,c,使得t*x-n*c=gcd(n,x)

并且根据裴属定理,不能用x凑出一个更小的密码比gcd(n,x)更小,

 

结论2:若x,y是密码,则gcd(x,y)是密码。

根据裴属定理,p*x+q*y=gcd(x,y)有整数解。

如果q是负数q=-q,那么就是p*x+(c*n-q)*y=gcd(x,y)+c*n*y

那么,就存在非负数p,q使得p*x+q*y=gcd(x,y) mod n

 

结论3:若x是所有密码中最小的那一个,那么,所有的密码就是x,2x,3x,...kx,并且x是n的约数。

反证。设x是最小的,y是另一个密码,若x不是y的约数,那么gcd(x,y)<x,根据结论二,那么gcd(x,y)就是一个更小的密码。矛盾。

所以,任意的y都是x的倍数。

由于对于一个密码z,根据结论1,gcd(n,z)也是密码,所以,最小的密码x是gcd(n,z)的约数,也就是n的约数。

 

所以,如果我们求出了满足条件的x,那么n/x就是答案。

我们密码数量最多,所以,x必须取最小的。

 

由于给了一个a[k]是密码,而x又是n的约数,所以x就一定是gcd(a[k],n)的约数。

并且,x不能是a[1~k-1]的约数,只要是,那么x就能凑出ai,与ai不是密码矛盾。

只要不存在这样的ai,那么x一定可以是最小的密码(裴属定理可以证明)。

 

所以,我们可以枚举gcd(a[k],n)的约数,然后排序。

从小到大枚举x,再暴力验证是否是a[i]的约数,第一个符合的x,n/x就是答案。

O(ksqrt(n))

 

代码:

 

#include
using namespace std;typedef long long ll;const int N=250000+5;ll n,k;ll gcd(ll a,ll b){ return b?gcd(b,a%b):a;}ll a[N],fac[N];int tot;int main(){ scanf("%lld%lld",&n,&k); for(int i=1;i<=k;i++)scanf("%lld",&a[i]); ll g=gcd(a[k],n); //cout<<"gg "<
<

 然而这个代码BZOJ过不去

还要优化。

发现,每次对一个因数判断所有的k个数代价太大。

一个a[i]影响的是哪些gcd(a[k],n)的因数呢?

一定是gcd(a[i],a[k])的因数。

所以,我们可以把所有的gcd(a[i],a[k])的因数找出来干掉。

然后直接再扫一遍,找到最小的没有被干掉的因数即可。

 

复杂度?O(sqrt(n)+klogn+玄学)

代码:

#pragma GCC optimize(2)#pragma GCC optimize(3)#include
#define ri register int #define numb (ch^'0')using namespace std;typedef long long ll;const int N=250000+5;void rd(ll &x){ x=0; char ch; while(!isdigit(ch=getchar())); for(x=numb;isdigit(ch=getchar());x=(x<<1)+(x<<3)+numb);}ll n,k;ll gcd(ll a,ll b){ return b?gcd(b,a%b):a;}ll a[N],fac[N];bool kil[N];int tot;int main(){ scanf("%lld%lld",&n,&k); for(ri i=1;i<=k;i++)scanf("%lld",&a[i]); ll g=gcd(a[k],n); for(ri i=1;i

 

转载于:https://www.cnblogs.com/Miracevin/p/9777115.html

你可能感兴趣的文章
TF-IDF与余弦相似性的应用(一):自动提取关键词
查看>>
javascript面向对象2
查看>>
限制容器对CPU的使用 - 每天5分钟玩转 Docker 容器技术(28)
查看>>
jquery 实现的一个 随机云标签网页背景
查看>>
RPC
查看>>
android广播事件处理broadcast receive
查看>>
在eclipse 里面 修改tomcat的配置--Server Locations
查看>>
网站 mvc url 路径 设置 为 *.html 的原因
查看>>
mybatis 开启使用 默认的 二级缓存
查看>>
docker 容器 创建和 使用
查看>>
SQLITE使用指南
查看>>
用Maven部署war包到远程Tomcat服务器
查看>>
android字体大小的设置
查看>>
2015.06.04 工作任务与心得
查看>>
icinga2使用587端口发邮件
查看>>
hpasmcli查看HP服务器内存状态
查看>>
极客工具
查看>>
【14】Python100例基础练习(1)
查看>>
boost bind使用指南
查看>>
oracle高速缓存机制
查看>>