[CF997C]Sky Full of Stars_二项式反演_等比数列_容斥原理
Sky Full of Stars
题目链接:http://codeforces.com/problemset/problem/997/C
数据范围:略。
题解:
首先考虑拟对象,如果至少有一行完全相等即可。
这个的答案就需要多步容斥:$\sum\limits_{i = 1} ^ n (-1)^{i + 1}\cdot 3 ^ i\cdot 3 ^ {n \cdot (n - i)}$。
那么至少有一列的答案跟这个一样。
把他俩加一起就是答案么?我们需要减去什么?
显然,需要减掉至少有一行且至少有一列的。
这个怎么弄?
是这样的,如果我们钦定了$i$行$j$列都必须相等之后,$x$行$y$列相等这种情况会被算$C_{x} ^ {i}\times C_{y} ^ {j}$次。
假设,$F_{(i,j)}$表示$i$行$j$列都相等的答案。
假设$A_{(i,j)}$表示$F_{(i,j)}$的容斥系数,发现当$A_{(i,j)} = (-1)^{i + j + 1}$时满足题意。
故此,至少有一列切至少有一行的答案是:
$\sum\limits_{i = 1} ^ n\sum\limits_{j = 1} ^ n (-1) ^ {i + j + 1} C_{n} ^ {i}C_{n}^{j} \cdot 3\cdot 3^{(n - i)\times (n - j)}$。
这个怎么看都是$O(n^2)$的对不对....
我们考虑把他转化转化:
$=\sum\limits_{i = 1} ^ n\sum\limits_{j = 1} ^ n (-1) ^ {i + j + 1} C_{n} ^ {i}C_{n}^{j} \cdot 3^{(n - i)\times (n - j) + 1}$
$=\sum\limits_{i = 1} ^ n\sum\limits_{j = 1} ^ n (-1) ^ {i + j + 1} C_{n} ^ {i}C_{n}^{j} \cdot 3^{n^2 - in - jn +ij + 1}$
$=\sum\limits_{i = 1} ^ n\sum\limits_{j = 1} ^ n (-1) ^ {i + j + 1} C_{n} ^ {i}C_{n}^{j} \cdot 3^{n^2}\cdot 3^{-in}\cdot 3^{-jn}\cdot 3^{ij}$
$=-3^{n^2}\cdot \sum\limits_{i = 1} ^ n (-1) ^ i\cdot C_{n} ^ {i}\cdot 3^{-in} \sum\limits_{j = 1} ^ n (-1) ^ j C_{n}^{j} \cdot 3^{-jn}\cdot 3^{ij}$
$=-3^{n^2}\cdot \sum\limits_{i = 1} ^ n (-1) ^ i\cdot C_{n} ^ {i}\cdot 3^{-in} \sum\limits_{j = 1} ^ n (-1) ^ j C_{n}^{j} \cdot (3^{-n})^j\cdot (3^i)^j$
$=-3^{n^2}\cdot \sum\limits_{i = 1} ^ n (-1) ^ i\cdot C_{n} ^ {i}\cdot 3^{-in} \sum\limits_{j = 1} ^ n C_{n}^{j} \cdot (-3^{i - n})^j$
发现,第二个$\sum$后面的式子类似于二项式反演,$(-3^{i - n})$相当于$(a-b)^n$中的$b$。
故此式子可以被化简成:
$=-3^{n^2}\cdot \sum\limits_{i = 1} ^ n (-1) ^ i\cdot C_{n} ^ {i}\cdot 3^{-in} \cdot (1 - (-3^{i - n}))^n$。
这个就能$O(n)$求了。
代码:
#include <bits/stdc++.h> #define setIO(s) freopen(s".in", "r", stdin), freopen(s".out", "w", stdout) #define N 1000010 using namespace std; typedef long long ll; const int mod = 998244353 ; int fac[N], inv[N]; int qpow(int x, int y) {
int ans = 1;
while (y) {
if (y & 1) {
ans = (ll)ans * x % mod;
}
y >>= 1;
x = (ll)x * x % mod;
}
return ans;
} inline int C(int x, int y) {
return (ll)fac[x] * inv[y] % mod * inv[x - y] % mod;
} int main() {
// setIO("c");
int n;
cin >> n ;
fac[0] = inv[0] = 1;
for (int i = 1; i <= n; i ++ ) {
fac[i] = (ll)fac[i - 1] * i % mod;
inv[i] = qpow(fac[i], mod - 2);
} int ans = 0;
for (int i = 1; i <= n; i ++ ) {
ans = (ans + (ll)qpow(3, ((ll)n * (n - i) % (mod - 1) + i) % (mod - 1))
* qpow(mod - 1, i + 1) % mod
* C(n, i) % mod) % mod;
}
ans = ans * 2 % mod; int mdl = 0;
for (int i = 0; i < n; i ++ ) {
int t = (mod - qpow(3, i)) % mod;
mdl = (mdl + (ll)C(n, i)
* qpow(mod - 1, i + 1) % mod
* ( ( (ll) qpow(t + 1, n) + mod - qpow(t, n) ) % mod) % mod) % mod;
}
ans = (ans + (ll)mdl * 3) % mod;
cout << ans << endl ;
fclose(stdin), fclose(stdout);
return 0;
}
最新文章
- JavaScript——之对象参数的引用传递
- Symbol not found for architecture arm64 错误
- spring多线程与定时任务
- sublime好看的字体
- cdoj 65 CD Making 水题
- MetadataType的使用,MVC的Model层数据验证
- 字符串分割与数组的分割 split()VSsplice()&;slice()
- UVA 1605 Building for UN
- linux进程解析--进程切换
- 老毛桃pe安装系统
- Java的JDK下StringBuffer与StringBuilder的区别
- H3C ER6300 + 两台 H3C S5120 组网举例
- _编程语言_C语言_C++_时间
- date 类型转为varchar
- Linux下chkconfig命令详解转载
- Codeforces 1091 Good Bye 2018
- [HihoCoder1596]Beautiful Sequence
- 图像边缘检測--OpenCV之cvCanny函数
- 如何用Python实现常见机器学习算法-3
- Go语言学习笔记三: 常量
热门文章
- Codeforces Round #564 比赛总结
- 工作流调度系统Azkaban的简介和使用
- [题解] [BZOJ4152] The Captain
- 解决ScrollView中Recyclerview显示不全,滑动不流畅的问题
- 走进JavaWeb技术世界4:Servlet 工作原理详解
- 图解如何利用Intellij IDEA进行代码重构
- EntityFramework Core Code First 已有数据库
- AngularJS开发中常用的写法,如:获取URL参数、路由跳转、$http、获取元素等
- intellj idea 常用快捷键
- Smarty的分页实现