Solution AT987

2022-04-03 19:00:27 阅读：162 来源： 互联网

标签：dbinom limits int sum Solution AT987 aligned mod

\(\mathtt{description}\)：

Link
\(t\) 组询问，每组给出 \(n,k\le 10^5\)，求 \(\begin{aligned}\sum\limits_{i=0}^{k}\dbinom{n}{i}\end{aligned}\)，
\(t\le 10^5\)。

\(\mathtt{Solution}\)：

\(n\to n+1\) 时：

\[\begin{aligned}\sum\limits_{i=0}^{k}\dbinom{n+1}{i}&=\sum\limits_{i=0}^{k}\begin{pmatrix}\dbinom{n}{i}+\dbinom{n}{i-1}\end{pmatrix}\\&=2\sum\limits_{i=0}^{k}\dbinom{n}{i}-\dbinom{n}{k}\end{aligned} \]

移项反解得 \(n\to n-1\) 时：

\[\begin{aligned}\sum\limits_{i=0}^{k}\dbinom{n-1}{i}=\dfrac{\sum\limits_{i=0}^{k}\dbinom{n}{i}+\dbinom{n}{k}}{2}\end{aligned} \]

\(k\to k+1\) 时：

\[\begin{aligned}\sum\limits_{i=0}^{k+1}\dbinom{n}{i}=\sum\limits_{i=0}^{k}\dbinom{n}{i}+\dbinom{n}{k+1}\end{aligned} \]

同理，\(k\to k-1\) 时:

\[\begin{aligned}\sum\limits_{i=0}^{k-1}\dbinom{n}{i}=\sum\limits_{i=0}^{k}\dbinom{n}{i}-\dbinom{n}{k+1}\end{aligned} \]

注意到上面四个转移式的余项都为单个组合数，而预处理阶乘及逆元计算组合数时 \(O(1)\) 的，所以上述四个转移也是 \(O(1)\) 的。

莫队即可，复杂度 \(O(n\sqrt n)\)。

#include <bits/stdc++.h>
#define int long long
namespace mystd {
	inline int read() {
	    char c = getchar();
	    int x = 0, f = 1;
	    while (c < '0' || c > '9') f = (c == '-') ? -1 : f, c = getchar();
	    while (c >= '0' && c <= '9') x = (x << 1) + (x << 3) + c - '0', c = getchar();
	    return x * f;
	}
	inline void write(int x) {
	    if (x < 0) x = ~(x - 1), putchar('-');
	    if (x > 9) write(x / 10);
	    putchar(x % 10 + '0');
	}
}
using namespace std;
using namespace mystd;

const int maxn = 2e5 + 200;// 不知道为啥要开大的空间，不然就过不去样例
const int mod = 1e9 + 7;
int n, l, r, b, c, mx, ans[maxn], fac[maxn], inv[maxn], ifac[maxn];
struct qry {
	int l, r, id;
	bool operator < (const qry &rhs) const {
		return l / b == rhs.l / b ? r < rhs.r : l < rhs.l;
	}
} q[maxn];

int C(int n, int m) { return fac[n] * ifac[m] % mod * ifac[n - m] % mod; }
void addr() { c = ((c * 2 % mod - C(r++, l)) % mod + mod) % mod; }
void delr() { c = inv[2] * ((c + C(--r, l)) % mod) % mod; }
void addl() { c = ((c - C(r, l--)) % mod + mod) % mod; }
void dell() { c = (c + C(r, ++l)) % mod; }

signed main() {
	n = read(), b = sqrt(n); 
	for (int i = 1; i <= n; i++) {
		q[i] = (qry) { read(), read(), i };
		swap(q[i].l, q[i].r);
		mx = max(mx, q[i].r);
	}
	sort(q + 1, q + n + 1);
	r = 1, c = 1;
	fac[0] = inv[1] = ifac[0] = 1;
	for (int i = 1; i <= mx; i++) {
		if (i >= 2) inv[i] = (mod - mod / i) * inv[mod % i] % mod;
		fac[i] = fac[i - 1] * i % mod;
		ifac[i] = ifac[i - 1] * inv[i] % mod;
	}
	for (int i = 1; i <= n; i++) {
		while (l > q[i].l) addl();
		while (l < q[i].l) dell();
		while (r < q[i].r) addr();
		while (r > q[i].r) delr();
		ans[q[i].id] = c;
	}
	for (int i = 1; i <= n; i++) write(ans[i]), puts("");
	return 0;
}

标签：dbinom,limits,int,sum,Solution,AT987,aligned,mod
来源： https://www.cnblogs.com/Ender32k/p/16096965.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

ICode9

Solution AT987