bzoj3661
网络流/贪心
网络流做法是对于每一列，如果一个兔子下一天继续可以存在，那么连一条容量为1的边，然后设立一个中转站，来控制可以换的数量，容量限制l。时限100s，能跑过去我的太慢了，一个点100s
正解是贪心。不会证明。每一只兔子，我们当然希望能够尽量连续多跑几天，给其他兔子留机会。每次把在的兔子和不在的兔子排序，具体程序，然后看前l个能不能交换。
巨慢无比的网络流
#include<bits/stdc++.h> using namespace std; const int N = 805 * 805, inf = 1 << 29; struct edge { int nxt, to, f; } e[N * 40]; int n, m, k, l, S1, T, S, cnt = 1; int q[N * 3 + 1621], d[N * 3 + 1621], head[N * 3 + 1621], iter[N * 3 + 1621]; int a[810][810], b[810][810]; #define id(i, j) (i - 1) * m + j inline void link(int u, int v, int f) { e[++cnt].nxt = head[u]; head[u] = cnt; e[cnt].to = v; e[cnt].f = f; } inline void ins(int u, int v, int f) { link(u, v, f); link(v, u, 0); } bool bfs() { int l = 1, r = 0; q[++r] = 0; memset(d, 0, sizeof(d)); d[0] = 1; while(l <= r) { int u = q[l++]; for(int i = head[u]; i; i = e[i].nxt) if(e[i].f && !d[e[i].to]) { d[e[i].to] = d[u] + 1; if(e[i].to == T) return 1; q[++r] = e[i].to; } } return 0; } int dfs(int u, int delta) { if(u == T) return delta; int ret = 0; for(int &i = iter[u]; i && delta; i = e[i].nxt) if(e[i].f && d[e[i].to] == d[u] + 1) { int x = dfs(e[i].to, min(delta, e[i].f)); e[i].f -= x; e[i ^ 1].f += x; ret += x; delta -= x; } d[u] = 0; return ret; } inline int dinic() { int ret = 0; while(bfs()) { for(int i = 0; i <= T; ++i) iter[i] = head[i]; ret += dfs(0, inf); } return ret; } void build() { T = n * m + N * 2 + m * 2 + 2; S1 = n * m + N * 2 + m * 2 + 1; ins(0, S1, k); for(int i = 1; i <= n; ++i) for(int j = 1; j <= m; ++j) if(a[i][j]) { ins(id(i, j), id(i, j) + N, 1); b[i][j] = cnt - 1; } for(int i = 1; i <= n; ++i) { if(a[i][1]) ins(S1, id(i, 1), 1); if(a[i][m]) ins(id(i, m) + N, T, 1); } for(int j = 2; j <= m; ++j) ins(n * m + N * 2 + j, n * m + N * 2 + j + m, l); for(int i = 1; i <= n; ++i) for(int j = 1; j <= m; ++j) if(a[i][j]) { if(j < m) if(a[i][j + 1]) ins(id(i, j) + N, id(i, j + 1), 1); if(j < m) ins(id(i, j) + N, n * m + N * 2 + j + 1, 1); if(j > 1) ins(n * m + N * 2 + j + m, id(i, j), 1); } } int main() { // freopen("input.txt", "r", stdin); scanf("%d%d%d%d", &n, &m, &k, &l); for(int i = 1; i <= n; ++i) { char s[810]; scanf("%s", s + 1); for(int j = 1; j <= m; ++j) a[i][j] = s[j] - '0'; } build(); int ans = dinic(); if(ans < k) { puts("-1"); return 0; } for(int j = 1; j <= m; ++j) { for(int i = 1; i <= n; ++i) if(b[i][j]) if(e[b[i][j]].f == 0) printf("%d ", i); puts(""); } // fclose(stdin); return 0; }
贪心
#include<bits/stdc++.h> using namespace std; const int N = 810 * 810, inf = 1 << 29; int n, m, k, l; typedef pair<int, int> PII; int pos[810][810], a[810][810], mark[810]; PII b[810][810], c[810][810]; bool cp1(PII x, PII y) { return x.first < y.first; } bool cp2(PII x, PII y) { return x.first > y.first; } int main() { scanf("%d%d%d%d", &n, &m, &k, &l); for(int i = 1; i <= n; ++i) { char s[810]; scanf("%s", s + 1); for(int j = 1; j <= m; ++j) a[i][j] = s[j] - '0'; } for(int j = 1; j <= m; ++j) for(int i = 1; i <= n; ++i) if(a[i][j]) { int Pos = j; b[j][i] = c[j][i] = make_pair(1, i); while(a[i][Pos + 1]) { ++Pos; ++b[j][i].first; ++c[j][i].first; } } sort(b[1] + 1, b[1] + n + 1, cp2); for(int i = 1; i <= k; ++i) if(b[1][i].first > 0) pos[1][++pos[1][0]] = b[1][i].second; if(pos[1][0] < k) { puts("-1"); return 0; } sort(pos[1] + 1, pos[1] + k + 1); for(int i = 1; i <= k; ++i) printf("%d ", pos[1][i]); puts(""); for(int j = 2; j <= m; ++j) { memset(mark, 0, sizeof(mark)); for(int i = 1; i <= k; ++i) mark[pos[j - 1][i]] = 1; for(int i = 1; i <= n; ++i) if(mark[i]) c[j][i].first = -1; else b[j][i].first = inf; sort(b[j] + 1, b[j] + n + 1, cp1); sort(c[j] + 1, c[j] + n + 1, cp2); for(int i = 1; i <= l; ++i) if(b[j][i].first < c[j][i].first) pos[j][i] = c[j][i].second; else pos[j][i] = b[j][i].second; for(int i = l + 1; i <= k; ++i) pos[j][i] = b[j][i].second; sort(pos[j] + 1, pos[j] + k + 1); for(int i = 1; i <= k; ++i) printf("%d ", pos[j][i]); puts(""); } return 0; }
