結果

問題 No.1731 Product of Subsequence
ユーザー tails
提出日時 2021-11-18 15:48:21
言語 C
(gcc 13.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
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 31
権限があれば一括ダウンロードができます
コンパイルメッセージ
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

ソースコード

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