結果

問題 No.42 貯金箱の溜息
コンテスト
ユーザー ciel
提出日時 2015-11-21 23:33:56
言語 C++11
(gcc 15.2.0 + boost 1.89.0)
コンパイル:
g++-15 -O2 -lm -std=gnu++11 -Wuninitialized -DONLINE_JUDGE -o a.out _filename_
実行:
./a.out
結果
AC  
実行時間 37 ms / 5,000 ms
コード長 961 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 390 ms
コンパイル使用メモリ 40,704 KB
実行使用メモリ 6,400 KB
最終ジャッジ日時 2026-04-04 06:35:55
合計ジャッジ時間 668 ms
ジャッジサーバーID
（参考情報） 		judge3_1 / judge1_0
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 3
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#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