BZOJ.4842.[NEERC2016]Delight for a Cat(费用流)

BZOJ


参考这儿

首先如果一个活动的时间满足条件,那么另一个活动也一定满足。还有就是这题就是费用流没有为什么。不妨假设最初所有时间都用来睡觉,那么我们要对每个\(k\)大小区间选出\([t2,k-t1]\)个时刻打游戏,同时要让选出的\(\sum s_i-e_i\)尽量小。

假如要求打游戏的时间恰好为\(t2\)怎么做呢?

对每个时刻建一个点\(i\),由\(i\)向\(i+k\)连边(不存在则直接连向汇点),容量\(1\)费用\(s_i-e_i\)。建一个虚点\(p\),\(p\)向\(1\sim k\)每个点连容量\(1\)费用\(0\)的边,源点\(S\)向\(p\)连容量\(t2\)费用\(0\)的边。

这样如果选了\(i\)就一定会选\(i+k\),所以每个\(k\)区间都会恰好选了\(t2\)个时刻。

那现在打游戏的时间可以提高到\(k-t1\)呢。

由\(i\)向\(i+1\)连边,容量\(k-t1-t2\)费用\(0\)。\(S\)向\(p\)的连边容量改成\(k-t1\)。

这样对于每个\(k\)区间,应该有的\(t2\)条打游戏的边和之前一样,不会少走。还可以多选\(k-t1-t2\)个时刻打游戏。

线性规划的做法实在是不想看惹...出了也做不出来


//980kb	4212ms
#include <queue>
#include <cstdio>
#include <cctype>
#include <cstring>
#include <algorithm>
#define gc() getchar()
typedef long long LL;
const int N=1007,M=3*N*2;
const LL INF=0x3f3f3f3f3f3f3f3f; int S,T,A[N],B[N],Enum,H[N],nxt[M],to[M],cap[M],len[M],cur[N];
bool vis[N];
LL Cost,dis[N]; inline int read()
{
int now=0;register char c=gc();
for(;!isdigit(c);c=gc());
for(;isdigit(c);now=now*10+c-48,c=gc());
return now;
}
inline void AE(int u,int v,int w,int c)
{
to[++Enum]=v, nxt[Enum]=H[u], H[u]=Enum, cap[Enum]=w, len[Enum]=c;
to[++Enum]=u, nxt[Enum]=H[v], H[v]=Enum, cap[Enum]=0, len[Enum]=-c;
}
bool SPFA()
{
static bool inq[N];
static std::queue<int> q;
memset(dis,0x3f,T+1<<3);//<<3!!!
dis[S]=0, q.push(S);
while(!q.empty())
{
int x=q.front();
q.pop(), inq[x]=0;
for(int i=H[x],v; i; i=nxt[i])
if(cap[i]&&dis[to[i]]>dis[x]+len[i])
dis[v=to[i]]=dis[x]+len[i], !inq[v]&&(q.push(v),inq[v]=1);
}
return dis[T]<INF;
}
bool DFS(int x)
{
if(x==T) return 1;
vis[x]=1;
for(int &i=cur[x],v; i; i=nxt[i])
if(!vis[v=to[i]]&&dis[v]==dis[x]+len[i]&&cap[i]&&DFS(v))//&&cap[i]!!! 有必要的啊
return --cap[i],++cap[i^1],Cost+=len[i],1;
return 0;
}
void MCMF()
{
while(SPFA())
{
memcpy(cur,H,T+1<<2), memset(vis,0,T+1);
while(DFS(S));
}
} int main()
{
static int e[N];
int n=read(),K=read(),t1=read(),t2=read(),P=n+1;
Enum=1, S=0, T=n+2; LL ans=0;
for(int i=1; i<=n; ++i) ans+=A[i]=read();
for(int i=1; i<=n; ++i) B[i]=read();
AE(S,P,K-t1,0);
for(int i=1; i<=K; ++i) AE(P,i,1,0);
for(int i=1,tmp=K-t1-t2; i<=n; ++i) AE(i,i+1>n?T:i+1,tmp,0), AE(i,i+K>n?T:i+K,1,A[i]-B[i]), e[i]=Enum;
MCMF(), printf("%lld\n",ans-Cost);
for(int i=1; i<=n; ++i) putchar(cap[e[i]]?'E':'S'); return 0;
}
上一篇:操作VCF卡片信息的第三方jar包:ez-vcard


下一篇:【BZOJ】1607: [Usaco2008 Dec]Patting Heads 轻拍牛头(特殊的技巧)