洛谷P3355 骑士共存问题

原题链接:传送门

这题的标签是

由于还没学网络流，我能想到的就是去用匈牙利算法去求最大独立集

 看这个图片，我们如果骑士在黄色格子那么它能攻击到的就都红色格子，反之就是骑士在红色格子那么他能攻击到的就都是黄色格子，那么我们可以根据此性质来把能互相攻击的格子建边。易证由于是骑士共存问题，所以两端点不能同时选取，有连边的两点必然属于二分图的两端，所以如果两个端点且有边不能同时被选取，因而必然能构成一个匹配。根据性质没有匹配到的点都可以被选取（反证法：没有匹配说明没有连边，若有连边则可以增加匹配数进而与最大匹配矛盾)。
最后问题转化为求最大独立集=n-最小点覆盖(最小点覆盖=最大匹配数)

首先建边问题我们可以选择把二维变一维地去处理，其他的就很朴素了。

AC代码：

#include <bits/stdc++.h>
#define x first
#define y second
using namespace std;
typedef pair<int, int> PII;
const int N = 505;
int n, m, k;
bool g[N][N], vis[N * N];
int cnt[N][N], match[N * N];
int h[N * N], e[4 * N * N], ne[4 * N * N], idx;

void add(int a, int b)
{
  e[idx] = b, ne[idx] = h[a], h[a] = idx++;
}

const int dx[8] = {-2, -1, 1, 2, 2, 1, -1, -2};
const int dy[8] = {1, 2, 2, 1, -1, -2, -2, -1};

inline int read()
{
  char ch = getchar();
  int s = 0, w = 1;
  while (ch < '0' || ch > '9')
  {
    if (ch == '-')
      w = -1;
    ch = getchar();
  }
  while (ch >= '0' && ch <= '9')
  {
    s = s * 10 + ch - '0';
    ch = getchar();
  }
  return s * w;
}

bool find(int x)
{
  for (int i = h[x]; i != -1; i = ne[i])
  {
    int j = e[i];
    if (!vis[j])
    {
      vis[j] = 1;
      if (!match[j] || find(match[j]))
      {
        match[j] = x;
        return true;
      }
    }
  }
  return false;
}

int main()
{
  n = read(), k = read();
  memset(h, -1, sizeof h);
  for (int i = 0; i < k; i++)
  {
    g[read()][read()] = 1;
  }
  int tt = 0;
  for (int i = 1; i <= n; i++)
  {
    for (int j = 1; j <= n; j++)
    {
      cnt[i][j] = ++tt;
    }
  }
  for (int i = 1; i <= n; i++)
  {
    for (int j = 1; j <= n; j++)
    {
      if (g[i][j]) continue;
      for (int k = 0; k < 8; k++)
      {
        int a = i + dx[k], b = j + dy[k];
        if (a < 1 || a > n || b < 1 || b > n) continue;
        if (g[a][b]) continue;
        add(cnt[i][j], cnt[a][b]), add(cnt[a][b], cnt[i][j]);
      }
    }
  }
  int ans = n * n - k;
  int res = 0;
  for (int i = 1; i <= tt; i++)
  {
    memset(vis, 0, sizeof vis);
    if (find(i)) res++;
  }
  printf("%d\n", ans - res / 2);
  return 0;
}