General
 
 
# Author Problem Lang Verdict Time Memory Sent Judged  
252391210 Practice:
DisconnectedGraph
1781F - 9 C++20 (GCC 13-64) Accepted 1747 ms 20020 KB 2024-03-20 07:40:08 2024-03-20 07:40:08
→ Source
#pragma GCC optimize("O2")
#include <bits/stdc++.h>
using namespace std;
#define _rep(i_,a_,b_) for(int i_ = (a_); i_ <= (b_); ++i_)
#define mid ((L+R) >> 1)
#define multiCase() int testCnt = in(); _rep(curCase,1,testCnt)
#ifdef ONLINE_JUDGE
#define debug(...) 0
#else
#define debug(...) fprintf(stderr, __VA_ARGS__), fflush(stderr)
#endif
using ll = long long;
using pii = pair<int,int>;

int in(void) { int x; scanf("%d", &x); return x; } ll inl(void) { ll x; scanf("%lld", &x); return x; }
void out(int x) { printf("%d ", x); } void outln(int x) { printf("%d\n", x); }
void out(ll x) { printf("%lld ", x); } void outln(ll x) { printf("%lld\n", x); }
template<typename T> void chkmax(T &a, const T &b) { a = max(a, b); } 
template<typename T> void chkmin(T &a, const T &b) { a = min(a, b); } 
const int p = 998244353, kN = 550, N = 500;
int fpm(int a, int b) {
	int r = 1;
	for(;b;b >>= 1, a = 1ll * a * a % p) if(b & 1) r = 1ll * r * a % p;
	return r;
}
int n = in(), q = 1ll * in() * fpm(10000, p - 2) % p;
int f[kN][2][4][kN], C[kN][kN], sum[kN][2][4][kN];
int main() { 
	C[0][0] = 1;
	_rep(i,1,N) {
		C[i][0] = C[i][i] = 1;
		_rep(j,1,i - 1) C[i][j] = (C[i - 1][j] + C[i - 1][j - 1]) % p; 
	}
	f[0][0][3][0] = 1;
	_rep(i,0,3)
		f[1][0][i][0] = q,
		f[1][1][i][1] = (1 - q + p) % p;
	auto init_sum = [](int l) {
		_rep(a,0,1) _rep(b,l <= 1 ? 0 : 1,3) _rep(j,0,n + 1) 
			sum[l][a][b][j] = (f[l][a][b][j] + (j == 0 ? 0 : sum[l][a][b][j - 1])) % p;
	};
	init_sum(0);
	init_sum(1);
	// debug("f %d\n", sum[0][0][3][0]);
	// debug("%d %d\n", f[1][0][2][0], sum[1][0][1][0]);
	_rep(i,2,n) {
		_rep(l,0,i - 1) {
			//首先考虑 x(x)x 形式
			_rep(nd,0,i)	
				f[i][0][1][nd] = (f[i][0][1][nd] + 
				1ll * f[i - l][0][0][nd] * (sum[l][0][3][nd] + sum[l][1][3][nd]) % p * C[i - 1][l]) % p,
				nd && (
				f[i][0][1][nd] = (f[i][0][1][nd] + 
				1ll * (f[l][0][3][nd] + f[l][1][3][nd]) * sum[i - l][0][0][nd - 1] % p * C[i - 1][l]) % p
				, 1);
		}
		_rep(l,0,i - 1) {
			_rep(nd,0,i)	
				f[i][0][2][nd] = (f[i][0][2][nd] + 
				1ll * f[i - l][0][1][nd] * (sum[l][0][3][nd + 1] + sum[l][1][3][nd + 1]) % p * C[i - 1][l]) % p,
				nd && (
				f[i][0][2][nd] = (f[i][0][2][nd] + 
				1ll * (f[l][0][3][nd + 1] + f[l][1][3][nd + 1]) * sum[i - l][0][1][nd - 1] % p * C[i - 1][l]) % p
				, 1);
		}
		_rep(l,0,i - 1) {
			_rep(nd,0,i)	
				f[i][0][3][nd] = (f[i][0][3][nd] + 
				1ll * f[i - l][0][2][nd] * (sum[l][0][3][nd] + sum[l][1][3][nd]) % p * C[i - 1][l]) % p,
				nd && (
				f[i][0][3][nd] = (f[i][0][3][nd] + 
				1ll * (f[l][0][3][nd] + f[l][1][3][nd]) * sum[i - l][0][2][nd - 1] % p * C[i - 1][l]) % p
				, 1);
		}
		_rep(l,0,i - 1) {
			//然后考虑 x)x(x 形式
			_rep(nd,0,i)	
				f[i][1][1][nd] = (f[i][1][1][nd] + 
				1ll * f[i - l][1][0][nd] * (sum[l][0][3][nd] + sum[l][1][3][nd]) % p * C[i - 1][l]) % p,
				nd && (
				f[i][1][1][nd] = (f[i][1][1][nd] + 
				1ll * (f[l][0][3][nd] + f[l][1][3][nd]) * sum[i - l][1][0][nd - 1] % p * C[i - 1][l]) % p
				, 1);
		}
		_rep(l,0,i - 1) {
			_rep(nd,0,i)	
				nd && (
				f[i][1][2][nd] = (f[i][1][2][nd] + 
				1ll * f[i - l][1][1][nd] * (sum[l][0][3][nd - 1] + sum[l][1][3][nd - 1]) % p * C[i - 1][l]) % p
				, 1),
				nd && (
				f[i][1][2][nd] = (f[i][1][2][nd] + 
				1ll * (f[l][0][3][nd - 1] + f[l][1][3][nd - 1]) * sum[i - l][1][1][nd - 1] % p * C[i - 1][l]) % p
				, 1);
		}
		_rep(l,0,i - 1) {
			_rep(nd,0,i)	
				f[i][1][3][nd] = (f[i][1][3][nd] + 
				1ll * f[i - l][1][2][nd] * (sum[l][0][3][nd] + sum[l][1][3][nd]) % p * C[i - 1][l]) % p,
				nd && (
				f[i][1][3][nd] = (f[i][1][3][nd] + 
				1ll * (f[l][0][3][nd] + f[l][1][3][nd]) * sum[i - l][1][2][nd - 1] % p * C[i - 1][l]) % p
				, 1);
		}
		init_sum(i);
		// debug("%d\n", f[i][0][0][0]);
		// debug("%d\n", f[i][0][1][0]);
		// debug("%d\n", f[i][0][2][0]);
		// debug("%d\n", f[i][0][3][0]);
	}
	// debug("%d %d\n", f[n][0][3][0], f[n][1][3][0]);
	int multiple = 1, len = 1;
	_rep(i,1,n) multiple = 1ll * multiple * fpm(len, p - 2) % p, len += 2;
	outln(1ll * (f[n][0][3][0] + f[n][1][3][0]) * multiple % p);
	return 0;
}

/* 
a list of keywords
clear empty push_back pop_back push pop top front back
emplace_back emplace push_front pop_front insert erase
find count set reset bitset map vector string multiset
first second iterator prev next deque multimap reverse
sort begin end list modify query init check calc prime
putchar getchar puts scanf printf max min swap replace
make_pair make_tuple numeric_limits auto function null
*/
?
Time: ? ms, memory: ? KB
Verdict: ?
Input
?
Participant's output
?
Jury's answer
?
Checker comment
?
Diagnostics
?
Click to see test details