BZOJ 1283 序列 费用流 网络流 线性规划

https://darkbzoj.cf/problem/1283

给出一个长度为N的正整数序列Ci，求一个子序列，使得原序列中任意长度为M的子串中被选出的元素不超过K(K,M<=100) 个，并且选出的元素之和最大。

http://www.cnblogs.com/137shoebills/p/8871648.html

↑和这道题一样

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<queue>
using namespace std;
const int maxn=1100;
int n,m,k,S,T,SS;
struct nod{
    int x,y,v,co,rev,next;
}e[10*maxn];
int head[maxn]={},tot=0;
int a[maxn]={};
void init(int x,int y,int v,int co){
    e[++tot].x=x;e[tot].y=y;e[tot].v=v;e[tot].co=co;e[tot].rev=tot+1;e[tot].next=head[x];head[x]=tot;
    e[++tot].x=y;e[tot].y=x;e[tot].v=0;e[tot].co=-co;e[tot].rev=tot-1;e[tot].next=head[y];head[y]=tot;
}
queue<int>q;
int vis[maxn]={},fa[maxn]={},dis[maxn]={};
bool SPFA(){
    memset(dis,31,sizeof(dis));
    memset(vis,0,sizeof(vis));
    memset(fa,0,sizeof(fa));
    q.push(S);vis[S]=1;dis[S]=0;
    while(!q.empty()){
        int x=q.front();q.pop();vis[x]=0;
        for(int i=head[x];i;i=e[i].next){
            if(!e[i].v)continue;
            if(dis[e[i].y]>dis[x]+e[i].co){
                dis[e[i].y]=dis[x]+e[i].co;
                fa[e[i].y]=i;
                if(!vis[e[i].y]){
                    q.push(e[i].y);vis[e[i].y]=1;
                }
            }
        }
    }
    return fa[T];
}
int doit(){
    int val=maxn,ans=0;
    for(int i=fa[T];i;i=fa[e[i].x])val=min(val,e[i].v);
    for(int i=fa[T];i;i=fa[e[i].x]){
        e[i].v-=val;e[e[i].rev].v+=val;ans+=e[i].co;
    }
    return ans;
}
int main(){
    scanf("%d%d%d",&n,&m,&k);
    for(int i=1;i<=n;i++)scanf("%d",&a[i]);
    T=n+1;S=T+1;SS=S+1;
    for(int i=1;i<=n;i++){
        init(i,i+m>n?T:i+m,1,-a[i]);
        init(i,i+1>n?T:i+1,k,0);
    }
    for(int i=1;i<=m;i++)init(SS,i,maxn,0);
    init(S,SS,k,0);
    int ans=0;
    while(SPFA())ans-=doit();
    printf("%d\n",ans);
    return 0;
}