結果
| 問題 |
No.1879 How many matchings?
|
| コンテスト | |
| ユーザー |
serufe_ni
|
| 提出日時 | 2022-03-22 20:26:48 |
| 言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
WA
|
| 実行時間 | - |
| コード長 | 1,142 bytes |
| コンパイル時間 | 1,521 ms |
| コンパイル使用メモリ | 167,780 KB |
| 実行使用メモリ | 6,824 KB |
| 最終ジャッジ日時 | 2024-10-10 14:40:13 |
| 合計ジャッジ時間 | 2,163 ms |
|
ジャッジサーバーID (参考情報) |
judge4 / judge3 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| other | AC * 8 WA * 7 |
ソースコード
#include <bits/stdc++.h>
using namespace std;
#define int long long
const int mod = 1000000007;
int n;
int martxB[2][2] = {1, 1, 1, 0};
int martxA[2] = {1, 1};
int ans[2];
void multiply(int x[][2], int y[][2]) {
int res[2][2];
memset(res, 0, sizeof res);
for (int i = 0; i < 2; i ++) {
for (int j = 0; j < 2; j ++) {
for (int k = 0; k < 2; k ++) {
res[i][j] = (x[i][k] * y[k][j] + res[i][j]) % mod;
}
}
}
memcpy(x,res,sizeof res);
}
void Qpow(int x[][2], int k) {
int res[2][2];
memset(res, 0, sizeof res);
for (int i = 0; i < 2; i ++) {
res[i][i] = 1;
}
while (k) {
if (k % 2) {
multiply(res, x);
}
multiply(x, x);
k /= 2;
}
memcpy(x, res, sizeof res);
}
signed main () {
cin >> n;
if (n <= 2) {
cout << 1;
return 0;
}
Qpow(martxB, n / 2 - 1);
for (int i = 0; i < 2; i ++) {
for (int k =0; k < 2; k ++) {
ans[i] = (ans[i] + martxB[i][k] * martxA[k] % mod ) % mod;
}
}
cout << ans[0] * (n % 2 ? n:1) << endl;
}
serufe_ni