結果
問題 | No.1073 無限すごろく |
ユーザー |
![]() |
提出日時 | 2020-06-05 23:07:25 |
言語 | C++11(廃止可能性あり) (gcc 13.3.0) |
結果 |
AC
|
実行時間 | 2 ms / 2,000 ms |
コード長 | 2,387 bytes |
コンパイル時間 | 1,339 ms |
コンパイル使用メモリ | 160,220 KB |
実行使用メモリ | 6,820 KB |
最終ジャッジ日時 | 2024-12-17 17:40:58 |
合計ジャッジ時間 | 2,388 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 30 |
ソースコード
#include "bits/stdc++.h" using namespace std; #define int long long #define FOR(i, a, b) for(int i=(a);i<(b);i++) #define RFOR(i, a, b) for(int i=(b-1);i>=(a);i--) #define REP(i, n) for(int i=0; i<(n); i++) #define RREP(i, n) for(int i=(n-1); i>=0; i--) #define REP1(i, n) for(int i=1; i<=(n); i++) #define RREP1(i, n) for(int i=(n); i>=1; i--) #define ALL(a) (a).begin(),(a).end() #define UNIQUE_SORT(l) sort(ALL(l)); l.erase(unique(ALL(l)), l.end()); #define CONTAIN(a, b) find(ALL(a), (b)) != (a).end() #define out(...) printf(__VA_ARGS__) #define chmax(a,b) a = max(a,b) #define chmin(a,b) a = min(a,b) #if DEBUG #define debug(...) printf(__VA_ARGS__) #else #define debug(...) /* ... */ #endif void solve(); signed main() { #if DEBUG std::ifstream in("input.txt"); std::cin.rdbuf(in.rdbuf()); #endif cin.tie(0); ios::sync_with_stdio(false); solve(); return 0; } /*================================*/ #if DEBUG #define SIZE 6 #else #define SIZE 6 #endif int N; const int MAX = 10; const int MOD = 1e9+7; long long fac[MAX], finv[MAX], inv[MAX]; void COMinit() { fac[0] = fac[1] = 1; finv[0] = finv[1] = 1; inv[1] = 1; for (int i = 2; i < MAX; i++){ fac[i] = fac[i - 1] * i % MOD; inv[i] = MOD - inv[MOD%i] * (MOD / i) % MOD; finv[i] = finv[i - 1] * inv[i] % MOD; } } long long COM(int n, int k){ if (n < k) return 0; if (n < 0 || k < 0) return 0; if (!fac[n]) COMinit(); return fac[n] * (finv[k] * finv[n - k] % MOD) % MOD; } int P(int n, int k) { if (n < k) return 0; if (n < 0 || k < 0) return 0; if (!fac[n]) COMinit(); return (fac[n] * finv[n-k]) % MOD; } void matmul(int a[SIZE][SIZE], int b[SIZE][SIZE]) { int tmp[SIZE][SIZE]; REP(x,SIZE)REP(y,SIZE) tmp[x][y] = 0; REP(x,SIZE)REP(y,SIZE)REP(z,SIZE) (tmp[x][y] += a[x][z]*b[z][y]%MOD) %= MOD; REP(x,SIZE)REP(y,SIZE) a[x][y]=tmp[x][y]; } void matmul(int a[SIZE][SIZE], int b[SIZE]) { int tmp[SIZE]; REP(i,SIZE) tmp[i]=0; REP(i,SIZE)REP(j,SIZE) tmp[i] += (a[i][j]*b[j]) % MOD; REP(i,SIZE) b[i]=tmp[i]%MOD; } int a[SIZE][SIZE]; int ans[SIZE]={0,0,0,0,0,1}; void solve() { cin>>N; COMinit(); REP(i,5) a[i][i+1]=1; REP(i,SIZE) { a[5][i] = inv[6]; } REP(i,61) { if (N&(1LL<<i)) matmul(a, ans); matmul(a, a); } cout << ans[5] << endl; }