#572351 (cLay) No.194 フィボナッチ数列の理解(1)

• 問題: No.194 フィボナッチ数列の理解(1)
• ユーザー: tails
• 提出日時: 2020-10-30 00:57:25
• 言語: cLay (20240714-1)
• 結果: AC
• 実行時間: 4 ms / 5,000 ms
• コード長: 538 bytes
• コンパイル時間: 3,719 ms
• コンパイル使用メモリ: 185,476 KB
• 実行使用メモリ: 5,376 KB
• 最終ジャッジ日時: 2024-07-05 14:39:58
• 合計ジャッジ時間: 4,376 ms
• ジャッジサーバーID （参考情報）: judge2 / judge5

入力	結果	実行時間; 実行使用メモリ
testcase_00	AC	1 ms; 5,248 KB
testcase_01	AC	2 ms; 5,376 KB
testcase_02	AC	1 ms; 5,376 KB
testcase_03	AC	2 ms; 5,376 KB
testcase_04	AC	2 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	1 ms; 5,376 KB
testcase_09	AC	1 ms; 5,376 KB
testcase_10	AC	1 ms; 5,376 KB
testcase_11	AC	2 ms; 5,376 KB
testcase_12	AC	1 ms; 5,376 KB
testcase_13	AC	1 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	1 ms; 5,376 KB
testcase_18	AC	1 ms; 5,376 KB
testcase_19	AC	2 ms; 5,376 KB
testcase_20	AC	2 ms; 5,376 KB
testcase_21	AC	4 ms; 5,376 KB
testcase_22	AC	1 ms; 5,376 KB
testcase_23	AC	2 ms; 5,376 KB
testcase_24	AC	3 ms; 5,376 KB
testcase_25	AC	3 ms; 5,376 KB
testcase_26	AC	2 ms; 5,376 KB
testcase_27	AC	3 ms; 5,376 KB
testcase_28	AC	2 ms; 5,376 KB
testcase_29	AC	4 ms; 5,376 KB
testcase_30	AC	1 ms; 5,376 KB
testcase_31	AC	2 ms; 5,376 KB
testcase_32	AC	2 ms; 5,376 KB
testcase_33	AC	1 ms; 5,376 KB
testcase_34	AC	1 ms; 5,376 KB
testcase_35	AC	2 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	2 ms; 5,376 KB

ソースコード

{
	ll@n,@k;
	Modint a[1d4+1];
	rd(a(n));

	if(n>30){
		Modint f,s;
		rep(i,n){
			f+=a[i];
		}
		s=f+f;
		a[n]=f;
		int j=0;
		rep(i,n+1,k){
			a[j]=f=f+f-a[j];
			if(++j>n) j=0;
			s+=f;
		}
		wt(f,s);
	}else{
		Polynomial<Modint> b,c,m;
		b.change(0,1);
		c.change(1,1);
		m.change(0,1);
		m.change(n,-2);
		m.change(n+1,1);
	
		for(int j=40;j>=0;--j){
			b*=b;
			b%=m;
			if(k-n&1l<<j){
				b*=c;
				b%=m;
			}
		}
	
		Modint x,y,z;
		rrep(i,n){
			x+=a[i];
			y-=x*b.c[0];
			b*=c;
			b%=m;
			z-=x*b.c[0];
		}
		wt(z-y,z+x);
	}
}