結果
問題 | No.579 3 x N グリッド上のサイクルのサイズ(hard) |
ユーザー |
![]() |
提出日時 | 2020-02-17 08:54:21 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 3 ms / 2,000 ms |
コード長 | 1,598 bytes |
コンパイル時間 | 1,297 ms |
コンパイル使用メモリ | 117,268 KB |
実行使用メモリ | 5,248 KB |
最終ジャッジ日時 | 2024-10-06 14:48:54 |
合計ジャッジ時間 | 3,471 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 80 |
ソースコード
#include <cstdio> #include <cstring> #include <iostream> #include <string> #include <cmath> #include <bitset> #include <vector> #include <map> #include <set> #include <queue> #include <deque> #include <algorithm> #include <complex> #include <unordered_map> #include <unordered_set> #include <random> #include <cassert> #include <fstream> #include <utility> #include <functional> #include <time.h> #include <stack> #include <array> #define popcount __builtin_popcount using namespace std; typedef long long int ll; typedef pair<int, int> P; const ll MOD=1e9+7; vector<vector<ll>> matrixmul(int l, int m, int n, vector<vector<ll>> a, vector<vector<ll>> b){ vector<vector<ll>> c(l, vector<ll>(n)); for(int i=0; i<l; i++){ for(int j=0; j<n; j++){ for(int k=0; k<m; k++){ (c[i][j]+=a[i][k]*b[k][j])%=MOD; } } } return c; } vector<vector<ll>> matrixpow(int n, vector<vector<ll>> a, ll k){ vector<vector<ll>> ap=a, ans(n, vector<ll>(n)); for(int i=0; i<n; i++) ans[i][i]=1; while(k){ if(k&1) ans=matrixmul(n, n, n, ap, ans); ap=matrixmul(n, n, n, ap, ap); k>>=1; } return ans; } ll a[]={32,316,2292,14422,84744,479004,2638328,14258574,75940592,399782668,84795558,786749020}; ll c[]={4,4,999999968,56,62,999999621,753,999999123,685,999999665,102,999999991}; int main() { ll n; cin>>n; vector<vector<ll>> mat(12, vector<ll>(12)); for(int i=0; i<11; i++) mat[i][i+1]=1; for(int i=0; i<12; i++) mat[11][i]=MOD-c[i]; auto matp=matrixpow(12, mat, n-1); ll ans=0; for(int i=0; i<12; i++) (ans+=matp[0][i]*a[i])%=MOD; cout<<ans<<endl; return 0; }