zju monthly contest 2012 aug.



#include <cstdio>
#include <cstring>
#include <map>
#include <algorithm>

using namespace std;
#pragma comment(linker, "/STACK:1024000000,1024000000")
const int maxn=40000;
const int inf=0x7fffffff;
const int s=0;

struct edge{
       int v,next,w;
}edge[1000000];
int head[2*maxn],cnt;//for sap

void addedge(int u, int v, int w)
{
     edge[cnt].v=v;
     edge[cnt].w=w;
     edge[cnt].next=head[u];
     head[u]=cnt++;
     edge[cnt].v=u;
     edge[cnt].w=0;
     edge[cnt].next=head[v];
     head[v]=cnt++;
}

int sap(int t)
{
    int pre[2*maxn],cur[2*maxn],dis[2*maxn],gap[2*maxn];
    int flow=0 , aug=inf ,u;
    bool flag;
    for (int i=0 ; i<=t ; ++i)
    {
        cur[i]=head[i];
        gap[i]=dis[i]=0;
    }
    gap[s]=t+1;
    u=pre[s]=s;
    while (dis[s]<=t)
    {
          flag=0 ;
          for (int &j=cur[u] ; ~j ; j=edge[j].next)
          {
              int v=edge[j].v;
              if (edge[j].w>0 && dis[u]==dis[v]+1)
              {
                   flag=1;
                   if(edge[j].w<aug)aug=edge[j].w;
                   pre[v]=u;
                   u=v;
                   if (u==t)
                   {
                       flow+=aug;
                       while (u!=s)
                       {
                             u=pre[u];
                             edge[cur[u]].w-=aug;
                             edge[cur[u]^1].w+=aug;
                       }
                       aug=inf;
                   }
                   break;
              }
          }
          if (flag)continue ;
          int mindis=t+1;
          for (int j=head[u]; ~j ; j=edge[j].next)
          {
              int v=edge[j].v;
              if (edge[j].w>0 && dis[v]<mindis)
              {
                 mindis=dis[v];
                 cur[u]=j;
              }
          }
          if(--gap[dis[u]]==0)break;
          gap[dis[u]=mindis+1]++;
          u=pre[u];
    }
    return flow;
}

void init ()
{
     memset (head , -1 , sizeof(head));
     cnt=0;
}

map <int , int >idx;

int ct[234];
int info[234][234];
int c[234];

int main ()
{
    int n;
    while (~scanf("%d", &n))
    {
        init();
        idx.clear();
        int Id=n;
        for (int i=1; i<=n; ++i)
        {
            int b; scanf("%d %d %d", ct+i, &b, c+i);
            addedge(0, i, c[i]);
            for (int j=0; j<ct[i]; ++j)
            {
                scanf("%d", &info[i][j]);
                if(idx[info[i][j]]==0)idx[info[i][j]]=++Id;
            }
        }
        for (int i=1; i<=n; ++i)
        {
            for (int j=0; j<ct[i]; ++j)
            addedge(i, idx[info[i][j]], 1);
        }
        for (int i=n+1; i<=Id; ++i) addedge(i, Id+1, 1);
        int m; scanf("%d", &m);
        addedge(0, m, ct[m]-c[m]);
        int ans=sap(Id+1);
        printf("%d\n", ans);
    }
    return 0;
}



  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值