結果

問題 No.997 Jumping Kangaroo
ユーザー kotatsugame
提出日時 2020-02-21 22:45:47
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 1,677 bytes
コンパイル時間 1,183 ms
コンパイル使用メモリ 73,820 KB
実行使用メモリ 5,248 KB
最終ジャッジ日時 2024-10-08 22:24:47
合計ジャッジ時間 2,989 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 19 WA * 6
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp:11:1: warning: ISO C++ forbids declaration of 'main' with no type [-Wreturn-type]
   11 | main()
      | ^~~~

ソースコード

diff # 

#pragma GCC optimize ("O3")
#pragma GCC target ("avx")
#include<iostream>
using namespace std;
int W,N;
long K;
int A[99];
long mod=1e9+7,D[160][160];
long X[2][160];
long E[160][160],tmp[160][160];
main()
{
	cin>>N>>W>>K;
	for(int i=0;i<N;i++)cin>>A[i];
	for(int k=0;k<2;k++)for(int i=0;i<2*W;i++)
	{
		for(int j=0;j<2;j++)for(int l=0;l<4*W;l++)X[j][l]=0;
		X[k][i]=1;
		for(int j=i;j<2*W;j++)for(int l=0;l<2;l++)
		{
			for(int I=0;I<N;I++)
			{
				if(j+A[I]<4*W)
				{
					int nl=l;
					nl+=(j+A[I])/W-j/W;
					if((j+A[I])%W==0)
					{
						if(nl<=2)nl=0;
						else continue;
					}
					if(nl<=1)(X[nl][j+A[I]]+=X[l][j])%=mod;
				}
			}
		}
		for(int j=0;j<2;j++)for(int l=0;l<2*W;l++)
		{
			D[j*2*W+l][k*2*W+i]=X[j][l+2*W];
		}
	}
	for(int i=0;i<4*W;i++)E[i][i]=1;
	long KK=K/2;
	while(KK>0)
	{
		if(KK&1)
		{
			for(int i=0;i<4*W;i++)for(int j=0;j<4*W;j++)tmp[i][j]=0;
			for(int i=0;i<4*W;i++)for(int j=0;j<4*W;j++)for(int k=0;k<4*W;k++)
			{
				tmp[i][j]+=E[i][k]*D[k][j];
			}
			for(int i=0;i<4*W;i++)for(int j=0;j<4*W;j++)E[i][j]=tmp[i][j]%mod;
		}
		KK>>=1;
		for(int i=0;i<4*W;i++)for(int j=0;j<4*W;j++)tmp[i][j]=0;
		for(int i=0;i<4*W;i++)for(int j=0;j<4*W;j++)for(int k=0;k<4*W;k++)
		{
			tmp[i][j]+=D[i][k]*D[k][j];
		}
		for(int i=0;i<4*W;i++)for(int j=0;j<4*W;j++)D[i][j]=tmp[i][j]%mod;
	}
	long ans[160]={};
	for(int i=0;i<4*W;i++)ans[i]=E[i][0];
	for(int j=0;j<2*W;j++)for(int i=0;i<2;i++)
	{
		for(int I=0;I<N;I++)
		{
			if(j+A[I]<2*W)
			{
				int nl=i;
				nl+=(j+A[I])/W-j/W;
				if((j+A[I])%W==0)
				{
					if(nl<=2)nl=0;
					else continue;
				}
				if(nl<2)(ans[nl*2*W+j+A[I]]+=ans[i*2*W+j])%=mod;
			}
		}
	}
	cout<<ans[K%2*W]<<endl;
}
0