usaco-camp

Short, concise solutions for problems from the USACO camp

  1. 1
  2. 2
  3. 3
  4. 4
  5. 5
  6. 6
  7. 7
  8. 8
  9. 9
  10. 10
  11. 11
  12. 12
  13. 13
  14. 14
  15. 15
  16. 16
  17. 17
  18. 18
  19. 19
  20. 20
  21. 21
  22. 22
  23. 23
  24. 24
  25. 25
  26. 26
  27. 27
  28. 28
  29. 29
  30. 30
  31. 31
  32. 32
  33. 33
  34. 34
  35. 35
  36. 36
  37. 37
  38. 38
  39. 39
  40. 40
  41. 41
  42. 42
  43. 43
  44. 44
  45. 45
  46. 46
  47. 47
  48. 48
  49. 49
  50. 50
  51. 51
  52. 52
  53. 53
  54. 54
  55. 55
  56. 56
  57. 57
  58. 58
  59. 59
  60. 60
  61. 61
  62. 62
  63. 63
  64. 64
  65. 65
  66. 66
  67. 67
  68. 68
  69. 69
  70. 70
  71. 71
  72. 72
  73. 73
  74. 74
  75. 75
  76. 76
  77. 77
  78. 78
  79. 79
  80. 80
  81. 81
  82. 82
  83. 83
  84. 84
  85. 85
  86. 86
  87. 87
  88. 88
  89. 89
  90. 90
  91. 91
#include <bits/stdc++.h>
using namespace std;
using ll = long long;
constexpr ll INF = 1e12;
 
int dx[] = { 2, 1, -1, -2, -2, -1, 1, 2 };
int dy[] = { 1, 2, 2, 1, -1, -2, -2, -1 };
char A[101][101];
 
template<int V, class T = ll> class max_flow {
	struct edge {
		int t, rev;
		T cap, f;
	};
	vector<edge> adj[V];
	int dist[V];
	int ptr[V];
	bool bfs(int s, int t) {
		memset(dist, -1, sizeof dist);
		dist[s] = 0;
		queue<int> q({ s });
		while (!q.empty() && dist[t] == -1) {
			int n = q.front();
			q.pop();
			for (auto& e : adj[n]) {
				if (dist[e.t] == -1 && e.cap != e.f) {
					dist[e.t] = dist[n] + 1;
					q.push(e.t);
				}
			}
		}
		return dist[t] != -1;
	}
	T augment(int n, T amt, int t) {
		if (n == t) return amt;
		for (; ptr[n] < adj[n].size(); ptr[n]++) {
			edge& e = adj[n][ptr[n]];
			if (dist[e.t] == dist[n] + 1 && e.cap != e.f) {
				T flow = augment(e.t, min(amt, e.cap - e.f), t);
				if (flow != 0) {
					e.f += flow;
					adj[e.t][e.rev].f -= flow;
					return flow;
				}
			}
		}
		return 0;
	}
public:
	void add(int u, int v, T cap = 1, T rcap=0) {
		adj[u].push_back({ v, (int) adj[v].size(), cap, 0 });
		adj[v].push_back({ u, (int) adj[u].size() - 1, rcap, 0 });
	}
	T calc(int s, int t) {
		T flow = 0;
		while (bfs(s, t)) {
			memset(ptr, 0, sizeof ptr);
			while (T df = augment(s, INF, t)) flow += df;
		}
		return flow;
	}
	void clear() {
		for (int n = 0; n < V; n++) adj[n].clear();
	}
};
 
int main() {
	int R, C; 
	cin >> R >> C;
	int ans = R * C;
	for (int i = 0; i < R; ++i) {
		for (int j = 0; j < C; ++j) cin >> A[i][j];
	}
 
	max_flow<101*101> MF;
	for (int i = 0; i < R; ++i) {
		for (int j = 0; j < C; ++j) {
			if (A[i][j] == '#') { --ans; continue; }
			if (i + j & 1) {
				MF.add(10000, i * C + j, (A[i][j] == 'C' ? INF : 1));
				for (int d = 0; d < 8; ++d) {
					int x = i + dx[d], y = j + dy[d];
					if (x < 0 || x >= R || y < 0 || y >= C || A[x][y] == '#') continue;
					MF.add(i * C + j, x * C + y, INF);
				}
			}
			else MF.add(i * C + j, 10001, (A[i][j] == 'C' ? INF : 1));
		}
	}
	cout << ans - MF.calc(10000, 10001) << '\n';
}