题目大意:一个集合A,包含了0~n-1这n个数。另有一个集合B,满足:

      1.B是A的子集。

      2.如果a、b均在B中,则(a+b)%n也在B中(a=b亦可)

    给出k个数ai,前k-1个不在B中,第k个在B中,问B最大有几个元素。

数据范围:1<=k<=250000,k<=n<=1e14;

(以下x,y均代表互不相关的整数)

考虑集合B。

假设B中有一元素v。则方程 v*x-n*y=c(即v*x与c对n同余)当且仅当gcd(v,n)|c 时有正整数解x,y。

那么显然c=gcd(v,n)时方程有正整数解x,y。则v∈B->gcd(v,n)∈B……①

又因为方程 v*x+u*y=gcd(u,v)必有整数解x、y,所以根据B的特点2可得 u,v∈B->gcd(u,v)∈B……②

根据题意,(gcd(v,n)*x)%n也必在B中……③

由①、②,设B中所有数的最大公约数与n的最大公约数为g,则g必在B内且为B中最小正整数,又根据③,g,2*g,3*g,4*g……均在B中。有因为g是所有元素的最大公约数,所以B中不会再有别的元素。

以上,我们只要求出满足条件的最小g即可。g的限制条件如下:ai%g!=0,gcd(n,ak)%g==0。

最自然的想法是O(sqrt(gcd(n,ak)))枚举g再对于每个g进行O(k)检验,显然会T。然后……我就不会做了……

根据机智的题解,到这一步我们有以下两种做法:

一、

对于1<=n<=1e14,n最多有17280个因子,那么我们发现g最多有17280个。将所有gcd(n,ai)去重后,个数也会降到17280以下。这样效率就基本OK了。

 #include <bits/stdc++.h>
#define ll long long
using namespace std;
const int sz=25e4+;
ll k,a[sz],n,q;
int l,len;
ll gcd(ll a,ll b){
return b?gcd(b,a%b):a;
}
bool ck(ll v){
for(int i=;i<=l;i++)
if(a[i]%v==)return ;
return ;
}
int main(){
scanf("%lld%lld",&n,&k);
for(int i=;i<=k;i++){
scanf("%lld",&a[i]);
a[i]=gcd(a[i],n);
}
sort(a+,a+k);
for(int i=;i<k;i++)
if(a[i]!=a[i-])a[++l]=a[i];
for(ll j=;j*j<=a[k];j++)if(a[k]%j==0){
if(ck(j)){
printf("%lld",n/j);
return ;
}
if(ck(a[k]/j))q=a[k]/j;
}
printf("%lld",n/q);
return ;
}

二、

在某数学一本通上看到的,但是它的代码和说明蒟蒻我都没看懂,感兴趣的大佬可以看看这篇,貌似和某数学一本通一模一样,看懂的大佬好心的话来教教我。

BZOJ 2277 Poi2011 Strongbox的更多相关文章

  1. bzoj 2277 [Poi2011]Strongbox 数论

    2277: [Poi2011]Strongbox Time Limit: 60 Sec  Memory Limit: 32 MBSubmit: 527  Solved: 231[Submit][Sta ...

  2. 【BZOJ】2277: [Poi2011]Strongbox

    题意 有一个密码箱,\(0\)到\(n-1\)中的某些整数是它的密码.如果\(a\)和\(b\)都是它的密码,那么\((a+b)%n\)也是它的密码(\(a,b\)可以相等).某人试了\(k\)次密码 ...

  3. bzoj2277 [Poi2011]Strongbox

    2277: [Poi2011]Strongbox Time Limit: 60 Sec  Memory Limit: 32 MBSubmit: 498  Solved: 218[Submit][Sta ...

  4. BZOJ 2530 Poi2011 Party 【枚举】

    BZOJ 2530 Poi2011 Party Description Byteasar intends to throw up a party. Naturally, he would like i ...

  5. [bzoj 2216] [Poi2011] Lightning Conductor

    [bzoj 2216] [Poi2011] Lightning Conductor Description 已知一个长度为n的序列a1,a2,-,an. 对于每个1<=i<=n,找到最小的 ...

  6. [poi2011]bzoj 2277 —— strongbox·[洛谷3518]

    ·问题描述· 有一个密码箱,0到n-1中的某些数是它的密码.且满足:如果a和b都是它的密码,那么(a+b)%n也是它的密码.某人试了k次密码,前k-1次都失败了,最后一次成功. 问:该密码箱最多有多少 ...

  7. BZOJ 2277 strongbox (gcd)

    题意 有一个密码箱,0到n-1中的某些整数是它的密码. 且满足,如果a和b都是它的密码,那么(a+b)%n也是它的密码(a,b可以相等) 某人试了k次密码,前k-1次都失败了,最后一次成功了. 问:该 ...

  8. [BZOJ 2212] [Poi2011] Tree Rotations 【线段树合并】

    题目链接:BZOJ - 2212 题目分析 子树 x 内的逆序对个数为 :x 左子树内的逆序对个数 + x 右子树内的逆序对个数 + 跨越 x 左子树与右子树的逆序对. 左右子树内部的逆序对与是否交换 ...

  9. [BZOJ 2350] [Poi2011] Party 【Special】

    题目链接: BZOJ - 2350 题目分析 因为存在一个 2/3 n 大小的团,所以不在这个团中的点最多 1/3 n 个. 牺牲一些团内的点,每次让一个团内的点与一个不在团内的点抵消删除,最多牺牲 ...

随机推荐

  1. freecodecamp数字转化成罗马数字

    做数字转罗马数字时,用到了贪心算法,很好用,记录一下,有时间系统的学一下 罗马数字的规则: 罗马数字网址 1 5 10 50 100 500 1000 I V X L C D M1 1当一个符号在一个 ...

  2. datatables弹窗报错信息屏蔽方法

    在使用datatables的时候,总是会弹出这样的warning: Error: DataTables warning: table id=data_table- Requested unknown ...

  3. MATLAB绘图功能(2) 二维底层绘图修饰

    文末源代码 部分源代码   % x=0:0.1:2*pi; % y=sin(x); % plot(x,y); % line对象 % h = line([-pi:0.01:pi],sin([-pi:0. ...

  4. 最全面的移动APP测试点

    随着互联网,大数据时代的不断推进,演化.移动开发领域得到普遍普及,APP开发如潮水般涌现.下面我将详细介绍app的测试点: 首先我们先熟悉app测试基本流程: 1.1流程图 1.2测试周期 测试周期可 ...

  5. Juploader 1.0 谷歌(chrome)浏览器中成功上传文件后返回信息异常

    在项目中使用了Juploader 1.0无刷新上传文件的js组件,在IE8以上没有问题,代码如下: function InitialUploadDirectly(OnUploadFunc, butto ...

  6. Container&injection

    容器(Container)就是组件和底层服务细节之间的接口.在web组件.企业级Bean等能够执行之前,它必须被装配为一个JavaEE模块,并部署在容器上. 在JavaEE5时代通过注解的方式注入(i ...

  7. Cassandra集群:一,搭建一个三节点的集群

    环境准备 JDK1.8 http://download.oracle.com/otn/java/jdk/8u171-b11/512cd62ec5174c3487ac17c61aaa89e8/jdk-8 ...

  8. quartz.net的使用

    Quartz.NET简介 Quartz.NET作业调度框架是伟大组织OpenSymphony开发的quartz scheduler项目的.net延伸移植版本. Quartz.NET是一个开源的作业调度 ...

  9. 为Linux技术学习推荐看的书籍—《Linux就该这么学》

    成长,是一种经历:经历,是一种人生的体验.人生的意义不在于我们拥有了什么,而在于从中我们体悟了什么.在这短短的三年,却在我的人生中弥足珍贵,在脑海中记忆犹新,在这大学三年里,我从一个莽撞少年成长为一名 ...

  10. hdu 6069 Counting Divisors(求因子的个数)

    Counting Divisors Time Limit: 10000/5000 MS (Java/Others)    Memory Limit: 524288/524288 K (Java/Oth ...