Codeforces Round #147 (Div. 2) E Build String(最小费用流)

题意:给你一个目标串str_t,你可以用下面的n个串(str[N][N]),每个串str[i]中选出limit[i]个字符,从第i个串中拿出一个字符消耗代价是i。问构成这个串最小的代价是多少,如果不能构成这个目标串str_t则输出-1。

设置一个超级源点st和超级汇点ed。超级源点到每个str[i]建立一条容量为limit[i],代价为0的边。将目标串拆成26个字母,根据目标串中有多少个这个字符,建立从这个字符到超级汇点的边,容量即有它的个数,代价也为0。为每个串str[i]中中每个字符建立从这个串str[i]到目标号对应字母的边,容量为1,代价即为i。

#include <iostream>
#include <cstdio>
#include <cstring>
#include <vector>
#include <string>
#include <set>
#include <map>
#include <cmath>
#include <algorithm>
#include <queue>
using namespace std;

const int MaxM=200001;
const int INF=(1<<30);

const int N=105;

struct Edge
{
	int u,v,pre,cap,cost;
	Edge(){}
	Edge(int u,int v,int pre,int cap,int cost) :
		u(u),v(v),pre(pre),cap(cap),cost(cost) {}
}edge[MaxM];

int head[MaxM],nEdge,re_flow;

void init()
{
	nEdge=0;
	memset(head,-1,sizeof(head));
}
void addEdge(int u,int v,int cap,int cost)
{
	edge[nEdge]=Edge(u,v,head[u],cap,cost);
	head[u]=nEdge++;
	edge[nEdge]=Edge(v,u,head[v],0,-cost);
	head[v]=nEdge++;
}
struct MinCostFlow
{
	queue<int> que;
	int vis[MaxM],pre[MaxM],dis[MaxM],pos[MaxM];
	int spfa(int s,int t,int n)
	{
		for(int i=0;i<=n;i++)
		{
			pre[i]=-1; vis[i]=0; dis[i]=INF;
		}

		que.push(s); pre[s]=s; dis[s]=0; vis[s]=1;
		while(!que.empty())
		{
			int u=que.front();

			que.pop(); vis[u]=0;

			for(int i=head[u];i!=-1;i=edge[i].pre)
			{
				int v=edge[i].v,cost=edge[i].cost;
				if(edge[i].cap>0&&dis[u]+cost<dis[v])
				{
					dis[v]=dis[u]+cost;

					pre[v]=u; pos[v]=i;

					if(!vis[v])
					{
						vis[v]=1;
						que.push(v);
					}
				}
			}
		}
		if(pre[t]!=-1&&dis[t]<INF) return 1;
		return 0;
	}
	void solve(int s,int t,int n,int &flow,int &cost)
	{
		flow=0,cost=0;
		while(spfa(s,t,n))
		{
			int mi=INF;
			for(int u=t;u!=s;u=pre[u]) mi=min(mi,edge[pos[u]].cap);
			flow+=mi;
			cost+=mi*dis[t];
			for(int u=t;u!=s;u=pre[u])
			{
				edge[ pos[u] ].cap-=mi;
				edge[ pos[u]^1 ].cap+=mi;
			}
		}
	}
}flow;


void input(int &st,int &ed)
{
    int n,len,limit[N],cnt[30]={0};
    char str_t[N],str[N][N];

    scanf("%s %d",str_t,&n);

    st=0;
    for(int i=1;i<=n;i++)
    {
        scanf("%s %d",str[i],&limit[i]);

        addEdge(st,i,limit[i],0);

        len=(int)strlen(str[i]);
        for(int j=0;j<len;j++)
        {
            addEdge(i,n+str[i][j]-'a'+1,1,i);
        }
    }

    re_flow=len=(int)strlen(str_t);

    for(int i=0;i<len;i++)
    {
        cnt[str_t[i]-'a']++;
    }
    ed=n+31;

    for(int i=0;i<30;i++)
    {
        if(cnt[i]) addEdge(n+i+1,ed,cnt[i],0);
    }
}
int main()
{
    init();
    int st,ed,mx_flow,cost;
    input(st,ed);
    flow.solve(st,ed,ed-st+1,mx_flow,cost);
    if(mx_flow!=re_flow) cout<<-1<<endl;
    else cout<<cost<<endl;
    return 0;
}



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

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值