結果
問題 | No.891 隣接3項間の漸化式 |
ユーザー | tarattata1 |
提出日時 | 2019-09-20 23:35:17 |
言語 | C++11 (gcc 11.4.0) |
結果 |
AC
|
実行時間 | 2 ms / 2,000 ms |
コード長 | 2,165 bytes |
コンパイル時間 | 824 ms |
コンパイル使用メモリ | 72,784 KB |
実行使用メモリ | 5,376 KB |
最終ジャッジ日時 | 2024-09-14 20:25:42 |
合計ジャッジ時間 | 1,889 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 2 ms
5,248 KB |
testcase_01 | AC | 2 ms
5,376 KB |
testcase_02 | AC | 1 ms
5,376 KB |
testcase_03 | AC | 1 ms
5,376 KB |
testcase_04 | AC | 1 ms
5,376 KB |
testcase_05 | AC | 2 ms
5,376 KB |
testcase_06 | AC | 1 ms
5,376 KB |
testcase_07 | AC | 2 ms
5,376 KB |
testcase_08 | AC | 2 ms
5,376 KB |
testcase_09 | AC | 1 ms
5,376 KB |
testcase_10 | AC | 1 ms
5,376 KB |
testcase_11 | AC | 1 ms
5,376 KB |
testcase_12 | AC | 1 ms
5,376 KB |
testcase_13 | AC | 2 ms
5,376 KB |
testcase_14 | AC | 1 ms
5,376 KB |
testcase_15 | AC | 2 ms
5,376 KB |
testcase_16 | AC | 2 ms
5,376 KB |
testcase_17 | AC | 2 ms
5,376 KB |
testcase_18 | AC | 2 ms
5,376 KB |
testcase_19 | AC | 2 ms
5,376 KB |
testcase_20 | AC | 2 ms
5,376 KB |
testcase_21 | AC | 1 ms
5,376 KB |
testcase_22 | AC | 1 ms
5,376 KB |
testcase_23 | AC | 1 ms
5,376 KB |
testcase_24 | AC | 2 ms
5,376 KB |
testcase_25 | AC | 2 ms
5,376 KB |
testcase_26 | AC | 1 ms
5,376 KB |
testcase_27 | AC | 2 ms
5,376 KB |
testcase_28 | AC | 2 ms
5,376 KB |
testcase_29 | AC | 2 ms
5,376 KB |
testcase_30 | AC | 2 ms
5,376 KB |
testcase_31 | AC | 1 ms
5,376 KB |
testcase_32 | AC | 1 ms
5,376 KB |
testcase_33 | AC | 2 ms
5,376 KB |
testcase_34 | AC | 1 ms
5,376 KB |
testcase_35 | AC | 1 ms
5,376 KB |
testcase_36 | AC | 2 ms
5,376 KB |
testcase_37 | AC | 2 ms
5,376 KB |
testcase_38 | AC | 2 ms
5,376 KB |
testcase_39 | AC | 1 ms
5,376 KB |
testcase_40 | AC | 1 ms
5,376 KB |
testcase_41 | AC | 2 ms
5,376 KB |
コンパイルメッセージ
main.cpp: In function ‘int main(int, char**)’: main.cpp:95:10: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result] 95 | scanf("%d%d%d", &a, &b, &n); | ~~~~~^~~~~~~~~~~~~~~~~~~~~~
ソースコード
#include <stdio.h> #include <string> #include <cstring> #include <stdlib.h> #include <math.h> #include <algorithm> #include <vector> #include <set> #include <map> #include <queue> #include <stack> #include <list> #include <iterator> #include <assert.h> #pragma warning(disable:4996) typedef long long ll; #define MIN(a, b) ((a)>(b)? (b): (a)) #define MAX(a, b) ((a)<(b)? (b): (a)) #define LINF 9223300000000000000 #define INF 2140000000 const long long MOD = 1000000007; using namespace std; ll mpow(ll x, ll n){ //x^n(mod M) ll ans = 1; while(n != 0){ if(n&1) ans = ans*x % MOD; x = x*x % MOD; n = n >> 1; } return ans; } ll minv(ll x){ return mpow( x, MOD-2 ); } void mtxprd( int p, int q, int r, const ll* A, const ll* B, ll* C ) { int i, j, k; for(i=0; i<p; i++) { for(k=0; k<r; k++) { C[i*r+k] = 0; for(j=0; j<q; j++) { C[i*r+k] = (C[i*r+k] + A[i*q+j] * B[j*r+k])%MOD; } } } return; } void mtxuni(int dim, ll* mtx){ int i,j; for(i=0; i<dim; i++) { for(j=0; j<dim; j++) { mtx[i*dim+j]=(i==j? 1: 0); } } return; } void mtxcpy(int dim, const ll* mtx, ll* mtx2){ int i,j; for(i=0; i<dim; i++) { for(j=0; j<dim; j++) { mtx2[i*dim+j]=mtx[i*dim+j]; } } return; } void mtxpow(int dim, const ll* mtx, ll n, ll* mtx2){ //x^n(mod M) vector<ll> mtx0(dim*dim); vector<ll> mtxtmp(dim*dim); mtxcpy( dim, mtx, &mtx0[0] ); mtxuni( dim, mtx2 ); while(n != 0){ if(n&1) { mtxprd( dim, dim, dim, mtx2, &mtx0[0], &mtxtmp[0] ); mtxcpy( dim, &mtxtmp[0], mtx2 ); } mtxprd( dim, dim, dim, &mtx0[0], &mtx0[0], &mtxtmp[0] ); mtxcpy( dim, &mtxtmp[0], &mtx0[0] ); n = n >> 1; } return; } int main(int argc, char* argv[]) { int a,b,n; scanf("%d%d%d", &a, &b, &n); ll A[2*2] = {a, b, 1, 0}; ll B[2*2]; mtxpow( 2, A, n, B); ll x0[2]={1,0}; ll x1[2]; mtxprd( 2, 2, 1, B, x0, x1 ); printf("%lld\n", x1[1]); return 0; }