POJ3275 [USACO07MAR]Ranking the Cows

Description

Each of Farmer John's N cows (1 ≤ N ≤ 1,000) produces milk at a different positive rate, and FJ would like to order his cows according to these rates from the fastest milk producer to the slowest.

FJ has already compared the milk output rate for M (1 ≤ M ≤ 10,000) pairs of cows. He wants to make a list of C additional pairs of cows such that, if he now compares those C pairs, he will definitely be able to deduce the correct ordering of all N cows. Please help him determine the minimum value of C for which such a list is possible.

Input

Line 1: Two space-separated integers: N and M
Lines 2..M+1: Two space-separated integers, respectively: X and Y. Both X and Y are in the range 1...N and describe a comparison where cow X was ranked higher than cow Y.

Output

Line 1: A single integer that is the minimum value of C.

Sample Input

5 5
2 1
1 5
2 3
1 4
3 4

Sample Output

3

Hint

From the information in the 5 test results, Farmer John knows that since cow 2 > cow 1 > cow 5 and cow 2 > cow 3 > cow 4, cow 2 has the highest rank. However, he needs to know whether cow 1 > cow 3 to determine the cow with the second highest rank. Also, he will need one more question to determine the ordering between cow 4 and cow 5. After that, he will need to know if cow 5 > cow 3 if cow 1 has higher rank than cow 3. He will have to ask three questions in order to be sure he has the rankings: "Is cow 1 > cow 3? Is cow 4 > cow 5? Is cow 5 > cow 3?"

 

第一想法用拓扑排序玩传递闭包，毕竟是稀疏图

#include<cstdio>
#include<cstdlib>
#include<stack>
#include<bitset>
#include<vector>
using namespace std;
const int N=1024;
stack<int> S;
bitset<N> a[N];
vector<int> G[N];
int n,m;
int d[N];
int read(){
    char c;int x=0,f=1;
    for(c=getchar();!isdigit(c);c=getchar())if(c=='-')f=-1;
    for(;isdigit(c);c=getchar())x=(x<<3)+(x<<1)+(c^48);
    return x*f;
}
void TopSort(){
    for(int i=0;i<n;++i)
        if(!d[i])S.push(i);
    while(!S.empty()){
        int u=S.top();S.pop();
        for(int i=0,i_end=G[u].size();i<i_end;++i){
            int v=G[u][i];
            a[v]|=a[u];
            --d[v];
            if(!d[v])S.push(v);
        }
    }
}
int main(){
    scanf("%d%d",&n,&m);
    while(m--){
        int u=read()-1,v=read()-1;
        G[u].push_back(v);
        ++d[v];
    }
    for(int i=0;i<n;++i)a[i][i]=1;
    TopSort();
    for(int i=0;i<n;++i)
        m+=a[i].count();
    printf("%d\n",n*(n-1)/2-m-1+n);
    return 0;
}

就当是对ZJOI2017Day2的交代吧