【hihoCoder 1315】Reachable Permutations（状压DP）

Description

ZZX手里有一个数列 a[1..n]，它是1, 2, …, n 的一个排列。

现在ZZX想对这个数列进行一些变换：每一次他可以选择一对 i, j，满足1 ≤ i < j ≤ n 且 a[i] > a[j]，然后将 a[i] 和 a[j] 交换。

如果一个排列 b[1..n] 可以由初始数列 a 经过若干次变换而得到，那么ZZX就称b是可到达的。

现在ZZX想知道，有多少种不同的可到达排列。

由于答案可能很大，你只需要输出答案模109+7后的结果。

Solution

考虑如果数列A中全部是0/1，那么交换之后得到的序列B中所有1都在序列A中对应1的右边或者位置相同。
对于排列，我们考虑将大于等于 k k <script type="math/tex" id="MathJax-Element-17">k</script>的数看做1，否则看做0，然后状压DP即可。

Source

/************************************************
 * Au: Hany01
 * Prob: line
 * Email: hany01@foxmail.com
************************************************/

#include<bits/stdc++.h>

using namespace std;

typedef long long LL;
typedef pair<int, int> PII;
#define rep(i, j) for (register int i = 0, i##_end_ = (j); i < i##_end_; ++ i)
#define For(i, j, k) for (register int i = (j), i##_end_ = (k); i <= i##_end_; ++ i)
#define Fordown(i, j, k) for (register int i = (j), i##_end_ = (k); i >= i##_end_; -- i)
#define Set(a, b) memset(a, b, sizeof(a))
#define Cpy(a, b) memcpy(a, b, sizeof(a))
#define fir first
#define sec second
#define pb(a) push_back(a)
#define mp(a, b) make_pair(a, b)
#define ALL(a) (a).begin(), (a).end()
#define SZ(a) ((int)(a).size())
#define INF (0x3f3f3f3f)
#define INF1 (2139062143)
#define Mod (1000000007)
#define debug(...) fprintf(stderr, __VA_ARGS__)
#define y1 wozenmezhemecaia

template <typename T> inline bool chkmax(T &a, T b) { return a < b ? a = b, 1 : 0; }
template <typename T> inline bool chkmin(T &a, T b) { return b < a ? a = b, 1 : 0; }

inline int read()
{
    register int _, __; register char c_;
    for (_ = 0, __ = 1, c_ = getchar(); c_ < '0' || c_ > '9'; c_ = getchar()) if (c_ == '-') __ = -1;
    for ( ; c_ >= '0' && c_ <= '9'; c_ = getchar()) _ = (_ << 1) + (_ << 3) + (c_ ^ 48);
    return _ * __;
}

inline void File()
{
#ifndef ONLINE_JUDGE
    freopen("line.in", "r", stdin);
    freopen("line.out", "w", stdout);
#endif
}

const int maxn = 22, maxk = (1 << 21);

int n, a[maxn], dp[maxk];

inline void check(int x)
{
    int cnt = __builtin_popcount(x), lmt = n - cnt + 1, t = 0, cur1 = 0, cur2 = 0;
    For(i, 1, n) if (a[i] >= lmt) t |= (1 << (i - 1));
    For(T, 1, cnt) {
        while (!(x & (1 << cur1))) ++ cur1;
        while (!(t & (1 << cur2))) ++ cur2;
        if (cur1 < cur2) return ;
        ++ cur1, ++ cur2;
    }
    rep(i, n) if (!(x & (1 << i))) (dp[x | (1 << i)] += dp[x]) %= Mod;
}

int main()
{
    File();

    n = read();
    For(i, 1, n) a[i] = read();

    dp[0] = 1;
    rep(i, 1 << n) if (dp[i]) check(i);

    printf("%d\n", dp[(1 << n) - 1]);

    return 0;
}
//客心已百念，孤游重千里。
//    -- 何逊《相送》