【费用流】[CodeForces - 717G]Underfail
题目大意
题目大概说给一个主串和几个有价值的模式串,某个模式串与主串匹配就能累加对应的价值,一个模式串可以在多个位置和主串匹配但同一个位置只能一次,此外主串各个字符最多可以用x次,问如何匹配使获得的价值最大。
分析
暴力匹配模式串在主串中的位置,然后在匹配区间的左端点和右端点+1的地方连一条边,容量为1,费用为匹配这个模式串的收益。
然后在相邻的两个位置之间连边,容量为1,费用为0。
跑一次最大费用流即可。
代码
#include<cstdio>
#include<algorithm>
#include<queue>
#include<cstring>
using namespace std;
#define MAXN 500
#define MAXM 100
#define INF 0x7f7f7f7f
queue<int>q;
char s[MAXN+10],t[MAXN+10];
int n,T,S,tot,dist[MAXN+10],m,x;
bool vis[MAXN+10];
void Read(int &x){
char c;
while(c=getchar(),c!=EOF)
if(c>='0'&&c<='9'){
x=c-'0';
while(c=getchar(),c>='0'&&c<='9')
x=x*10+c-'0';
ungetc(c,stdin);
return;
}
}
struct node{
int v,cap,wt;
node *next,*back;
}*adj[MAXN+10],edge[MAXN*MAXM*2+10],*ecnt=edge,*pre[MAXN+10];
inline void addedge(int u,int v,int cap,int wt){
node *p=++ecnt;
p->v=v;
p->wt=wt;
p->cap=cap;
p->next=adj[u];
adj[u]=p;
p=p->back=++ecnt;
p->v=u;
p->wt=-wt;
p->cap=0;
p->next=adj[v];
adj[v]=p;
p->back=ecnt-1;
}
void read(){
Read(n);
scanf("%s",s+1);
int i,j,len,p;
bool flag=0;
Read(m);
while(m--){
scanf("%s",t+1);
Read(p);
len=strlen(t+1);
for(i=1;i+len-1<=n;i++){
flag=0;
for(j=1;j<=len;j++)
if(s[i+j-1]!=t[j]){
flag=1;
break;
}
if(!flag)
addedge(i,i+len,1,-p);
}
}
Read(x);
for(i=1;i<n;i++)
addedge(i,i+1,x,0);
T=n+1,S=tot=T+1;
addedge(S,1,x,0);
addedge(n,T,x,0);
}
bool spfa(){
int u;
int i;
for(i=1;i<=tot;i++)
dist[i]=INF;
dist[S]=0;
q.push(S);
while(!q.empty()){
u=q.front();
q.pop();
vis[u]=0;
for(node *p=adj[u];p;p=p->next){
if(p->cap&&dist[p->v]>dist[u]+p->wt){
pre[p->v]=p;
dist[p->v]=dist[u]+p->wt;
if(!vis[p->v])
vis[p->v]=1,q.push(p->v);
}
}
}
return dist[T]!=INF;
}
int mcmf(){
int cost=0,delta;
node *p;
while(spfa()){
delta=INF;
for(p=pre[T];p;p=pre[p->back->v])
delta=min(delta,p->cap);
for(p=pre[T];p;p=pre[p->back->v]){
p->cap-=delta;
p->back->cap+=delta;
}
cost+=dist[T]*delta;
}
return cost;
}
int main()
{
read();
printf("%d\n",-mcmf());
}