首页 > 其他分享 >AtCoder Regular Contest 132 F Takahashi The Strongest

AtCoder Regular Contest 132 F Takahashi The Strongest

时间:2023-05-24 13:15:16浏览次数:51  
标签:AtCoder typedef limits Contest sum long 132 otimes

洛谷传送门

AtCoder 传送门

没见过这种在新运算下做卷积的题,感觉挺新奇的。

考虑 Takahashi 成为绝对赢家的必要条件,发现前提是 Aoki 和 Snuke 出的要相同。不妨将每种策略映射到一个四进制数(\(P \to 1, R \to 2, S \to 3\)),定义运算 \(x \otimes y = \begin{cases} x & x = y \\ 0 & x \ne y \end{cases}\),设 \(a_i, b_i\) 分别为 Aoki 和 Snuke 每种策略的数量,那么实际上我们希望先求:

\[c_i = \sum\limits_{j \otimes k = i} a_j b_k \]

这个东西是可以 FWT 算的。先算 \(A_i = \sum\limits_{j \otimes i = i} a_j, B_i = \sum\limits_{j \otimes i = i} b_j\),那么 \(C_i = A_i B_i\),最后 \(C \to c\) IFWT 回去。定义一种新的集合从属关系 \(0 \subset 1, 0 \subset 2, 0 \subset 3\) 且 \(1,2,3\) 互不相交,那么 FWT 的过程相当于做一个超集和。

那么得到了 \(c_i\) 之后,我们还希望求 \(d_i = \sum\limits_{j \otimes i \ne 0} c_j\)。正着算不好算,考虑变成求 \(d_i = \sum\limits_{j \otimes i = 0} c_j\)。考虑设计一个 dp,\(f_{k,i}\) 表示考虑了前 \(k\) 位,\(i\) 表示前 \(k\) 位是已经确定的策略,后面的位和 \(d_i = \sum\limits_{j \otimes i = 0} c_j\) 中的 \(j\) 后面的位相同,这样的方案数。转移枚举在当前位填 \(1/2/3\),并且不能与 \(j\) 在这个位相同。

这样就做完了,时间复杂度 \(O(k 4^k)\)。

code
// Problem: F - Takahashi The Strongest
// Contest: AtCoder - AtCoder Regular Contest 132
// URL: https://atcoder.jp/contests/arc132/tasks/arc132_f
// Memory Limit: 1024 MB
// Time Limit: 5000 ms
// 
// Powered by CP Editor (https://cpeditor.org)

#include <bits/stdc++.h>
#define pb emplace_back
#define fst first
#define scd second
#define mems(a, x) memset((a), (x), sizeof(a))

using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef double db;
typedef long double ldb;
typedef pair<ll, ll> pii;

const int maxn = 15;
const int maxm = (1 << 24) + 50;

int n, m, K, all;
ll f[maxm], g[maxm], h[2][maxm];

inline void FWT(ll *f, ll op) {
	for (int k = 1; (k << 2) <= all; k <<= 2) {
		for (int i = 0; i < all; i += (k << 2)) {
			for (int j = 0; j < k; ++j) {
				f[i + j] += (f[i + j + k] + f[i + j + k * 2] + f[i + j + k * 3]) * op;
			}
		}
	}
}

void solve() {
	scanf("%d%d%d", &K, &n, &m);
	all = (1 << (K * 2));
	for (int _ = 0; _ < n; ++_) {
		char s[15];
		scanf("%s", s);
		int x = 0;
		for (int i = 0; i < K; ++i) {
			if (s[K - i - 1] == 'R') {
				x += (1 << (i * 2));
			} else if (s[K - i - 1] == 'S') {
				x += (2 << (i * 2));
			} else {
				x += (3 << (i * 2));
			}
		}
		++f[x];
	}
	for (int _ = 0; _ < m; ++_) {
		char s[15];
		scanf("%s", s);
		int x = 0;
		for (int i = 0; i < K; ++i) {
			if (s[K - i - 1] == 'R') {
				x += (1 << (i * 2));
			} else if (s[K - i - 1] == 'S') {
				x += (2 << (i * 2));
			} else {
				x += (3 << (i * 2));
			}
		}
		++g[x];
	}
	FWT(f, 1);
	FWT(g, 1);
	for (int i = 0; i < all; ++i) {
		f[i] *= g[i];
	}
	FWT(f, -1);
	for (int i = 0; i < all; ++i) {
		h[0][i] = f[i];
	}
	for (int i = 0, o = 0; i < K; ++i, o ^= 1) {
		mems(h[o ^ 1], 0);
		for (int S = 0; S < all; ++S) {
			if (!h[o][S]) {
				continue;
			}
			for (int j = 1; j <= 3; ++j) {
				int T = S - (((S >> (i * 2)) & 3) << (i * 2)) + (j << (i * 2));
				if (S != T) {
					h[o ^ 1][T] += h[o][S];
				}
			}
		}
	}
	for (int S = 0; S < all; ++S) {
		bool flag = 1;
		for (int i = 0; i < K; ++i) {
			if (((S >> (i * 2)) & 3) == 0) {
				flag = 0;
				break;
			}
		}
		if (flag) {
			printf("%lld\n", 1LL * n * m - h[K & 1][S]);
		}
	}
}

int main() {
	int T = 1;
	// scanf("%d", &T);
	while (T--) {
		solve();
	}
	return 0;
}

标签:AtCoder,typedef,limits,Contest,sum,long,132,otimes
From: https://www.cnblogs.com/zltzlt-blog/p/17427975.html

相关文章

  • AtCoder Beginner Contest 296
    AtCoderBeginnerContest296D题意给出n和m,问\(1\leqi,j\leqn\),使得\(ij\geqm\),求出这个乘积的最小值思路这两个乘数至少有一个在\([1,\sqrt{m}]\),枚举代码voidsolve(){ cin>>n>>m; intx=sqrt(m); if(n>=m){cout<<m<<endl;return;} if(x*x==m)......
  • AtCoder Regular Contest 139 C One Three Nine
    洛谷传送门AtCoder传送门闲话:做这场的B用时跟C差不多不会直接构造,因此这是一个无脑做法。考虑对于\(\forallx\in[1,n],y\in[1,m],(x+3y,3x+y)\)看成一个点,那么选择的\((x,y)\)要满足,不存在一行或一列有超过\(1\)个点。这启发我们对于合法的点\((a......
  • Atcoder 选做
    [ARC103F]DistanceSums(构造,重心)首先显然\(D_i\)的最小值被重心取到,不妨以重心为根。对于一条边连接的两个点\(x,y\),不妨设这条边\(x\)侧的点数为\(siz_x\),\(y\)侧为\(siz_y\)。那么\(D_y=D_x+siz_x-siz_y=D_x+siz_x-(n-siz_x)=D_x+2\timessiz_x-n\)。那么......
  • Atcoder Grand Contest 060 D - Same Descent Set
    先推式子。设\(f(S)\)表示decent集合恰好为\(S\)的排列个数,\(g(S)\)表示\(S\)是\(p\)的decent集合的一个子集的排列\(p\)个数,\(g'(\{a_1,a_2,\cdots,a_k\})=\dfrac{n!}{a_1!(a_2-a_1)!(a_3-a_2)!\cdots(a_k-a_{k-1})!(n-a_k)!}\),那么有:\[\begin{aligned}ans=&\......
  • Atcoder Beginner Contest ABC302 题解
    代码见此:https://atcoder.jp/contests/abc302/submissions?f.Task=&f.LanguageName=&f.Status=&f.User=frequenter。AAttackhttps://atcoder.jp/contests/abc302/tasks/abc302_a直接计算a/b,有余数的话答案加一。BFindSnukehttps://atcoder.jp/contests/abc302/tasks/abc......
  • AtCoder Regular Contest 132 D Between Two Binary Strings
    洛谷传送门AtCoder传送门提供一个dp思路。下文设串长为\(n\),串中\(1\)个数为\(m\)。考虑如何求\(d(s,t)\)。设\(s\)的\(1\)位置分别为\(a_1,a_2,...,a_m\),\(t\)的\(1\)位置分别为\(b_1,b_2,...,b_m\)。那么\(d(s,t)=\sum\limits_{i=1}^m|a_i-b......
  • 2023 Xian Jiaotong University Programming Contest
    A.大水题#include<bits/stdc++.h>#include<ext/rope>#include<ext/pb_ds/assoc_container.hpp>usingnamespacestd;usingnamespace__gnu_cxx;usingnamespace__gnu_pbds;#definefifirst#definesesecond#definelcu<<1#define......
  • AtCoder Beginner Contest 302 Ex Ball Collector
    洛谷传送门AtCoder传送门考虑如果只询问一次怎么做。连边\((a_i,b_i)\),对于每个连通块分别考虑。这是ARC111B,如果一个连通块是树,肯定有一个点不能被选;否则有环,一定能构造一种方案,使得每个点都被选。那么现在对于每个点都要求,考虑dfs时合并当前的\((a_u,b_u)\),并且使用......
  • 【题解】Atcoder ABC302 F,G,Ex
    完全不会G和Ex,这些套路还是要积累一下的。F.MergeSet题目描述:给定\(n\)个集合,每次可以合并两个有交的集合,问最少多少次合并可以让\(1\)和\(m\)位于同一个集合中。题目分析:一开始将题读成了将\([1,m]\)位于同一个集合中,然后就感觉这个题相当不可做。因为集合的元......
  • AtCoder Regular Contest 130 E Increasing Minimum
    这题太神仙了吧!感觉还不是很懂,但是尽力理一下思路。设\(t_x\)为最大的\(j\)使得\(i_j=x\),不存在则\(t_x=0\)。设\(1\simn\)的数按照\(t\)从小到大排序后是\(p_1,p_2,...,p_n\),设\(q_i\)为\(i\)在\(p\)中的排名,即\(q_{p_i}=i\)。发现正着不好考虑,......