P4376 [USACO18OPEN]Milking Order

链接：https://www.luogu.org/problemnew/show/P4376

题意：中文题面，自读

思路：因为是从上到下连续取x个观测状态，那么就可以二分x，得到一个最大的x。如何去判断取前x个观测条件满不满足条件呢，只需要用下拓扑判断一下有没有环就可以了，有环显然是矛盾的（不合法的）。

复杂度差不多是(n+edge)logmlogn，差不多是几e7。一开始在别的oj交，很轻松就过了，然后洛谷里就一直超时。后来发现二分check部分的拓扑排序根本不需要优先队列，判断有没有环就可以，不需要顺序。只有最后找到x最大值后，这个时候用优先队列的拓扑跑个答案出来就好了，复杂度少了个log。

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef pair<int,int > pii;
typedef pair<char,char> pcc;
vector<int>ve[50005];

int head[100005],cnt;
int n,m;
struct node{
    int v,nxt;
}edge[200005];

void addedge(int u,int v){
    edge[++cnt].v=v;
    edge[cnt].nxt=head[u];
    head[u]=cnt;
}
bool vis[100005];
int du[100005];

vector<int>ans;
bool topu(){
    int sz=0;
    queue<int>q;
    for(int i=1;i<=n;i++){
        if(!du[i]){
            sz++;
            q.push(i);
        }
    }
    while(!q.empty()){
        int u=q.front();
        q.pop();
        for(int i=head[u];~i;i=edge[i].nxt){
            int v=edge[i].v;
            du[v]--;
            if(!du[v])q.push(v),vis[v]=true,sz++;
        }
    }
    return sz==n;
}

bool judge(int m){
    memset(head,-1,sizeof(head));
    memset(du,0,sizeof(du));
    cnt=0;
    for(int i=1;i<=m;i++){
        if(ve[i].size()==0)continue;
        int u=ve[i][0];
        for(int j=1;j<ve[i].size();j++){
            int v=ve[i][j];
            addedge(u,v);
            du[v]++;
            u=v;
        }
    }
    return topu();
}

void getans(int x){
    memset(head,-1,sizeof(head));
    memset(du,0,sizeof(du));
    cnt=0;
    for(int i=1;i<=x;i++){
        if(ve[i].size()==0)continue;
        int u=ve[i][0];
        for(int j=1;j<ve[i].size();j++){
            int v=ve[i][j];
            addedge(u,v);
            du[v]++;
            u=v;
        }
    }
    priority_queue<int,vector<int>,greater<int> >q;
    for(int i=1;i<=n;i++){
        if(!du[i]){
            q.push(i);
        }
    }
    while(!q.empty()){
        int u=q.top();
        q.pop();
        ans.push_back(u);
        for(int i=head[u];~i;i=edge[i].nxt){
            int v=edge[i].v;
            du[v]--;
            if(!du[v])q.push(v);
        }
    }
}

template <typename _Tp> inline _Tp read(_Tp&x){
    char c11=getchar(),ob=0;x=0;
    while(c11^'-'&&!isdigit(c11))c11=getchar();if(c11=='-')c11=getchar(),ob=1;
    while(isdigit(c11))x=x*10+c11-'0',c11=getchar();if(ob)x=-x;return x;
}

int main(){

    read(n);read(m);
    for(int i=1,x,y;i<=m;i++){
        read(x);
        for(int j=1;j<=x;j++){
            read(y);
            ve[i].push_back(y);
        }
    }
    int x=0;
    int l=0,r=m;
    while(l<=r){
        int mid=(l+r)>>1;
        if(judge(mid))x=mid,l=mid+1;
        else r=mid-1;
    }
    getans(x);
    int ed=ans.size()-1;
    for(int i=0;i<=ed;i++){
        printf("%d%c",ans[i],i==ed?'\n':' ');
    }
    return 0;
}