UVA 12433 Rent a Car

在网上查貌似没有这道题的题解，大致写了一下

# include <string.h>
# include <algorithm>
# include <queue>
using namespace std;
const int INF=0x3f3f3f3f;
const int MAXN=205;

/*转化为最小费用最大流问题：
假设有5天的租车需求，虚拟出2*n+2 即 12个节点，0为源点，12为汇点。
源点到1 2 3 4 5流量为r[i],费用为0。6 7 8 9 10到汇点流量为r[i-n]，费用为0。此题为一个检验能否满流且求满流花费最小的问题
虚拟第2*n+1个节点为买车途径，源点到2*n+1节点花费为p[i],流量为c[i]，多重边。对于每一个i+n节点，其来源有两个，一个是
2*n+1节点，即购买新车，一个是之前的车辆送去维修后的可用车辆。连接i和i+1节点，流量为INF，花费为0。同时根据维修的天数
连接i和d[j]+i+1+n节点，花费为s[i]，流量为INF。*/

//采用MCMF最小费用最大流算法，每次找花费最小的增广路。

struct Edge
{
    int to,next,cap,flow,cost;
} edge[100010];
int head[MAXN],cnt=0;
int pre[MAXN],dis[MAXN];
int r[55],c[55],p[55],d[55],s[55],n,cost;
bool vis[MAXN];

void addedge(int u,int v,int cap,int cost)
{
    edge[cnt].to=v;
    edge[cnt].cap=cap;
    edge[cnt].cost=cost;
    edge[cnt].flow=0; //记录当前流量
    edge[cnt].next=head[u];
    head[u]=cnt++;
    edge[cnt].to=u;
    edge[cnt].cap=0;
    edge[cnt].cost=-cost;
    edge[cnt].flow=0;
    edge[cnt].next=head[v];
    head[v]=cnt++;
}

bool spfa(int s,int t)
{
    queue<int>q;
    for(int i=0; i<=2*n+2; i++)
    {
        dis[i]=INF;
        vis[i]=0;
        pre[i]=-1;
    }
    dis[s]=0;
    vis[s]=1;
    q.push(s);
    while(!q.empty())
    {
        int u=q.front();
        q.pop();
        vis[u]=0;
        for(int i=head[u]; i!=-1; i=edge[i].next)
        {
            int v=edge[i].to;
            if(edge[i].cap>edge[i].flow && dis[v]>dis[u]+edge[i].cost)
            {
                dis[v]=dis[u]+edge[i].cost;
                pre[v]=i;
                if(!vis[v])
                {
                    vis[v]=1;
                    q.push(v);
                }
            }
        }
    }
    //return dis[t]!=INF;
    if(pre[t]==-1)return 0;
    else return 1;
}
int MCMF(int s,int t)
{
    int flow=0;
    cost=0;
    while(spfa(s,t))
    {
        int Min=INF;
        for(int i=pre[t]; i!=-1; i=pre[edge[i^1].to]) //求增广的流量
        {
            if(Min>edge[i].cap-edge[i].flow)
                Min=edge[i].cap-edge[i].flow;
        }
        for(int i=pre[t]; i!=-1; i=pre[edge[i^1].to])
        {
            edge[i].flow+=Min;
            edge[i^1].flow-=Min;
            //printf("%d  %d\n",i,i^1);
            cost+=edge[i].cost*Min;
        }
        flow+=Min;
    }
    return flow;
}


int main ()
{
    int u,j,t,C,R,ans;
    scanf("%d",&t);
    for(int ca=1; ca<=t; ca++)
    {
        ans=0;cnt=0;
        scanf("%d%d%d",&n,&C,&R);
        for(int i=1; i<=n; i++)scanf("%d",&r[i]),ans+=r[i];
        for(int i=1; i<=C; i++)scanf("%d%d",&c[i],&p[i]);
        for(int i=1; i<=R; i++)scanf("%d%d",&d[i],&s[i]);
        memset(head,-1,sizeof(head));
        for(int i=1; i<=n; i++)
            addedge(0,i,r[i],0);
        for(int i=n+1; i<=2*n; i++)
            addedge(i,(n+1)*2,r[i-n],0);
        for(int i=1; i<=C; i++)
            addedge(0,2*n+1,c[i],p[i]);
        for(int i=n+1; i<=2*n; i++)
            addedge(2*n+1,i,INF,0);
        for(int i=1; i<n; i++)
            addedge(i,i+1,INF,0);
        for(int i=1; i<=n; i++)
        {
            for(int j=1; j<=R; j++)
            {
                if(i+d[j]+1<=n) addedge(i,i+d[j]+1+n,INF,s[j]);
            }
        }
        int f=MCMF(0,2*n+2);
        printf("Case %d: ",ca);
        if(f!=ans) printf("impossible\n");
        else printf("%d\n",cost);
    }
    return 0;
}