【组合数学】【P4996】 咕咕咕

Description

小 F 注意到,自己总是在某些情况下会产生歉意。每当他要检查自己的任务表来决定下一项任务的时候,如果当前他干了某些事情,但是没干另一些事情,那么他就会产生一定量的歉意——比如,无论他今天看没看比赛,只要没有补完月赛的锅,他都会在选择任务的时候产生 1 点歉意。小 F 完成所有任务后,他这一天的歉意值等于他每次选择任务时的歉意之和。

过高的歉意值让小 F 感到不安。现在,小 F 告诉你他还有 n 项任务,并告诉你在 m 种情况中的一种 statei 的情况下,小 F 会产生 ai 点歉意。请你帮忙计算一下,小 F 在那一天所有可能的完成所有任务方式的歉意值之和是多少。

由于答案可能很大,你只需要输出答案对 998244353 取模即可。

Input

输入一行两个整数 n,m,表示有 n 项任务,在 m 种情况中下小 F 会产生歉意值。

输入接下来 m 行,每行有一个长度为 n0/1 串 和 statei 一个歉意值 aistatei,j0/1 表示第 j 项任务此时没做 / 已经做了。

详情请参考样例和样例解释。

Output

输出一行一个整数,表示小 F 在那一天所有可能的完成任务方式的歉意值之和对 998244353 取模的结果。

Hint

1  n  20

1  m  min(2n,105)

1  ai  105

Solution

考虑状压,发现只能枚举子集转移,复杂度 O(3n),显然过不去

接着发现对答案产生贡献的状态并不多,考虑是否可以直接统计这些状态对答案产生的贡献。

考虑一个非常显然的事情,如果状态 A1 的个数与状态 B1 的个数相同,那么从空集转到状态 AB 的方案数是相同的,同时从这两个状态转到全集的方案数也是相同的。

根据乘法原理,从空集经过状态 S 再到全集的方案数为 cnts1 × cnts0,其中 cntx 为将 x0 变成 1 的方案数

于是答案即为 i=1m cntstatei1 × cntstatei0 × ai

考虑怎么求 cnt 数组。显然可以 O(3n) 状压在本机跑一会然后打表。

考虑正常一点的做法。发现 cnt 的组合意义即为从几个数中选出几个当 1 然后剩下的 0 一次性变成 1 的方案数

于是 cnti = j=1icntij × Cij

于是递推一下组合数,递推一下cnt就做完了。

Code

#include <cstdio>
#ifdef ONLINE_JUDGE
#define freopen(a, b, c)
#endif
#define rg register
#define ci const int
#define cl const long long

typedef long long int ll;

namespace IPT {
	const int L = 1000000;
	char buf[L], *front=buf, *end=buf;
	char GetChar() {
		if (front == end) {
			end = buf + fread(front = buf, 1, L, stdin);
			if (front == end) return -1;
		}
		return *(front++);
	}
}

template <typename T>
inline void qr(T &x) {
	rg char ch = IPT::GetChar(), lst = ' ';
	while ((ch > '9') || (ch < '0')) lst = ch, ch=IPT::GetChar();
	while ((ch >= '0') && (ch <= '9')) x = (x << 1) + (x << 3) + (ch ^ 48), ch = IPT::GetChar();
	if (lst == '-') x = -x;
}

template <typename T>
inline void ReadDb(T &x) {
	rg char ch = IPT::GetChar(), lst = ' ';
	while ((ch > '9') || (ch < '0')) lst = ch, ch = IPT::GetChar();
	while ((ch >= '0') && (ch <= '9')) x = x * 10 + (ch ^ 48), ch = IPT::GetChar();
	if (ch == '.') {
		ch = IPT::GetChar();
		double base = 1;
		while ((ch >= '0') && (ch <= '9')) x += (ch ^ 48) * ((base *= 0.1)), ch = IPT::GetChar();
	}
	if (lst == '-') x = -x;
}

namespace OPT {
	char buf[120];
}

template <typename T>
inline void qw(T x, const char aft, const bool pt) {
	if (x < 0) {x = -x, putchar('-');}
	rg int top=0;
	do {OPT::buf[++top] = x % 10 + '0';} while (x /= 10);
	while (top) putchar(OPT::buf[top--]);
	if (pt) putchar(aft);
}

const int maxt = 25;
const int MOD = 998244353;

int n, m, ans;
int C[maxt][maxt], frog[maxt];

void beginning();

int main() {
	freopen("1.in", "r", stdin);
	qr(n); qr(m);
	beginning();
	while (m--) {
		char ch = '-'; int _cnt = 0;
		do {ch = IPT::GetChar();} while ((ch != '0') && (ch != '1'));
		do {_cnt += ch - '0', ch = IPT::GetChar();} while ((ch == '0') || (ch == '1'));
		int _a = 0; qr(_a);
		ans = (ans + 1ll * frog[_cnt] * _a % MOD * frog[n - _cnt]) % MOD;
	}
	qw(ans, '\n', true);
	return 0;
}

void beginning() {
	C[0][0] = 1;
	for (rg int i = 1; i <= n; ++i) {
		C[i][0] = 1;
		for (rg int j = 1; j <= i; ++j) C[i][j] = (C[i - 1][j - 1] + C[i - 1][j]) % MOD;
	}
	frog[0] = 1;
	for (rg int i = 1; i <= n; ++i) {
		for (rg int j = 0; j <= i; ++j) {
			frog[i] = (frog[i] + 1ll * frog[i - j] * C[i][j]) % MOD;
		}
	}
}

Summary

当对答案产生贡献的状态很少时,考虑直接计算每个状态的贡献。

posted @   一扶苏一  阅读(279)  评论(0编辑  收藏  举报
编辑推荐:
· 从 HTTP 原因短语缺失研究 HTTP/2 和 HTTP/3 的设计差异
· AI与.NET技术实操系列:向量存储与相似性搜索在 .NET 中的实现
· 基于Microsoft.Extensions.AI核心库实现RAG应用
· Linux系列:如何用heaptrack跟踪.NET程序的非托管内存泄露
· 开发者必知的日志记录最佳实践
阅读排行:
· winform 绘制太阳,地球,月球 运作规律
· TypeScript + Deepseek 打造卜卦网站:技术与玄学的结合
· AI 智能体引爆开源社区「GitHub 热点速览」
· Manus的开源复刻OpenManus初探
· 写一个简单的SQL生成工具
点击右上角即可分享
微信分享提示