CF1151F Sonya and Informatics 题解
Descirption
给定一个 01 序列,每次随机交换两个不同位置的数,求交换
对于
Solution
定义
定义
- 两端内部互换:
- 前面的
与后面的 互换: - 前面的
与后面的 互换: - 前面的
与后面的 互换: - 前面的
与后面的 互换:
整理一下:
然后矩阵快速幂求出
时间复杂度
Code
#include <bits/stdc++.h>
using namespace std;
namespace ModInt {
int inverse(int a, int m) {
int r = 1; while (a > 1) { int t = m / a; r = 1LL * r * (m - t) % m, a = m - t * a; } return r; }
#define MI ModInt::Mint<mod>
template <int mod> struct Mint {
typedef long long LL; int x;
Mint(int x = 0) : x(normal(x)) {}
Mint(LL x) : x(normal(x % mod)) {}
int normal(int x) { if (x >= mod) x -= mod; if (x < 0) x += mod; return x; }
int val() const { return x; }
Mint inv() const { assert(x != 0); return Mint(inverse(x, mod)); }
Mint operator - () const { return Mint(mod - x); }
Mint operator += (const Mint &rhs) { x = normal(x + rhs.x); return *this; }
Mint operator -= (const Mint &rhs) { x = normal(x - rhs.x); return *this; }
Mint operator *= (const Mint &rhs) { x = 1LL * x * rhs.x % mod; return *this; }
Mint operator /= (const Mint &rhs) { return *this *= rhs.inv(); }
friend Mint operator + (const Mint &lhs, const Mint &rhs) { Mint res = lhs; res += rhs; return res; }
friend Mint operator - (const Mint &lhs, const Mint &rhs) { Mint res = lhs; res -= rhs; return res; }
friend Mint operator * (const Mint &lhs, const Mint &rhs) { Mint res = lhs; res *= rhs; return res; }
friend Mint operator / (const Mint &lhs, const Mint &rhs) { Mint res = lhs; res /= rhs; return res; }
friend istream &operator >> (istream &is, Mint &a) { LL v; is >> v; a = Mint(v); return is; }
friend ostream &operator << (ostream &os, const Mint &a) { return os << a.x; }
};
}
namespace Milkcat {
typedef long long LL;
const int N = 105, K = 1e5 + 5, mod = 1e9 + 7;
LL n, k, qwq, cnt, a[N];
struct matrix {
MI a[N][N];
matrix() { memset(a, 0, sizeof a); }
matrix operator * (const matrix& x) const {
matrix c;
for (int i = 0; i <= cnt; i ++)
for (int k = 0; k <= cnt; k ++)
for (int j = 0; j <= cnt; j ++)
c.a[i][j] += a[i][k] * x.a[k][j];
return c;
}
} bas, mat;
MI qpow(MI b, LL k) { MI r = 1; for (; k; b *= b, k >>= 1) if (k & 1) r *= b; return r; }
matrix qpow(matrix b, LL k) { matrix r = b; k --; for (; k; b = b * b, k >>= 1) if (k & 1) r = r * b; return r; }
MI calc(MI x) { return x * (x - 1) / 2; }
int main() {
cin >> n >> k;
for (int i = 1; i <= n; i ++)
cin >> a[i], cnt += !a[i];
for (int i = 1; i <= cnt; i ++) qwq += a[i];
mat.a[qwq][0] = 1;
for (int i = 0; i <= cnt; i ++) {
bas.a[i][i] = calc(cnt) + calc(n - cnt) + (cnt - i) * i + i * (n - cnt - i);
if (i > 0) bas.a[i][i - 1] = (cnt - i + 1) * (n - cnt - i + 1);
if (i < cnt) bas.a[i][i + 1] = (i + 1) * (i + 1);
}
cout << (qpow(bas, k) * mat).a[0][0] / qpow(calc(n), k) << '\n';
return 0;
}
}
int main() {
freopen("randsort.in", "r", stdin);
freopen("randsort.out", "w", stdout);
int T = 1;
while (T --) Milkcat::main();
return 0;
}
/*
010->001->001
010->010->001
010->100->001
*/
【推荐】还在用 ECharts 开发大屏?试试这款永久免费的开源 BI 工具!
【推荐】国内首个AI IDE,深度理解中文开发场景,立即下载体验Trae
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· 历时 8 年,我冲上开源榜前 8 了!
· 分享一个我遇到过的“量子力学”级别的BUG。
· 物流快递公司核心技术能力-海量大数据处理技术
· 四大AI编程工具组合测评
· 关于能否用DeepSeek做危险的事情,DeepSeek本身给出了答案