[BZOJ 2277]Poi2011 Strongbox

蒟蒻被神题虐暴了。。。。。

感觉之前的数论都白学了

考虑可行的密码x，y，那么(x+y)%n也可行。那么也就是说k*gcd(x,y)%n都可行，我们权且称gcd(x.y)为基数p。

我们已经有了一个可行的密码t，那么p就是t的约数。

再者我们有m-1个不可行的密码x，也就是说p不能满足k*p%n=x，这个可以用扩欧判断一下。

最后我们还需要看到t必须是n的约数，因为如果不是这样的话，必定会有0<k*p%n<p,那么p应为gcd(k*p%n,p)

这样我们得到一个算法：从小到大枚举所有约数，判断所有不可行的密码。找到最小的p，答案为n/p，复杂度为O(sqrt(n)m)

超时。。。。。。

表算似乎很高级，于是就在暴力上改了改：对于每个不可行的密码，判断时有用的只是他与n的约数，取其gcd，这样相同的数太多了，去重排序，会少很多

然后水过去了。。。。。

MAIN上的速度（国外评测及好慢。。。。）

附上代码：

#include <cstdio>
#include <cstring>
#include <algorithm>

using namespace std;
typedef long long LL;
int m,N,i,j;
LL n,a[250005],b[250005],p[3000005],ans;

LL gcd(LL a,LL b){
  while (b) b^=a^=b^=a%=b;
  return a;
}

bool Judge(LL x){
  for (int i=m-1;i>0;i--){
  	if (b[i]<x) return 1;
    if (b[i]%x==0) return 0;
  }
  return 1;
}

int main(){
  //freopen("sej.in","r",stdin);
  //freopen("sej.out","w",stdout);
  scanf("%lld%d",&n,&m);
  for (i=1;i<=m;i++){
  	scanf("%lld",&a[i]);
  	if (a[i]==0) a[i]=n;
  }
  for (i=1;i<=a[m]/i;i++)
  if (a[m]%i==0){
  	if (n%i==0) p[++N]=i;
  	if (a[m]/i!=i && n%(a[m]/i)==0)
  	  p[++N]=a[m]/i;
  }
  sort(p+1,p+N+1);
  for (i=1;i<m;i++)
    b[i]=gcd(n,a[i]);
  sort(b+1,b+m);
  m=unique(b+1,b+m)-b;
  for (i=1;i<=N;i++)
    if (Judge(p[i]))
	  {ans=n/p[i];break;}
  printf("%lld\n",ans);
  return 0;
}