BZOJ1978 [BeiJing2010]取数游戏 建图+拓扑序

很容易想到在可以满足L≤gcd(i,j)的i,j之间连边 但是暴力n2枚举肯定过不了啦 考虑枚举每一个大于L的公因数 我们对于每个数根号n枚举他的因子 对每一个因子建一个vector来存 然后每一个vector里的i和i+1连边 yy一下 这样连就保证了L≤gcd(i,j)的点(i,j)可以到达通过几条边到达 然后我们就在图上跑拓扑序啦

#include<bits/stdc++.h>
using namespace std;
const int N=1e6+5;
inline int read(){
int x=0,f=1;char ch=getchar();
while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
return x*f;
}
int n,l,ans,tot,f[50005],head[50005],in[50005];
queue<int>q;
vector<int>b
;
struct Egde{
int v,nxt;
}e[N*2];
void add(int u,int v){
e[++tot].v=v,e[tot].nxt=head[u],head[u]=tot,in[v]++;
}
void topsort(){
for(int i=1;i<=n;i++) if(!in[i]) f[i]=1,q.push(i);
while(!q.empty()){
int x=q.front();q.pop();
ans=max(ans,f[x]);
for(int i=head[x];i;i=e[i].nxt){
int j=e[i].v;
f[j]=max(f[j],f[x]+1);
if(!--in[j]) q.push(j);
}
}
}
int main(){
#ifdef Devil_Gary
freopen("in.txt","r",stdin);
#endif
n=read(),l=read();
for(int p=1,x,i;p<=n;p++){
x=read();
for(i=2;i*i<x;i++) if(x%i==0) b[i].push_back(p),b[x/i].push_back(p);
if(i*i==x) b[i].push_back(p);b[x].push_back(p);
}
for(int i=l;i<N;i++){
if(b[i].size()>1)
for(int j=0;j<b[i].size()-1;j++)
add(b[i][j],b[i][j+1]);
}
topsort();
printf("%d",ans);
}