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

• 問題: No.194 フィボナッチ数列の理解(1)
• ユーザー: tails
• 提出日時: 2020-10-30 01:19:18
• 言語: cLay (20240714-1)
• 結果: AC
• 実行時間: 5 ms / 5,000 ms
• コード長: 774 bytes
• コンパイル時間: 3,870 ms
• コンパイル使用メモリ: 187,108 KB
• 実行使用メモリ: 5,376 KB
• 最終ジャッジ日時: 2024-07-05 14:40:59
• 合計ジャッジ時間: 4,926 ms
• ジャッジサーバーID （参考情報）: judge3 / judge1

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

ソースコード

#include <sys/mman.h>

Modint a[1d4+1];
{
	char*rp=static_cast<char*>(mmap(0,24l<<10,1,2,0,0));
	ll n=0; {int c;while(c=*rp++-48,c>=0)n=n*10+c;}
	ll k=0; {int c;while(c=*rp++-48,c>=0)k=k*10+c;}

	if(n>30){
		Modint f,s;
		rep(i,n){
			int v=*rp-48;
			rp+=2;
			f+=a[i]=v;
		}
		a[n]=f;
		s=f+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{
		rep(i,n){
			int v=*rp-48;
			rp+=2;
			a[i]=v;
		}

		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);
	}
}