結果

問題 No.1310 量子アニーリング
コンテスト
ユーザー tails
提出日時 2020-12-08 10:53:55
言語 C
(gcc 13.3.0)
結果
AC  
実行時間 5 ms / 2,000 ms
コード長 745 bytes
コンパイル時間 157 ms
コンパイル使用メモリ 30,592 KB
実行使用メモリ 6,820 KB
最終ジャッジ日時 2024-10-10 08:38:33
合計ジャッジ時間 1,183 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 21
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.c:37:1: warning: return type defaults to 'int' [-Wimplicit-int]
   37 | main(){
      | ^~~~
main.c: In function 'main':
main.c:40:9: warning: implicit declaration of function 'scanf' [-Wimplicit-function-declaration]
   40 |         scanf("%d",&n);
      |         ^~~~~
main.c:1:1: note: include '<stdio.h>' or provide a declaration of 'scanf'
  +++ |+#include <stdio.h>
    1 | #pragma GCC optimize("Ofast")
main.c:40:9: warning: incompatible implicit declaration of built-in function 'scanf' [-Wbuiltin-declaration-mismatch]
   40 |         scanf("%d",&n);
      |         ^~~~~
main.c:40:9: note: include '<stdio.h>' or provide a declaration of 'scanf'
main.c:51:9: warning: implicit declaration of function 'printf' [-Wimplicit-function-declaration]
   51 |         printf("%d",(int)((y+z)%MD*fac[n]*2%MD));
      |         ^~~~~~
main.c:51:9: note: include '<stdio.h>' or provide a declaration of 'printf'
main.c:51:9: warning: incompatible implicit declaration of built-in function 'printf' [-Wbuiltin-declaration-mismatch]
main.c:51:9: note: include '<stdio.h>' or provide a declaration of 'printf'

ソースコード

diff # 

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

#define MD 998244353

int inverse(int a){
	int b=MD;
	int u=1;
	int v=0;
	int s,t;
	while(b){
		t=a/b;
		s=b; b=a-t*b; a=s;
		s=v; v=u-t*v; u=s;
	}
	if(u<0){
		u+=MD;
	}
	return u;
}

#define M 200001

int fac[M],ifac[M];

void mkfac(){
	fac[0]=1;
	for(int i=1;i<M;i++){
		fac[i]=(long)fac[i-1]*i%MD;
	}
	ifac[M-1]=inverse(fac[M-1]);
	for(int i=M;--i;){
		ifac[i-1]=(long)ifac[i]*i%MD;
	}
}

main(){
	mkfac();
	int n,h;
	scanf("%d",&n);
	h=(n>>1)+1;
	long y=0,z=0;
	for(int i=0;i<h;i+=2){
		y=(y*16+(long)ifac[i]*(long)ifac[n-i])%MD;
	}
	y*=1<<(n&3);
	for(int i=n&~1;i>=h;i-=2){
		z=(z*16+(long)ifac[i]*(long)ifac[n-i])%MD;
	}
	z*=16>>(n&3);
	printf("%d",(int)((y+z)%MD*fac[n]*2%MD));
}
0