HDU - 6437：Videos （裸的费用流）

。。。懒得说什么了

#include<bits/stdc++.h>
using namespace std;
const int maxn=;
const int inf=<<;int To[maxn],Laxt[maxn],Next[maxn],cap[maxn],cost[maxn];
int S,T,cnt,dis[maxn],ans;
bool inq[maxn],vis[maxn];
deque<int>q;
void add(int u,int v,int c,int cc)
{
    Next[++cnt]=Laxt[u];Laxt[u]=cnt; To[cnt]=v;cap[cnt]=c;cost[cnt]=-cc;
    Next[++cnt]=Laxt[v];Laxt[v]=cnt; To[cnt]=u;cap[cnt]=;cost[cnt]=cc;
}
bool spfa()
{
    for(int i=;i<=T;i++) inq[i]=;
    for(int i=;i<=T;i++) dis[i]=inf;
    inq[T]=; dis[T]=; q.push_back(T);
    while(!q.empty())
    {
        int u=q.front(); q.pop_front();
        inq[u]=;
        for(int i=Laxt[u];i;i=Next[i])
        {
            int v=To[i];
            if(cap[i^]&&dis[v]>dis[u]-cost[i])
            {
                dis[v]=dis[u]-cost[i];
                if(!inq[u]){
                    inq[v]=;
                    if(q.empty()||dis[v]>dis[q.front()]) q.push_back(v);
                    else q.push_front(v);
                }
            }
        }
    }
    return dis[S]<inf;
}
int dfs(int u,int flow)
{
    vis[u]=;
    if(u==T||flow==) return flow;
    int tmp,delta=;
    for(int i=Laxt[u];i;i=Next[i])
    {
        int v=To[i];
        if((!vis[v])&&cap[i]&&dis[v]==dis[u]-cost[i])
        {
            tmp=dfs(v,min(cap[i],flow-delta));
            delta+=tmp; cap[i]-=tmp; cap[i^]+=tmp;
        }
    }
    return delta;
}
int s[maxn],t[maxn],w[maxn],c[maxn];
int main()
{
    int N,M,K,W,C,i,j;
    scanf("%d",&C);
    while(C--){
        scanf("%d%d%d%d",&N,&M,&K,&W);
        cnt=;  ans=; T=M+M+;
        for(int i=;i<=T;i++)  vis[i]=Laxt[i]=;
        for(int i=;i<=M;i++) scanf("%d%d%d%d",&s[i],&t[i],&w[i],&c[i]);
        for(int i=;i<=M;i++)
         for(int j=;j<=M;j++)
          if(i!=j&&t[i]<=s[j]) add(i+M,j,,c[i]==c[j]?-W:);
        add(S,T-,K,); //限制N人
        for(int i=;i<=M;i++) add(T-,i,,); //收益
        for(int i=;i<=M;i++) add(i,i+M,,w[i]);
        for(int i=;i<=M;i++) add(i+M,T,,);
        while(spfa()){
           vis[T]=;
           while(vis[T]){
              for(i=;i<=T;i++) vis[i]=;
              ans-=dis[S]*dfs(S,N);
           }
        }
        printf("%d\n",ans);
    }
    return ;
}