[ZJOI2016]小星星

嘟嘟嘟


這題剛開始把\(17!\)想成了\(2 ^ {17}\)，覺得全排暴力枚舉就行，而後寫一半發現好像只能過\(n \leqslant 10\)的點。
討論正解以前，咱們先想狀壓dp，畢竟這個數據範圍就像狀壓。\(dp[i][j][S]\)表示點\(i\)所在子樹中，\(i\)對應\(j\)，子樹的對應狀況爲\(S\)時的方案數。轉移的時候枚舉子樹\(v\)和對應編號\(k\)，而後由於編號不能重，咱們枚舉\(S\)的補集的子集，而後暴力合併兩個集合。
用子集dp的話，複雜度能達到\(O(n ^ 3 3 ^ {n})\)，仍然會超時。


而後題解就說，若是不考慮編號不能重的限制，咱們就能夠省去第三維了。
如今加上了限制，就能夠用容斥。
考慮若是有一些點的編號被用了好多遍，那必定說明有一些點沒被用到。那咱們乾脆直接刪掉這些點而後再dp。因而\(O(2 ^ n)\)枚舉哪些點被刪掉，而後根據個數奇偶性容斥一下便可。
至關於咱們要求的是圖中剛好有\(n\)個點都被映射的方案數=至少有\(n\)個點被映射的方案數-至少有\(n - 1\)個點被映射的方案數+至少有\(n - 2\)個點被映射-至少有\(n - 3\)個點……
總的來講，咱們容斥的是圖上被刪去的點，即哪些點沒有被映射，而不是樹上哪些點不參與dp。（這裏我糾結了很久）

#include<cstdio>
#include<iostream>
#include<cmath>
#include<algorithm>
#include<cstring>
#include<cstdlib>
#include<cctype>
#include<vector>
#include<stack>
#include<queue>
#include<assert.h>
using namespace std;
#define enter puts("") 
#define space putchar(' ')
#define Mem(a, x) memset(a, x, sizeof(a))
#define In inline
typedef long long ll;
typedef double db;
const int INF = 0x3f3f3f3f;
const db eps = 1e-8;
const int maxn = 20;
In ll read()
{
  ll ans = 0;
  char ch = getchar(), last = ' ';
  while(!isdigit(ch)) last = ch, ch = getchar();
  while(isdigit(ch)) ans = (ans << 1) + (ans << 3) + ch - '0', ch = getchar();
  if(last == '-') ans = -ans;
  return ans;
}
In void write(ll x)
{
  if(x < 0) x = -x, putchar('-');
  if(x >= 10) write(x / 10);
  putchar(x % 10 + '0');
}
In void MYFILE()
{
#ifndef mrclr
  freopen("star3.in", "r", stdin);
  freopen("ha.out", "w", stdout);
#endif
}

int n, m, G[maxn][maxn];
struct Edge
{
  int nxt, to;
}e[maxn << 1];
int head[maxn], ecnt = -1;
In void addEdge(int x, int y)
{
  e[++ecnt] = (Edge){head[x], y};
  head[x] = ecnt;
}

int pos[maxn];
In bool dfs0(int now, int _f)
{
  for(int i = head[now], v; ~i; i = e[i].nxt)
    {
      if((v = e[i].to) == _f) continue;
      if(!G[pos[now]][pos[v]]) return 0;
      if(!dfs0(v, now)) return 0;
    }
  return 1;
}

In void work0()
{
  for(int i = 1; i <= n; ++i) pos[i] = i;
  int ans = 0;
  do
    {
      ans += dfs0(1, 0);
    }while(next_permutation(pos + 1, pos + n + 1));
  write(ans), enter;
}

bool vis[maxn];
ll dp[maxn][maxn], ans;
In void dfs(int now, int _f)
{
  for(int i = 1; i <= n; ++i) dp[now][i] = 1;
  for(int i = head[now], v; ~i; i = e[i].nxt)
    {
      if((v = e[i].to) == _f) continue;
      dfs(v, now);
      for(int j = 1; j <= n; ++j)
	{
	  ll sum = 0;
	  for(int k = 1; k <= n; ++k)
	    sum += dp[v][k] * (G[j][k] & vis[j] & vis[k]);
	  dp[now][j] *= sum;
	}
    }
}

int main()
{
  //MYFILE();
  Mem(head, -1);
  n = read(), m = read();
  for(int i = 1; i <= m; ++i)
    {
      int x = read(), y = read();
      G[x][y] = G[y][x] = 1;
    }
  for(int i = 1; i < n; ++i)
    {
      int x = read(), y = read();
      addEdge(x, y), addEdge(y, x);
    }
  if(n <= 10) {work0(); return 0;}  //開局寫了個n!暴力……
  for(int S = 1; S < (1 << n); ++S)
    {
      int tot = n;
      for(int i = 1; i <= n; ++i)
	if((S >> (i - 1)) & 1) vis[i] = 1, --tot;
	else vis[i] = 0;
      dfs(1, 0);
      ll tp = 0;
      for(int i = 1; i <= n; ++i) tp += dp[1][i];
      if(tot & 1) ans -= tp;
      else ans += tp;
    }
  write(ans), enter;
  return 0;
}