1.タイトルリンク:
HDU-3811
2.トピックの主なアイデア:
最初の行の2つの整数n、m。
次のm行では、各行に2つの整数aiとbiが含まれており、番号aiがbiであることを示しています。
mの制限の少なくとも1つを満たす1からnまでの順列の数を見つけます。
3.分析:
最も暴力的なアプローチは、バイナリ列挙によって満たされるべき制限を列挙し、許容範囲の概念を使用して計算することですが、これは間違いなくTLEです...
答えの反対のイベントを考えてみましょう。1〜nの配置は、mの制限のいずれも満たしていません。
なお、dp [i] [j]は、iの数が配置されていることを意味し、その数の使用状態は、mの制限を満たさないスキームの数です。
詳細については、コードを参照してください。
4.コードの実装:
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int M = (int)17;
const int inf = 0x3f3f3f3f;
const ll mod = (ll)1e9 + 7;
ll fac[M + 5];
int one[1<<M];
bool ban[M][M];
ll dp[2][1<<M];
int cal(int n)
{
int cnt = 0;
while(n)
{
cnt += (n & 1);
n >>= 1;
}
return cnt;
}
int main()
{
fac[0] = 1; for(int i = 1; i <= 17; ++i) fac[i] = fac[i - 1] * i;
for(int i = 0; i < (1<<17); ++i) one[i] = cal(i);
int T; scanf("%d", &T);
for(int ca = 1; ca <= T; ++ca)
{
memset(ban, 0, sizeof(ban));
int n, m; scanf("%d %d", &n, &m);
for(int i = 0, a, b; i < m; ++i) scanf("%d %d", &a, &b), ban[--a][--b] = 1;
for(int i = 0; i < n; ++i) dp[0 & 1][1<<i] = !ban[0][i];
for(int i = 1; i < n; ++i)
{
for(int j = 0; j < (1<<n); ++j)
{
if(one[j] != i + 1) continue;
dp[i & 1][j] = 0;
for(int k = 0; k < n; ++k)
{
if(!(j & (1<<k)) || ban[i][k]) continue;
dp[i & 1][j] += dp[i - 1 & 1][j^(1<<k)];
}
}
}
printf("Case %d: %lld\n", ca, fac[n] - dp[n - 1 & 1][(1<<n) - 1]);
}
return 0;
}