当前位置: 首页 > news >正文

NOIP20251009E

\(2n\) 个糖果,第 \(i\) 个糖果的美味值为 \(a_i\)

\(n\) 天中,每天选择恰好两个糖果一起吃,一天的开心度为两个糖果美味值之和。

\(f(i)\) 为所有 \(\frac{(2n)!}{n!2^n}\) 种不同方案中,第 \(i\) 小的开心度之和。

\(i = 1 \sim n\),输出 \(f(i) \pmod 998244353\)

\(1 \le n \le 100,1 \le a_i \le 10^9\)


称一种吃糖果的方案为一个匹配,例如对于 \(n = 2\) 的所有匹配为 \(S_1 = \{(a_1, a_2), (a_3,a_4)\}, S_2 = \{(a_1,a_3),(a_2,a_4)\},S_3=\{(a_1,a_4),(a_2,a_3)\}\)

\(b_i\) 为某个匹配中和是第 \(i\) 小的二元组的和,例如对于 \(\{(1, 3), (2, 6)\}\)\(b_1 = 1 + 3\)

\(f(i)\) 为所有匹配的 \(b_i\) 的和,即 \(f(i) = \sum\limits_{b_i}b_i\)

拆贡献,\(b_i = \sum\limits_{x=1}^V[b_i\ge x]\)\(b_i \ge x\) 等价于这个匹配至少有 \(n-x+1\)\(b\) 满足 \(b \ge x\)

枚举 \(x\)\(b\)

\(F_i\) 为满足 \(i = \sum\limits_{j=1}^n[b_j\ge x]\) 的匹配的个数,则 \(\sum\limits_{j=i}^nF_j\) 可以贡献到 \(b_i\)

直接求 \(F\) 不好求。设 \(G_i\) 为钦定了 \(i\)\(b\) 满足 \(b \ge x\) 的方案数,则有

\[G_i = \sum\limits_{j=i}^n\binom{j}{i}F_j \]

根据二项式反演可得

\[F_i = \sum\limits_{j=i}^n(-1)^{j-i}\binom{j}{i}G_j \]

考虑求 \(G\),将 \(a\) 升序排序。

\(a_i\)\(a_j(j<i)\) 可以被钦定,当且仅当 \(a_i + a_j \ge x\)\(a_j \le a_i\),即 \(x - a_i \le a_j \le a_i\),设对于 \(i\) 满足这个限制的 \(j\) 的取值范围为 \([l_i,r_i]\)。因为 \(a_i \le a_{i+1}\),所以 \([l_i,r_i]\subseteq[l_{i+1},r_{i+1}]\)

\(f_{i,j}\) 为考虑了 \(i\) 个糖果,有 \(j\) 对二元组满足其和 \(\ge x\) 的方案数。

选择 \(a_i\) 与前面的一个糖果:\(f_{i,j} += f_{i-1,j-1} \times (r_i - l_i + 1 - 2(j-1))\)

不选 \(a_i\)\(f_{i,j} += f_{i-1,j}\)

\(G_i = f_{2n,i} \times \frac{(2n - 2i)!}{(n - i)!2^{n-i}}\),表示先选 \(i\) 个二元组,剩下的随便。

以上可以在 \(O(n^2V)\) 的时间复杂度内求出每个 \(f(i)\)

但是 \(V \le 10^9\),考虑 \(b\) 的值,因为 \(b = a_i + a_j\),则值不同的 \(b\) 只有 \(O(n^2)\) 个,将这些值作为 \(x\),再乘上出现的次数就行了,时间复杂度 \(O(n^4)\)

#define M 998244353ll
inline int MOD(ll x) { return x >= M ? x - M : x; }
inline void ADD(ll& x, ll y) { x = MOD(x + y); }bool _st;const int N = 202;
int cp, n, n2;
ll a[N], c[N * N], C[N][N], fac[N], ifac[N], pw[N], ipw[N];int l[N];
ll f[N][N], F[N], G[N], ans[N];
void solve(int x) {rep(i, 1, n2) l[i] = ::std::lower_bound(a + 1, a + 1 + n2, c[x] - a[i]) - a;memset(f, 0, sizeof(f));f[0][0] = 1;rep(i, 1, n2) rep(j, 0, n){if((j << 1) > i) break;ADD(f[i][j], f[i - 1][j]);int val = i - l[i] - 2 * (j - 1);if(j && val > 0)ADD(f[i][j], f[i - 1][j - 1] * val % M);}rep(i, 1, n) G[i] = ((f[n2][i] * fac[n2 - 2 * i]) % M * (ifac[n - i] * ipw[n - i] % M)) % M;rep(i, 1, n){F[i] = 0;rep(j, i, n)if((i - j) & 1) ADD(F[i], M - C[j][i] * G[j] % M);else ADD(F[i], C[j][i] * G[j] % M);}per(i, 1, n) ADD(F[i], F[i + 1]);rep(i, 1, n) ADD(ans[i], F[n - i + 1] * (c[x] - c[x - 1]) % M);
}bool _ed;
int main() { FIN ::std::cerr << (&_ed - &_st) / 1024.0 / 1024 << ::std::endl;Rep(i, 0, N){C[i][0] = 1;rep(j, 1, i) C[i][j] = MOD(C[i - 1][j] + C[i - 1][j - 1]);}fac[0] = pw[0] = 1;Rep(i, 1, N)fac[i] = fac[i - 1] * i % M,pw[i] = MOD(pw[i - 1] + pw[i - 1]);ifac[N - 1] = cksm<M>(fac[N - 1], M - 2);ipw[N - 1] = cksm<M>(pw[N - 1], M - 2);per(i, 0, N - 2)ifac[i] = ifac[i + 1] * (i + 1) % M,ipw[i] = MOD(ipw[i + 1] + ipw[i + 1]);qr(n);n2 = n << 1;rep(i, 1, n2) qr(a[i]);::std::sort(a + 1, a + 1 + n2);rep(i, 1, n2) rep(j, i + 1, n2) c[++cp] = a[i] + a[j];::std::sort(c + 1, c + 1 + cp);cp = ::std::unique(c + 1, c + 1 + cp) - c - 1;rep(i, 1, cp) solve(i);rep(i, 1, n) qw(ans[i]), ps;return 0;
}
http://www.hskmm.com/?act=detail&tid=28905

相关文章:

  • 2025七水硫酸锌订做厂家推荐榜:品质保证与客户信赖之选
  • 2025螺杆泵厂家最新推荐榜:高效稳定与优质服务的行业首选!
  • 2025南通婚纱摄影最新推荐榜:创意拍摄与贴心服务的完美结合
  • 语义slam - MKT
  • 尝试茶叶数据集
  • 2025氧化镁供应厂家推荐:松辽镁业高纯度优质选择!
  • 2025硅藻土订制厂家口碑推荐:品质卓越与专业服务的双重保障
  • 20251011
  • 2025数控滚齿机源头厂家推荐榜:高精度与高效能的首选!
  • (第四次)回归与决策树
  • 2025数控滚齿机订做厂家推荐:吉莱特智能装备,精准高效品质
  • 2025机械加工实力厂家推荐:鑫铭机械专业制造,品质卓越首选
  • 高考语文做法
  • 2025机械加工优质厂家推荐榜:技术精湛与高效服务的行业先锋
  • P10960 SUBSTRACT 个人题解
  • 牛客网刷题
  • 2025新型千斤顶厂家推荐:柳州市联桥科技,品质卓越服务到位
  • 2025深圳网站建设推荐:华企网络专业定制,助力企业线上腾飞
  • 2025石头纸设备批发厂家推荐鼎浩包装,环保高效生产首选!
  • 2025液压阀块供货厂家最新推荐榜:品质卓越与高效服务的行业
  • 2025年PP鱼池优质厂家推荐:超众渔业机械,环保耐用首选!
  • centos安装atop工具,检测服务器情况
  • 完整教程:MongoDB Ops Manager部署
  • 2025医疗器械微弧氧化优质厂家推荐,华源漆业技术领先服务到
  • 2025气柱袋优质厂家推荐:戈尔德包装,防护包装解决方案专家
  • 【网络协议】SSL与TLS的关系 - 教程
  • 2025深圳网站建设公司最新推荐榜:创意设计与专业服务引领者
  • 2025年安全光栅厂家最新推荐榜:精准防护与高效性能的工业首
  • 2025七水硫酸锌实力厂家推荐:安通环保科技,品质卓越信赖之
  • 20232306刘博2025-2026-1《网络与系统攻防技术》实验一实验报告