結果
問題 |
No.1750 ラムドスウイルスの感染拡大-hard
|
ユーザー |
![]() |
提出日時 | 2021-11-19 21:45:00 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 271 ms / 2,000 ms |
コード長 | 2,636 bytes |
コンパイル時間 | 2,109 ms |
コンパイル使用メモリ | 184,252 KB |
実行使用メモリ | 6,820 KB |
最終ジャッジ日時 | 2024-12-31 22:37:37 |
合計ジャッジ時間 | 5,877 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 30 |
ソースコード
#include <bits/stdc++.h> #define fi first #define se second #define pb push_back #define sz(a) (int)a.size() #define all(a) a.begin(),a.end() #define rep(i,n) for(int i=0;i<n;i++) #define crep(i,x,n) for(int i=x;i<n;i++) #define drep(i,n) for(int i=n-1;i>=0;i--) #define vec(...) vector<__VA_ARGS__> #define _3HaBFkZ ios_base::sync_with_stdio(false),cin.tie(0),cout.tie(0) using namespace std; typedef long long ll; typedef long double ld; using pii=pair<int,int>; using vi=vector<int>; //snuke's modular int template <ll mod> struct modularint{ ll x; modularint(ll x=0):x(x%mod){} modularint& operator+=(const modularint a){ if ((x += a.x) >= mod) x -= mod; return *this; } modularint& operator-=(const modularint a){ if ((x += mod-a.x) >= mod) x -= mod; return *this; } modularint& operator*=(const modularint a){ (x *= a.x) %= mod; return *this; } modularint operator+(const modularint a)const{ modularint res(*this); return res+=a; } modularint operator-(const modularint a)const{ modularint res(*this); return res-=a; } modularint operator*(const modularint a)const{ modularint res(*this); return res*=a; } modularint pow(ll n)const{ modularint res=1,x(*this); while(n){ if(n&1)res*=x; x*=x; n>>=1; } return res; } modularint inv()const{ return pow(mod-2); } }; using mint=modularint<998244353>; template <typename T> struct matrix{ vec(vec(T)) a; //initilize matrix here matrix(){} matrix(int h,int w){ a.clear(); a.resize(h,vec(T)(w)); } matrix(vec(vec(T)) nea){ a=nea; } //outer vector size here int size()const{ return a.size(); } const vector<T>& operator[](int i)const{ return a[i]; } vector<T>& operator[](int i){ return a[i]; } matrix<T>& operator *=(const matrix<T>& rhs){ int h=a.size(), w=rhs[0].size(), c=rhs.size(); matrix<T> res(h,w); rep(i,h){ rep(j,w){ rep(k,c){ res[i][j] += a[i][k] * rhs[k][j]; } } } this->a = res.a; return *this; } matrix<T> operator *(const matrix<T>& rhs){ return (matrix<T>(*this) *= rhs); } }; template <typename T> matrix<T> idenmat(int n){ matrix<T> res(n,n); rep(i,n) res.a[i][i] = 1; return res; }; template <typename T> matrix<T> pow(matrix<T> mat,ll n){ matrix<T> res=idenmat<T>(mat.size()); while(n){ if(n&1) res*=mat; mat*=mat; n>>=1; } return res; }; int main(){ _3HaBFkZ; ll n,m,t; cin>>n>>m>>t; vec(vec(mint)) wys(n,vec(mint)(n,0)); rep(i,m){ int u,v; cin>>u>>v; wys[u][v]=wys[u][v]+1; wys[v][u]=wys[v][u]+1; } matrix<mint> mat(wys); matrix<mint> now=idenmat<mint>(n); now=now*pow(mat,t); cout<<now[0][0].x<<"\n"; // return 0; }