BZOJ1283 序列（费用流）

　　不妨看做是先用k个指针指向被选择的前k个元素，然后每次将选中当前第一个元素的指针移到最后，并且需要满足位置变化量>=m。显然这样可以构造出所有的合法方案。那么可以以此建立费用流模型，以一条流量k费用0的链将所有点串起来，再由每个位置向该位置+m连流量1费用为该元素权值的边，最大费用流即可。

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include<algorithm>
using namespace std;
int read()
{
    int x=,f=;char c=getchar();
    while (c<''||c>'') {if (c=='-') f=-;c=getchar();}
    while (c>=''&&c<='') x=(x<<)+(x<<)+(c^),c=getchar();
    return x*f;
}
#define N 1010
#define S 0
#define T 1001
int n,m,k,p[N],t=-,ans=;
int d[N],q[N],pre[N];
bool flag[N];
struct data{int to,nxt,cap,flow,cost;
}edge[N<<];
void addedge(int x,int y,int z,int c)
{
    t++;edge[t].to=y,edge[t].nxt=p[x],edge[t].cap=z,edge[t].flow=,edge[t].cost=c,p[x]=t;
    t++;edge[t].to=x,edge[t].nxt=p[y],edge[t].cap=,edge[t].flow=,edge[t].cost=-c,p[y]=t;
}
int inc(int &x){x++;if (x>n+) x-=n+;return x;}
bool spfa()
{
    memset(d,,sizeof(d));d[S]=;
    memset(flag,,sizeof(flag));
    int head=,tail=;q[]=S;
    do
    {
        int x=q[inc(head)];flag[x]=;
        for (int i=p[x];~i;i=edge[i].nxt)
        if (d[x]+edge[i].cost<d[edge[i].to]&&edge[i].flow<edge[i].cap)
        {
            d[edge[i].to]=d[x]+edge[i].cost;
            pre[edge[i].to]=i;
            if (!flag[edge[i].to]) q[inc(tail)]=edge[i].to,flag[edge[i].to]=;
        }
    }while (head!=tail);
    return d[T]<=;
}
void ekspfa()
{
    while (spfa())
    {
        int v=n;
        for (int i=T;i!=S;i=edge[pre[i]^].to)
        v=min(v,edge[pre[i]].cap-edge[pre[i]].flow);
        for (int i=T;i!=S;i=edge[pre[i]^].to)
        ans-=v*edge[pre[i]].cost,edge[pre[i]].flow+=v,edge[pre[i]^].flow-=v;
    }
}
int main()
{
#ifndef ONLINE_JUDGE
    freopen("bzoj1283.in","r",stdin);
    freopen("bzoj1283.out","w",stdout);
    const char LL[]="%I64d\n";
#else
    const char LL[]="%lld\n";
#endif
    n=read(),m=read(),k=read();
    memset(p,,sizeof(p));
    for (int i=;i<=n;i++)
    {
        int x=read();
        addedge(i-,i,k,);
        addedge(i,i+m>n?T:i+m,,-x);
    }
    addedge(n,T,k,);
    ekspfa();
    cout<<ans;
    return ;
}