結果

問題 No.42 貯金箱の溜息
ユーザー cielciel
提出日時 2015-11-21 23:33:56
言語 C++11
(gcc 11.4.0)
結果
AC  
実行時間 47 ms / 5,000 ms
コード長 961 bytes
コンパイル時間 194 ms
コンパイル使用メモリ 27,600 KB
実行使用メモリ 4,352 KB
最終ジャッジ日時 2023-10-11 18:11:13
合計ジャッジ時間 969 ms
ジャッジサーバーID
（参考情報） 		judge15 / judge13
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 34 ms
4,352 KB
testcase_01 AC 47 ms
4,348 KB
testcase_02 AC 47 ms
4,352 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <stdio.h>
#define M 3000
#define MOD 1000000009
long long pow_binary_mod(long long x,long long y,long long mod){
	long long z=1;
	for(;y;y>>=1){
		if((y&1)!=0)z=z*x%mod;
		x=x*x%mod;
	}
	return z;
}
int Z[6]={1,5,10,50,100,500};
int A[M+1],REV[13];
int main(){
	for(int i=-6;i<=6;i++)REV[6+i]=pow_binary_mod((i+MOD)%MOD,MOD-2,MOD);
	REV[6]=A[0]=1;
	for(int z=0;z<6;z++)for(int i=0;i<=M-Z[z];i++)A[i+Z[z]]+=A[i];
	int T;
	for(scanf("%d",&T);T--;){
		long long m,x,r,s=0,total=1;
		scanf("%lld",&m);
		if(m<=3000)s=A[m];
		else{
			x=m/500,r=m%500;
			for(int l=0;l<6;++l)total=( (x-l+MOD)%MOD*total%MOD*REV[l+1] )%MOD;
			for(int k=0;k<6;k++){
				long long p=A[k*500+r];
				for(int l=0;l<6;l++)if(l!=k)p=(x-l+MOD)%MOD*p%MOD*pow_binary_mod(k-l+MOD,MOD-2,MOD)%MOD;
				//p=p*total%MOD*pow_binary_mod((x-k+MOD)%MOD,MOD-2,MOD)%MOD;
				s=(s+p)%MOD;
				total=total*REV[k+6+1]%MOD*pow_binary_mod(REV[k+1],MOD-2,MOD)%MOD;
			}
		}
		printf("%lld\n",s);
	}
}
0