結果

問題 No.1750 ラムドスウイルスの感染拡大-hard
ユーザー MasKoaTSMasKoaTS
提出日時 2021-11-20 18:11:49
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 377 ms / 2,000 ms
コード長 1,651 bytes
コンパイル時間 4,231 ms
コンパイル使用メモリ 263,640 KB
実行使用メモリ 4,380 KB
最終ジャッジ日時 2023-09-02 14:04:26
合計ジャッジ時間 9,761 ms
ジャッジサーバーID
（参考情報） 		judge12 / judge13
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
4,380 KB
testcase_01 AC 2 ms
4,376 KB
testcase_02 AC 1 ms
4,380 KB
testcase_03 AC 1 ms
4,380 KB
testcase_04 AC 10 ms
4,380 KB
testcase_05 AC 1 ms
4,376 KB
testcase_06 AC 1 ms
4,376 KB
testcase_07 AC 1 ms
4,380 KB
testcase_08 AC 66 ms
4,376 KB
testcase_09 AC 67 ms
4,380 KB
testcase_10 AC 67 ms
4,376 KB
testcase_11 AC 51 ms
4,376 KB
testcase_12 AC 71 ms
4,380 KB
testcase_13 AC 66 ms
4,376 KB
testcase_14 AC 349 ms
4,376 KB
testcase_15 AC 353 ms
4,376 KB
testcase_16 AC 360 ms
4,380 KB
testcase_17 AC 365 ms
4,380 KB
testcase_18 AC 377 ms
4,376 KB
testcase_19 AC 364 ms
4,380 KB
testcase_20 AC 242 ms
4,376 KB
testcase_21 AC 282 ms
4,380 KB
testcase_22 AC 50 ms
4,376 KB
testcase_23 AC 346 ms
4,376 KB
testcase_24 AC 38 ms
4,376 KB
testcase_25 AC 91 ms
4,376 KB
testcase_26 AC 29 ms
4,380 KB
testcase_27 AC 2 ms
4,376 KB
testcase_28 AC 3 ms
4,376 KB
testcase_29 AC 2 ms
4,380 KB
testcase_30 AC 48 ms
4,380 KB
testcase_31 AC 44 ms
4,376 KB
testcase_32 AC 44 ms
4,376 KB
testcase_33 AC 43 ms
4,376 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <math.h>
#include <bits/stdc++.h>
#include <atcoder/all>
using namespace std;
using namespace atcoder;
#define rep(i, l, n) for (int i = (l); i < (n); i++)
#define max(p, q) ((p) > (q) ? (p) : (q))
#define min(p, q) ((p) < (q) ? (p) : (q))
#define all(x) x.begin(), x.end()
#define fi first
#define se second
#define lb lower_bound;
#define ub upper_bound;
#define lbi(A, x) (ll)(lb(all(A), x) - A.begin())
#define ubi(A, x) (ll)(ub(all(A), x) - A.begin())
using ll = long long;
using P = pair<int, int>;
template <class T>
using V = vector<T>;
template <class T>
using VV = V<V<T> >;
template <class T>
const ll INF = 1000000000000000001;
const int inf = 1001001001;
const ll mod = 1000000007;
const ll MOD = 998244353;
inline V<int> dtois(string& s) { V<int> vec = {}; for (auto& e : s) { vec.push_back(e - '0'); } return vec; }

VV<ll> matrix_prod(const VV<ll>& A, const VV<ll>& B) {
	int ac = A.size();	int al = A[0].size();
	int bc = B.size();	int bl = B[0].size();
	VV<ll> res(ac, V<ll>(bl));
	rep(i, 0, ac) {
		rep(j, 0, bl) {
			rep(k, 0, al) {
				res[i][j] += A[i][k] * B[k][j];
				res[i][j] %= MOD;
			}
		}
	}
	return res;
}

VV<ll> matrix_exp(VV<ll>A, ll n) {
	int ac = A.size();
	VV<ll> res(ac, V<ll>(ac));
	rep(i, 0, ac) {
		res[i][i] = 1;
	}
	while (n) {
		if (n & 1) {
			res = matrix_prod(A, res);
		}
		A = matrix_prod(A, A);
		n >>= 1;
	}
	return res;
}

int main(void) {
	ll n, m, t;	cin >> n >> m >> t;
	VV<ll> A(n, V<ll>(n));
	rep(i, 0, m) {
		int a, b;	cin >> a >> b;
		A[a][b] = A[b][a] = 1;
	}
	VV<ll> dp(n, V<ll>(1));	dp[0][0] = 1;
	dp = matrix_prod(matrix_exp(A, t), dp);

	cout << dp[0][0] << endl;
	return 0;
}
0