結果

問題 No.1731 Product of Subsequence
ユーザー tailstails
提出日時 2021-11-18 15:48:21
言語 C
(gcc 12.3.0)
結果
AC  
実行時間 97 ms / 2,000 ms
コード長 1,359 bytes
コンパイル時間 913 ms
コンパイル使用メモリ 31,804 KB
実行使用メモリ 4,384 KB
最終ジャッジ日時 2023-08-07 23:40:13
合計ジャッジ時間 3,010 ms
ジャッジサーバーID
（参考情報） 		judge15 / judge11
このコードへのチャレンジ（β）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 97 ms
4,376 KB
testcase_01 AC 1 ms
4,380 KB
testcase_02 AC 0 ms
4,376 KB
testcase_03 AC 1 ms
4,380 KB
testcase_04 AC 1 ms
4,380 KB
testcase_05 AC 0 ms
4,376 KB
testcase_06 AC 1 ms
4,376 KB
testcase_07 AC 1 ms
4,380 KB
testcase_08 AC 5 ms
4,380 KB
testcase_09 AC 1 ms
4,376 KB
testcase_10 AC 96 ms
4,380 KB
testcase_11 AC 82 ms
4,380 KB
testcase_12 AC 1 ms
4,376 KB
testcase_13 AC 2 ms
4,376 KB
testcase_14 AC 80 ms
4,376 KB
testcase_15 AC 0 ms
4,380 KB
testcase_16 AC 0 ms
4,376 KB
testcase_17 AC 0 ms
4,380 KB
testcase_18 AC 94 ms
4,380 KB
testcase_19 AC 2 ms
4,380 KB
testcase_20 AC 80 ms
4,376 KB
testcase_21 AC 97 ms
4,384 KB
testcase_22 AC 74 ms
4,380 KB
testcase_23 AC 1 ms
4,376 KB
testcase_24 AC 3 ms
4,380 KB
testcase_25 AC 2 ms
4,376 KB
testcase_26 AC 17 ms
4,384 KB
testcase_27 AC 23 ms
4,380 KB
testcase_28 AC 48 ms
4,376 KB
testcase_29 AC 20 ms
4,380 KB
testcase_30 AC 92 ms
4,380 KB
testcase_31 AC 74 ms
4,380 KB
testcase_32 AC 1 ms
4,380 KB
testcase_33 AC 5 ms
4,380 KB
testcase_34 AC 6 ms
4,376 KB
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.c:53:1: 警告: 戻り値の型をデフォルトの ‘int’ にします [-Wimplicit-int]
   53 | main(){
      | ^~~~
main.c: 関数 ‘main’ 内:
main.c:86:17: 警告: 関数 ‘write’ の暗黙的な宣言です [-Wimplicit-function-declaration]
   86 |                 write(1,wp,wbuf+sizeof wbuf-wp);
      |                 ^~~~~
main.c:88:9: 警告: implicit declaration of function ‘_exit’; did you mean ‘_Exit’? [-Wimplicit-function-declaration]
   88 |         _exit(0);
      |         ^~~~~
      |         _Exit

ソースコード

diff # 

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

char*mmap();

#define rd_skip() while(*rp++>=48)
#define rd(v) long v=0;{long _c;while(_c=*rp++-48,_c>=0)v=v*10+_c;}
#define wt(v) {long _z=v;do*--wp=_z%10+48;while(_z/=10);}

long gcd(long a,long b){
	if(a==0||b==0){
		return a+b;
	}
	long x=__builtin_ctzl(a);
	long y=__builtin_ctzl(b);
	a>>=x;
	b>>=y;
	long c;
	while(c=a-b){
		c>>=__builtin_ctzl(c);
		if(c<0){
			b=-c;
		}else{
			a=c;
		}
	}
	return a<<(x<=y?x:y);
}

#define H_BITS 13
#define H_LEN (1<<H_BITS)
#define H_MASK (H_LEN-1)
int h[H_LEN];

long geth(int v){
//	long l=v&H_MASK;
	long l=v+(v>>H_BITS)&H_MASK;
	while(h[l]!=0&&h[l]!=v){
		l=l+1&H_MASK;
	}
	h[l]=v;
	return l;
}

int divn;
int divs[2000];
int divh[2000];

int dp[H_LEN];

#define MD 1000000007

main(){
	char*rp=mmap(0l,1l<<25,1,2,0,0ll);
	rd_skip();
	rd(k);
	for(long i=1;i*i<=k;++i){
		if(k%i==0){
			divs[divn++]=i;
		}
	}
	for(long j=divn-(divs[divn-1]*divs[divn-1]==k);--j>=0;){
		divs[divn++]=k/divs[j];
	}

	dp[geth(1)]=1;

	while(*rp){
		rd(a);
		a=gcd(k,a);
		for(long j=divn;--j>=0;){
			long c=divs[j];
			long ch=geth(c);
			long b=gcd(k,a*c);
			long bh=geth(b);
			long v=dp[bh]+dp[ch];
			v-=v<MD?0:MD;
			dp[bh]=v;
		}
	}

	{
		long z=dp[geth(k)]-(k==1);
		char wbuf[64],*wp=wbuf+sizeof wbuf;
		wt(z);
		write(1,wp,wbuf+sizeof wbuf-wp);
	}
	_exit(0);

}
0