結果

問題 No.2788 4-33 Hard
コンテスト
ユーザー tails
提出日時 2024-06-17 19:52:19
言語 C90
(gcc 12.3.0)
結果
AC  
実行時間 1 ms / 2,000 ms
コード長 1,240 bytes
コンパイル時間 364 ms
コンパイル使用メモリ 27,264 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-06-17 19:52:22
合計ジャッジ時間 2,364 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 52
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.c: In function ‘main’:
main.c:8:59: warning: implicit declaration of function ‘write’ [-Wimplicit-function-declaration]
    8 | #define wt1(v) ({char wbuf[64],*wp=wbuf+sizeof wbuf;wt(v);write(1,wp,wbuf+sizeof wbuf-wp);})
      |                                                           ^~~~~
main.c:50:9: note: in expansion of macro ‘wt1’
   50 |         wt1(z);
      |         ^~~
main.c:51:9: warning: implicit declaration of function ‘_exit’ [-Wimplicit-function-declaration]
   51 |         _exit(0);
      |         ^~~~~
main.c:51:9: warning: incompatible implicit declaration of built-in function ‘_exit’ [-Wbuiltin-declaration-mismatch]

ソースコード

diff # 

#pragma GCC optimize("Ofast")
#pragma GCC target("avx2")

#define rd_init() char*rp=({char*mmap();mmap(0l,1l<<25,1,2,0,0ll);})
#define rd_skip() while(*rp++>=48)
#define rd() ({long _v=0,_c;while(_c=*rp++-48,_c>=0)_v=_v*10+_c;_v;})
#define wt(v) ({ulong _z=v;do*--wp=_z%10+48;while(_z/=10);})
#define wt1(v) ({char wbuf[64],*wp=wbuf+sizeof wbuf;wt(v);write(1,wp,wbuf+sizeof wbuf-wp);})
#define rep(v,e) for(typeof(e)v=0;v<e;++v)
#define reps(v,s,e) for(typeof(e)v=s;v<e;++v)
#define rrep(v,e) for(typeof(e)v=e;--v>=0;)
#define rreps(v,s,e) for(typeof(e)v=e;--v>=s;)
#define MD 998244353

typedef unsigned long ulong;

unsigned const v[10]={0,1,499122177,332748118,748683265,598946612,166374059,855638017,873463809,443664157};
unsigned d[9][5][34];

int main(){
	d[0][0][0]=1;
	rd_init();
	rep(a,5){
		rep(b,34){
			ulong o=(ulong)rd()%MD;
			ulong c[9];
			c[0]=1;
			rep(i,8){
				c[i+1]=c[i]*((o-i+MD)%MD)%MD*v[i+1]%MD;
			}
			rreps(i,1,9){
				reps(j,1,i+1){
					rrep(x,5-a*j){
						rrep(y,34-b*j){
							d[i][x+a*j][y+b*j]=(d[i][x+a*j][y+b*j]+d[i-j][x][y]*c[j])%MD;
						}
					}
				}
			}
		}
	}
	ulong z=0;
	rep(a,5){
		ulong o=(ulong)rd()%MD;
		z=(z+d[8][4-a][33]*o)%MD;
		rep(b,33){
			rd_skip();
		}
	}
	wt1(z);
	_exit(0);
}
0