[线段树 & 前缀 优化建图 二分 2-SAT] CF Gym100159 facebook-hacker-cup-2012 I. Unfriending

二分答案

对于一个点，肯定有一个区间里的点和它不能共存

这个可以用线段树优化建图，每个friend list最多只能删一个，用前缀优化建图

还要一些针对随机数据什么的trick才能过…比如二分中的那个特判

#include <cstdio>#include <iostream>#include <algorithm>#include <vector>#include <cstring>using namespace std;const int N=5000010;int T,n,m,cnt,tot;int G[N],vis[N],dfn[N],low[N],s[N],g[N],vst[N],ig,tp,it,imax;struct edge{    int t,nx;}E[N<<1];pair<int,int> a[N];inline char nc(){    static char buf[100000],*p1=buf,*p2=buf;    return p1==p2&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++;}inline void rea(int &x){    char c=nc(); x=0;    for(;c>'9'||c<'0';c=nc());for(;c>='0'&&c<='9';x=x*10+c-'0',c=nc());}vector<pair<int,int> > b[N];inline void addedge(int x,int y){    //cout<<x<<' '<<y<<endl;    E[++cnt].t=y; E[cnt].nx=G[x]; G[x]=cnt;}void Build(int g,int l,int r){    if(l==r) return imax=max(imax,g),addedge(tot+g,a[l].second<<1|1);    int mid=l+r>>1;    Build(g<<1,l,mid); Build(g<<1|1,mid+1,r);    addedge(tot+g,tot+(g<<1)); addedge(tot+g,tot+(g<<1|1));}void link(int g,int L,int R,int l,int r,int x){    if(L==l && r==R) return addedge(x,tot+g);    int mid=L+R>>1;    if(r<=mid) link(g<<1,L,mid,l,r,x);    else if(l>mid) link(g<<1|1,mid+1,R,l,r,x);    else link(g<<1,L,mid,l,mid,x),link(g<<1|1,mid+1,R,mid+1,r,x);}void tarjan(int x){    dfn[x]=low[x]=++it; vis[x]=1; s[++tp]=x;    for(int i=G[x];i;i=E[i].nx){        if(!vis[E[i].t]) tarjan(E[i].t);        if(vis[E[i].t]!=2) low[x]=min(low[x],low[E[i].t]);    }    if(dfn[x]==low[x]){        int k; ++ig;        do{ k=s[tp--]; vis[k]=2; g[k]=ig; }while(k!=x && tp);    }}inline bool check(int x){    memset(G,0,sizeof(G));    memset(vis,0,sizeof(vis));    cnt=tp=ig=it=imax=0; tot=(n-1<<1|1)+1;    Build(1,0,n-1);    int cc=0;    for(int i=0;i<n;i++)        if(!vst[i]) addedge(i<<1|1,i<<1);    for(int i=0,j=0;i<n;i++){        while(a[i].first-a[j].first>=x) j++;        if(j<i) link(1,0,n-1,j,i-1,a[i].second<<1),cc+=i-j;    }    if(cc>=1200000) return false;    for(int i=n-1,j=n-1;~i;i--){        while(a[j].first-a[i].first>=x) j--;        if(j>i) link(1,0,n-1,i+1,j,a[i].second<<1);    }    tot+=imax+1;    for(int i=1;i<=m;i++){        addedge(b[i][0].second<<1|1,tot);        addedge(tot+1,b[i][0].second<<1);        for(int j=1;j<b[i].size();j++){            addedge(tot+(j-1<<1),tot+(j<<1));            addedge(tot+(j<<1|1),tot+(j-1<<1|1));            addedge(b[i][j].second<<1|1,tot+(j<<1));            addedge(tot+(j<<1|1),b[i][j].second<<1);            addedge(b[i][j].second<<1|1,tot+(j-1<<1|1));            addedge(tot+(j-1<<1),b[i][j].second<<1);        }        tot+=(b[i].size()<<1|1)+1;    }    for(int i=0;i<=tot;i++)        if(!vis[i]) tarjan(i);    for(int i=0;i<n;i++)        if(g[i<<1]==g[i<<1|1]) return false;    return true;}int main(){    rea(T); int CASE=0;    while(T--){        memset(vst,0,sizeof(vst));        rea(n); rea(m); rea(a[0].first);        int x,y,p; rea(x); rea(y); rea(p);        for(int i=1;i<n;i++)             a[i].first=(1LL*a[i-1].first*x+y)%p;        for(int i=0;i<n;i++)             a[i].second=i;        for(int i=1;i<=m;i++){            int size,x,y; rea(size);            b[i].resize(size);            rea(b[i][0].second); rea(x); rea(y);            b[i][0].second%=n;            for(int j=1;j<size;j++)                 b[i][j].second=(1LL*b[i][j-1].second*x+y)%n;            for(int j=0;j<size;j++)                 b[i][j].first=a[b[i][j].second].first,vst[b[i][j].second]=1;            sort(b[i].begin(),b[i].end());            b[i].resize(unique(b[i].begin(),b[i].end())-b[i].begin());        }        sort(a,a+n);        int l=0,r=1e9,mid,ans=-1;        while(l<=r){            mid=l+r>>1;            if(l<10 && r>10) mid=10;            check(mid)?l=(ans=mid)+1:r=mid-1;        }        printf("Case #%d: %d\n",++CASE,ans);    }    return 0;}