結果

問題 No.2310 [Cherry 5th Tune A] Against Regret
コンテスト
ユーザー kwm_t
提出日時 2023-05-20 13:43:46
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 778 ms / 6,000 ms
コード長 1,779 bytes
コンパイル時間 3,877 ms
コンパイル使用メモリ 261,792 KB
最終ジャッジ日時 2025-02-13 03:43:39
ジャッジサーバーID
（参考情報） 		judge5 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 30
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
#include <atcoder/all>
using namespace std;
using namespace atcoder;
//using mint = modint1000000007;
//const int mod = 1000000007;
using mint = modint998244353;
const int mod = 998244353;
//const int INF = 1e9;
//const long long LINF = 1e18;
#define rep(i, n) for (int i = 0; i < (n); ++i)
#define rep2(i,l,r)for(int i=(l);i<(r);++i)
#define rrep(i, n) for (int i = (n-1); i >= 0; --i)
#define rrep2(i,l,r)for(int i=(r-1);i>=(l);--i)
#define all(x) (x).begin(),(x).end()
#define allR(x) (x).rbegin(),(x).rend()
#define endl "\n"
#define P pair<int,int>
template<typename A, typename B> inline bool chmax(A & a, const B & b) { if (a < b) { a = b; return true; } return false; }
template<typename A, typename B> inline bool chmin(A & a, const B & b) { if (a > b) { a = b; return true; } return false; }

int main() {
	ios::sync_with_stdio(false);
	cin.tie(nullptr);
	int n; cin >> n;
	vector x(n + 1, vector<mint>(n + 1));
	rep(i, n + 1)rep(j, n + 1) {
		int _x; cin >> _x;
		x[i][j] = _x;
	}
	vector dp(n + 1, vector<mint>(n + 1));
	rep(i, n + 1) {
		dp[i][i] = 1;
		rep2(j, i, n + 1)rep2(k, j + 1, n + 1) dp[i][k] += dp[i][j] * x[j][k];
	}
	int q; cin >> q;
	while (q--) {
		int k; cin >> k;
		vector a(k, vector<int>(3));
		rep(i, k) rep(j, 3)cin >> a[i][j];
		sort(all(a));
		int st = 0;
		int ed = k * 2 + 1;
		int sz = ed + 1;
		vector nx(sz, vector<mint>(sz));
		nx[st][ed] = dp[0][n];
		rep(i, k) {
			nx[st][i * 2 + 1] += dp[0][a[i][0]];
			nx[i * 2 + 1][i * 2 + 2] += a[i][2];
			nx[i * 2 + 2][ed] += dp[a[i][1]][n];
		}
		rep(i, k)rep(j, k) nx[i * 2 + 2][j * 2 + 1] += dp[a[i][1]][a[j][0]];// *a[i][2];
		vector<mint>ans(sz);
		ans[0] = 1;
		rep(i, sz) rep(j, sz) ans[j] += ans[i] * nx[i][j];
		cout << ans.back().val() << endl;
	}
	return 0;
}
0