結果

問題 No.321 (P,Q)-サンタと街の子供たち
ユーザー tails
提出日時 2021-07-06 14:32:55
言語 C
(gcc 13.3.0)
結果
TLE  
実行時間 -
コード長 1,157 bytes
コンパイル時間 406 ms
コンパイル使用メモリ 32,256 KB
実行使用メモリ 13,752 KB
最終ジャッジ日時 2024-07-01 12:08:30
合計ジャッジ時間 5,043 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 1 TLE * 1 -- * 2
other AC * 1 -- * 40
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.c:26:1: warning: return type defaults to 'int' [-Wimplicit-int]
   26 | main(){
      | ^~~~
main.c: In function 'main':
main.c:75:9: warning: implicit declaration of function 'write' [-Wimplicit-function-declaration]
   75 |         write(1,wp,wbuf+sizeof wbuf-wp);
      |         ^~~~~
main.c:76:9: warning: implicit declaration of function '_exit'; did you mean '_Exit'? [-Wimplicit-function-declaration]
   76 |         _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;{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){
	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