結果

問題 No.321 (P,Q)-サンタと街の子供たち
コンテスト
ユーザー 👑 tails
提出日時 2021-07-06 14:34:13
言語 C90(gcc12)
(gcc 12.4.0)
コンパイル:
gcc-12 -O2 -std=c90 -DONLINE_JUDGE -o a.out _filename_ -lm
実行:
./a.out
結果
AC  
実行時間 3 ms / 2,000 ms
コード長 1,191 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 136 ms
コンパイル使用メモリ 34,236 KB
実行使用メモリ 7,844 KB
最終ジャッジ日時 2026-03-08 16:27:58
合計ジャッジ時間 2,008 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 41
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

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

char*mmap();
#define rd_skip() while(*rp++>=48)
#define rd(v) long v=0;{int _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);
}

main(){
	char*rp=mmap(0l,1l<<25,1,2,0,0ll);
	rd(p);
	rd(q);
	rd_skip();
	long r=0;
	if(p+q==0){
		while(*rp){
			long x=*rp==48;
			rd_skip();
			long y=*rp==48;
			rd_skip();
			r+=x&y;
		}
	}else{
		long g=gcd(p,q);
		long b=__builtin_ctzl(g);
		unsigned long t,e;
		{
			long x=g>>b;
			t=~0ul/x;
			e=x;
			e*=2-x*e;
			e*=2-x*e;
			e*=2-x*e;
			e*=2-x*e;
			e*=2-x*e;
		}
		p=(p>>b)*e;
		q=(q>>b)*e;
		long m=(1<<b)-1;
		long f=!(p&q&1);
		while(*rp){
			rp+=*rp==45;
			rd(x);
			rp+=*rp==45;
			rd(y);
			if(!((x|y)&m)){
				x=(x>>b)*e;
				y=(y>>b)*e;
				if(x<=t&&y<=t){
					r+=f|~(x^y)&1;
				}
			}
		}
	}

	char wbuf[64],*wp=wbuf+sizeof wbuf;
	wt(r);
	write(1,wp,wbuf+sizeof wbuf-wp);
	_exit(0);
}
0