4842: [Neerc2016]Delight for a Cat_1283: 序列
题目大意:ls是一个特别堕落的小朋友,对于n个连续的小时,他将要么睡觉要么打隔膜,一个小时内他不能既睡觉也打隔膜,因此一个小时内他只能选择睡觉或者打隔膜,当然他也必须选择睡觉或打隔膜,对于每一个小时,他选择睡觉或打隔膜的愉悦值是不同的,对于第i个小时,睡觉的愉悦值为si,打隔膜的愉悦值为ei,同时又有一个奥妙重重的规定:对于任意一段连续的k小时,ls必须至少有t1时间在睡觉,t2时间在打隔膜。那么ls想让他获得的愉悦值尽量大,他该如何选择呢?
数据范围:n,k(1<=k<=n<=1000),t1,t2(0<=t1,t2<=k;t1+t2<=k)
题解:
发现每段区间都必须满足$0$的个数有一个上限,$1$的个数有一个上限。
把那个关于上限的不等式变成等式之后,单纯形即可。
我真的是这么以为的......
妈的卡常。
学了一下怎么用费用流维护这个东西,还挺好用的。
代码:
#include <bits/stdc++.h>
#define N 2010
#define M 100010
using namespace std;
typedef long long ll;
ll ans, all;
int n, k, t1, t2, tot, S, T;
int inq[N], to[M], nxt[M], flow[M], pe[N], pv[N], head[N], pos[N], A[N], B[N];
ll cost[M], dis[N];
queue<int> q;
char *p1, *p2, buf[100000];
#define nc() (p1 == p2 && (p2 = (p1 = buf) + fread(buf, 1, 100000, stdin), p1 == p2) ? EOF : *p1 ++ )
int rd() {
int x = 0, f = 1;
char c = nc();
while (c < 48) {
if (c == '-')
f = -1;
c = nc();
}
while (c > 47) {
x = (((x << 2) + x) << 1) + (c ^ 48), c = nc();
}
return x * f;
}
inline void add(int a, int b, int c, int d) {
to[tot] = b;
cost[tot] = c;
flow[tot] = d;
nxt[tot] = head[a];
head[a] = tot ++ ;
to[tot] = a;
cost[tot] = -c;
flow[tot] = 0;
nxt[tot] = head[b];
head[b] = tot ++ ;
}
bool bfs() {
memset(dis, 0x80, sizeof dis);
q.push(S);
dis[S] = 0;
int u, i;
while(!q.empty()) {
u = q.front();
q.pop();
inq[u] = 0;
for(i = head[u]; i != -1; i = nxt[i]) {
if(dis[to[i]] < dis[u] + cost[i] && flow[i]) {
dis[to[i]] = dis[u] + cost[i];
pv[to[i]] = u;
pe[to[i]] = i;
if(!inq[to[i]])
inq[to[i]] = 1, q.push(to[i]);
}
}
}
return dis[T]>(ll)0x8080808080808080ll;
}
int main() {
n = rd(), k = rd(), t1 = rd(), t2 = rd();
S = 0, T = ((n - k + 1) << 1) + 2;
memset(head, -1, sizeof head);
int i;
for (i = 1; i <= n; i ++ )
A[i] = rd();
for (i = 1; i <= n; i ++ )
B[i] = rd(), all += B[i];
for (i = 1; i <= n; i ++ ) {
pos[i] = tot + 1;
add(max(1, ((i - k) << 1) + 1), min((i << 1) + 1, ((n - k + 1) << 1) + 1), A[i] - B[i], 1);
}
add(S, 1, 0, t1);
for (i = 1; i <= n - k + 1; i ++ ) {
add(i << 1, (i << 1) - 1, 0, 1 << 30);
add(i << 1, (i << 1) + 1, 0, 1 << 30);
add(S, i << 1, 0, k - t1 - t2);
if(i != n - k + 1)
add((i << 1) + 1, T, 0, k - t1 - t2);
}
add(((n - k + 1) << 1) + 1, T, 0, k - t2);
while (bfs()) {
int mf = 1 << 30;
for (i = T; i != S; i = pv[i])
mf = min(mf, flow[pe[i]]);
ans += dis[T] * mf;
for (i = T; i != S; i = pv[i])
flow[pe[i]] -= mf, flow[pe[i] ^ 1] += mf;
}
printf("%lld\n", all + ans);
for (i = 1; i <= n; i ++ ) {
if(flow[pos[i]])
printf("S");
else
printf("E");
}
return 0;
}
小结:跟志愿者招募多倍经验,没放在题目。