BZOJ1283 序列

这个算是经典建模裸题了吧?然而还是没秒掉

使得每个长度为m的子串中选不超过k个,可以视为选k个子序列,每次选的时候相邻两个的位置差>m

那么每个点连后一个点费用0流量INF,连往后数m个点(如果超数组长度了就连T)费用a[i]流量1,S连1号点流量k费用0,跑最大费用流即可

#include<iostream>
#include<cstring>
#include<ctime>
#include<cmath>
#include<algorithm>
#include<iomanip>
#include<cstdlib>
#include<cstdio>
#include<map>
#include<bitset>
#include<set>
#include<stack>
#include<vector>
#include<queue>
using namespace std;
#define MAXN 1010
#define MAXM 10010
#define ll long long
#define eps 1e-8
#define MOD 1000000007
#define INF 1000000000
struct vec{
	int to;
	int fro;
	int f;
	int v;
};
vec mp[MAXM];
int tai[MAXN],cnt=1;
int answ;
int s,t;
int dis[MAXN];
int rd[MAXN];
int q[MAXN],hd,tl;
bool iq[MAXN];
inline void be(int x,int y,int f,int v){
	mp[++cnt].to=y;
	mp[cnt].fro=tai[x];
	tai[x]=cnt;
	mp[cnt].f=f;
	mp[cnt].v=v;
}
inline void bse(int x,int y,int f,int v){
	be(x,y,f,v);
	be(y,x,0,-v);
}
bool spfa(){
	int i,x,y;
	memset(dis,0xef,sizeof(dis));
	dis[s]=0;
	hd=tl=0;
	iq[q[tl++]=s]=1;
	while(hd!=tl){
		iq[x=q[(hd%=MAXN)++]]=0;
		for(i=tai[x];i;i=mp[i].fro){
			y=mp[i].to;
			if(dis[y]<dis[x]+mp[i].v&&mp[i].f){
				dis[y]=dis[x]+mp[i].v;
				rd[y]=i;
				if(!iq[y]){
					iq[q[(tl%=MAXN)++]=y]=1;
				}
			}
		}
	}
	if(dis[t]<0){
		return 0;
	}
	int flow=INF;
	for(i=t;i!=s;i=mp[rd[i]^1].to){
		flow=min(flow,mp[rd[i]].f);
	}
	answ+=flow*dis[t];
	for(i=t;i!=s;i=mp[rd[i]^1].to){
		mp[rd[i]].f-=flow;
		mp[rd[i]^1].f+=flow;
	}
	return 1;
}
int n,m,k;
int main(){
	int i,x;
	scanf("%d%d%d",&n,&m,&k);
	for(i=1;i<=n;i++){
		scanf("%d",&x);
		bse(i,min(n+1,i+m),1,x);
		bse(i,i+1,INF,0);
	}
	s=n+2;
	t=n+1;
	bse(s,1,k,0);
	while(spfa()){
		
	}
	printf("%d\n",answ);
	return 0;
}

/*

*/


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

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值