BZOJ5323：[JXOI2018]游戏

最新推荐文章于 2020-04-13 15:13:36 发布

Cyhlnj

最新推荐文章于 2020-04-13 15:13:36 发布

阅读量319

点赞数 2

CC 4.0 BY-SA版权

分类专栏：就是数学(计数等)

本文链接：https://blog.youkuaiyun.com/oi_Konnyaku/article/details/87775412

就是数学(计数等) 专栏收录该内容

29 篇文章

订阅专栏

本文探讨了埃氏筛法在筛选质数中的高效应用，并结合组合数学原理，解决了一个特定的算法问题。通过枚举和组合的方式，详细介绍了如何计算特定范围内数字的选择方案数量。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

传送门
不难发现，所有不能被其他数筛掉的数是一定要选的，只有选了这些数字才能结束
假设有 $m$ 个，枚举结束时间 $x$ ，答案就是 $\sum \binom{x-1}{m-1}m!(n-m)!x$
埃氏筛法即可求出 $m$

# include <bits/stdc++.h>
using namespace std;
typedef long long ll;

const int maxn(1e7 + 5);
const int mod(1e9 + 7);

inline void Inc(int &x, const int y) {
    x = x + y >= mod ? x + y - mod : x + y;
}

int l, r, n, m, ans, fac[maxn], inv[maxn];
bitset <maxn> vis;

int main() {
	int i, j, k;
	scanf("%d%d", &l, &r), n = r - l + 1;
	if (l == 1) {
		for (k = i = 1; i < n; ++i) k = (ll)k * i % mod;
		ans = (ll)n * (n + 1) / 2 % mod * k % mod;
		return printf("%d\n", ans), 0;
	}
	for (i = l; i <= r; ++i)
		if (!vis[i]) for (++m, j = i; j <= r; j += i) vis[j] = 1;
	inv[0] = inv[1] = fac[0] = fac[1] = 1;
	for (i = 2; i <= n; ++i) inv[i] = (ll)(mod - mod / i) * inv[mod % i] % mod;
	for (i = 2; i <= n; ++i) fac[i] = (ll)fac[i - 1] * i % mod, inv[i] = (ll)inv[i] * inv[i - 1] % mod;
	for (i = m; i <= n; ++i) Inc(ans, (ll)fac[i] * inv[i - m] % mod);
	ans = (ll)ans * m % mod * fac[n - m] % mod, printf("%d\n", ans);
	return 0;
}