結果
問題 | 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 inttemplate <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 herematrix(){}matrix(int h,int w){a.clear();a.resize(h,vec(T)(w));}matrix(vec(vec(T)) nea){a=nea;}//outer vector size hereint 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;}