【JZOJ 4699】 Password

Description

Analysis

这题比赛时的思路是对的，但是少考虑了一些情况，估了100得了0 其实可以统计一下A数组的每个数出现了多少次，按数的大小排好序。 你可以计算出当前这个数的出现次数，用n^2的时间再计算当前数对后面的数的影响。 至于计算当前的数，解个方程就好了。 具体细节看代码。

Code

#include<cstdio> #include<cmath> #include<algorithm> #define fo(i,a,b) for(int i=a;i<=b;i++) using namespace std; typedef long long ll; const int N=1010,M=100010; int n,m,g0; ll a[N*N],b[N],ans[N],bz[M]; struct lyd { ll x,y,z; }c[M]; bool cmp(int a,int b){return a>b;} int gcd(int x,int y) { if(x%y==0) return y; else return gcd(y,x%y); } int find(int x) { int l=1,r=m,mid; while(l<r) { mid=(l+r)>>1; if(x<c[mid].x) l=mid+1; else r=mid; } return l; } int main() { scanf("%d",&n); fo(i,1,n*n) scanf("%lld",&a[i]); sort(a+1,a+n*n+1,cmp); fo(i,1,n*n) if(a[i]!=a[i-1]) c[++m].x=a[i],c[m].y=1; else c[m].y++; fo(i,1,m) { ll B=2*bz[i],C=-c[i].y; ll times=(-B+sqrt(B*B-4*C))/2; if(times<=0) continue; c[i].z=times; fo(j,1,i-1) { if(c[j].z<=0) continue; int t=gcd(c[i].x,c[j].x); int pos=find(t); bz[pos]+=c[i].z*c[j].z; } fo(j,i+1,m) if(c[i].x%c[j].x==0) bz[j]+=c[i].z; fo(j,1,times) ans[++ans[0]]=c[i].x; } sort(ans+1,ans+ans[0]+1,cmp); fo(i,1,ans[0]) printf("%lld ",ans[i]); return 0; }