poj 2699(最大流dinic终极模版)

题意:

有N个参赛者,两两比赛,赢一场比赛得一分,得分最高的为king,此外,如果一名参赛者比赛中打败了所有得分比他高的参赛者,那么他也是一个king.....题目给出一个非递减的得分序列,表示这N个参赛者的得分情况,按照比赛的得分结果,求最多可能有多少个king.

思路:

以每名参赛者和每场比赛都为一点,所有参赛者的点连接源点,容量为各个分数,所有比赛连接一汇点,容量为1;然后从king的数目num枚举,从最大(N)开始,即每个人都为king,判断这种情况可不可能符合分数序列.

建图如下:得分不是最高(如i)的要成为King就得打败所有比他分高者(如j),所以他必须赢得那场比赛,连接那场比赛代表的点(v),即i->v,容量为1,而j->v为0. 剩下没有确定胜负的比赛则双方都可以赢,即i->v,j->v,容量为1;

然后求最大流,如果等于比赛数目,即表示这种情况有可能发生,则得到最多King值.否则,就减少num值,再判断,直至结果出来或者只有一个得分最高的为king.

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
using namespace std;
const int maxn=550;
const int inf=1<<20;
//******************************************//
struct Node
{
    int v,next;
    int val;
}edge[maxn*maxn*2];

int level[maxn];//顶点的层次
int head[maxn];
int que[maxn*10];//BFS中用于遍历的顶点，DFS求增广中记录边
int out[10*maxn];//DFS用于几乎定点的分支
int ind;

void init()
{
    ind=0;
    memset(head,-1,sizeof(head));
}
void addedge(int x,int y,int z)
{
    edge[ind].v=y;
    edge[ind].val=z;
    edge[ind].next=head[x];
    head[x]=ind++;
    edge[ind].v=x;
    edge[ind].val=0;
    edge[ind].next=head[y];
    head[y]=ind++;
}
int dinic(int n,int source,int sink)
{
    int ret=0;
    int h=0,r=0;
    while(1)//DFS
    {
        int i;
        for(i=0;i<=n;++i)
            level[i]=0;
        h=0,r=0;
        level[source]=1;
        que[0]=source;
        while(h<=r)//BFS
        {
            int  t=que[h++];
            for(i=head[t];i!=-1;i=edge[i].next)
                {
                    if(edge[i].val&&level[edge[i].v]==0)
                        {
                          level[edge[i].v]=level[t]+1;
                          que[++r]=edge[i].v;
                        }
                }
        }

        if(level[sink]==0)break;//找不到汇点
        for(i=0;i<=n;++i)
          out[i]=head[i];

        int  q=-1;
        while(1)
        {
            if(q<0)
            {
                int cur=out[source];
                for(;cur!=-1;cur=edge[cur].next)
                {
                    if(edge[cur].val&&out[edge[cur].v]!=-1&&level[edge[cur].v]==2)
                    {
                        break;
                    }
                }
                if(cur>=0)
                {
                    que[++q]=cur;
                    out[source]=edge[cur].next;
                }
                else
                    break;
            }
            int  u=edge[que[q]].v;
            if(u==sink)//一条增广路
            {
                int  dd=inf;
                int  index=-1;
                for(i=0;i<=q;i++)
                {
                    if(dd>edge[que[i]].val)
                    {
                        dd=edge[que[i]].val;
                        index=i;
                    }
                }
                ret+=dd;
                for(i=0;i<=q;i++)
                {
                    edge[que[i]].val-=dd;
                    edge[que[i]^1].val+=dd;
                }
                for(i=0;i<=q;i++)
                {
                    if(edge[que[i]].val==0)
                    {
                        q=index-1;
                        break;
                    }
                }
            }
            else
            {
                int cur=out[u];
                for(;cur!=-1;cur=edge[cur].next)
                {
                    if(edge[cur].val&&out[edge[cur].v]!=-1&&level[u]+1==level[edge[cur].v])
                    {
                        break;
                    }
                }
                if(cur!=-1)
                {
                    que[++q]=cur;
                    out[u]=edge[cur].next;
                }
                else
                {
                    out[u]=-1;
                    q--;
                }
            }
        }
    }
    return ret;
}
//**************************************************//


int n,m;
int score[maxn];
char str[200];
string x[maxn];
void input()  //恶心的输入
{
    int i;
    char c;
    for(i = 1; i <= 10; i++)
        x[i] = "";
    n = 0;
    while(1)
    {
        bool flag = true;
        while((c = getchar()) != ' ')
        {
            if (c == '\n')
            {
                n++;
                for(i = 1; i <= n; i++)
                    sscanf(x[i].c_str(),"%d",&score[i]);
                return;
            }
            if (flag)
            {
                flag = false;
                n++;
            }
            x[n] += c;
        }
    }
}

int main()
{
    int i, j, k, ca, ans;
    scanf("%d",&ca);
    getchar();
    while(ca--)
    {
        input();
        n--;
        m = n*(n-1)/2;
        int source = 0;
        int sink = n+m+1;
        ans = 0;
        for(k = 0; k <= n; k++)
        {
            init();
            for(i = 1; i <= n; i++)
                addedge(source,i,score[i]);
            int num = 0;
            for(i = 1; i <= n; i++)
            {
                for(j = i+1; j <= n; j++)
                {
                    num++;
                    addedge(num+n,sink,1);
                    if(i > n-k && score[j] > score[i])
                        addedge(i,num+n,1);
                    else
                    {
                        addedge(i,num+n,1);
                        addedge(j,num+n,1);
                    }
                }
            }
            if(dinic(n+m+2,source,sink) == m)
                ans = k;
            else
                break;
        }
        printf("%d\n",ans);
    }
    return 0;
}