1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
|
#include <bits/stdc++.h>
inline char read() {
static const int IN_LEN = 1000000;
static char buf[IN_LEN], *s, *t;
if (s == t) {
t = (s = buf) + fread(buf, 1, IN_LEN, stdin);
if (s == t) return -1;
}
return *s++;
}
template<class T>
inline void read(T &x) {
static bool iosig;
static char c;
for (iosig = false, c = read(); !isdigit(c); c = read()) {
if (c == '-') iosig = true;
if (c == -1) return;
}
for (x = 0; isdigit(c); c = read())
x = (x + (x << 2) << 1) + (c ^ '0');
if (iosig) x = -x;
}
const int OUT_LEN = 1000000;
char obuf[OUT_LEN], *oh = obuf;
inline void print(char c) {
if (oh == obuf + OUT_LEN) fwrite(obuf, 1, OUT_LEN, stdout), oh = obuf;
*oh++ = c;
}
template<class T>
inline void print(T x) {
static int buf[30], cnt;
if (x == 0) {
print('0');
} else {
if (x < 0) print('-'), x = -x;
for (cnt = 0; x; x /= 10) buf[++cnt] = x % 10 + 48;
while (cnt) print((char)buf[cnt--]);
}
}
inline void flush() {
fwrite(obuf, 1, oh - obuf, stdout);
}
#define long long long
const int MAXN = 1010;
struct Node {
int v, f, w, index;
Node(int v, int f, int w, int index) : v(v), f(f), w(w), index(index) {}
};
std::vector<Node> edge[MAXN];
inline void addEdge(int u, int v, int f, int w) {
edge[u].push_back(Node(v, f, w, edge[v].size()));
edge[v].push_back(Node(u, 0, -w, edge[u].size() - 1));
}
typedef std::pair<int, long> Pair;
const long INF = 9187201950435737471;
inline Pair minCostMaxFlow(int s, int t, int n) {
Pair ans(0, 0);
while (true) {
static bool vis[MAXN];
static long dis[MAXN];
static int prev[MAXN], pree[MAXN];
std::queue<int> q;
memset(dis, 127, sizeof(long) * (n + 1));
memset(vis, 0, sizeof(bool) * (n + 1));
q.push(s), dis[s] = 0;
while (!q.empty()) {
register int u = q.front();
q.pop();
vis[u] = false;
for (register int i = 0; i < edge[u].size(); i++) {
Node *e = &edge[u][i];
if (e->f && dis[u] + e->w < dis[e->v]) {
dis[e->v] = dis[u] + e->w;
prev[e->v] = u, pree[e->v] = i;
if (!vis[e->v]) q.push(e->v), vis[e->v] = true;
}
}
}
if (dis[t] == INF) break;
register int flow = INT_MAX;
for (register int i = t; i != s; i = prev[i])
flow = std::min(flow, edge[prev[i]][pree[i]].f);
ans.first += flow, ans.second += flow * dis[t];
for (register int i = t; i != s; i = prev[i]) {
Node *e = &edge[prev[i]][pree[i]];
e->f -= flow, edge[e->v][e->index].f += flow;
}
}
return ans;
}
int main() {
#ifndef ONLINE_JUDGE
freopen("in.in", "r", stdin);
#endif
register int n, k, ms, me;
read(n), read(k), read(ms), read(me);
register int S = 0, T = n + 1;
register long ans = 0;
static int s[MAXN], e[MAXN];
for (register int i = 1; i <= n; i++) read(s[i]), ans += s[i];
for (register int i = 1; i <= n; i++) read(e[i]);
addEdge(S, n - k + 2, k - ms, 0), addEdge(1, T, k - ms, 0);
for (register int i = 2, r = n - k + 2; i <= r; i++)
addEdge(i, i - 1, k - ms - me, 0);
for (register int i = 1; i <= n; i++)
addEdge(std::min(i + 1, n - k + 2), std::max(1, i - k + 1), 1, s[i] - e[i]);
print(ans - minCostMaxFlow(S, T, T + 1).second);
flush();
return 0;
}
|
