#717274 (C) No.1731 Product of Subsequence

• 問題: No.1731 Product of Subsequence
• ユーザー: tails
• 提出日時: 2021-11-18 15:48:21
• 言語: C (gcc 12.3.0)
• 結果: AC
• 実行時間: 96 ms / 2,000 ms
• コード長: 1,359 bytes
• コンパイル時間: 464 ms
• コンパイル使用メモリ: 33,024 KB
• 実行使用メモリ: 5,248 KB
• 最終ジャッジ日時: 2024-11-08 05:12:16
• 合計ジャッジ時間: 2,954 ms
• ジャッジサーバーID （参考情報）: judge3 / judge4

入力	結果	実行時間; 実行使用メモリ
testcase_00	AC	96 ms; 5,248 KB
testcase_01	AC	1 ms; 5,248 KB
testcase_02	AC	1 ms; 5,248 KB
testcase_03	AC	2 ms; 5,248 KB
testcase_04	AC	1 ms; 5,248 KB
testcase_05	AC	1 ms; 5,248 KB
testcase_06	AC	1 ms; 5,248 KB
testcase_07	AC	1 ms; 5,248 KB
testcase_08	AC	5 ms; 5,248 KB
testcase_09	AC	1 ms; 5,248 KB
testcase_10	AC	95 ms; 5,248 KB
testcase_11	AC	82 ms; 5,248 KB
testcase_12	AC	2 ms; 5,248 KB
testcase_13	AC	3 ms; 5,248 KB
testcase_14	AC	79 ms; 5,248 KB
testcase_15	AC	1 ms; 5,248 KB
testcase_16	AC	1 ms; 5,248 KB
testcase_17	AC	1 ms; 5,248 KB
testcase_18	AC	89 ms; 5,248 KB
testcase_19	AC	3 ms; 5,248 KB
testcase_20	AC	77 ms; 5,248 KB
testcase_21	AC	95 ms; 5,248 KB
testcase_22	AC	68 ms; 5,248 KB
testcase_23	AC	1 ms; 5,248 KB
testcase_24	AC	3 ms; 5,248 KB
testcase_25	AC	2 ms; 5,248 KB
testcase_26	AC	17 ms; 5,248 KB
testcase_27	AC	24 ms; 5,248 KB
testcase_28	AC	48 ms; 5,248 KB
testcase_29	AC	20 ms; 5,248 KB
testcase_30	AC	91 ms; 5,248 KB
testcase_31	AC	69 ms; 5,248 KB
testcase_32	AC	1 ms; 5,248 KB
testcase_33	AC	5 ms; 5,248 KB
testcase_34	AC	6 ms; 5,248 KB

コンパイルメッセージ

main.c:53:1: warning: return type defaults to 'int' [-Wimplicit-int]
   53 | main(){
      | ^~~~
main.c: In function 'main':
main.c:86:17: warning: implicit declaration of function 'write' [-Wimplicit-function-declaration]
   86 |                 write(1,wp,wbuf+sizeof wbuf-wp);
      |                 ^~~~~
main.c:88:9: warning: implicit declaration of function '_exit'; did you mean '_Exit'? [-Wimplicit-function-declaration]
   88 |         _exit(0);
      |         ^~~~~
      |         _Exit

ソースコード

#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);

}