結果
| 問題 | 
                            No.321 (P,Q)-サンタと街の子供たち
                             | 
                    
| コンテスト | |
| ユーザー | 
                             tails
                         | 
                    
| 提出日時 | 2021-07-06 14:34:13 | 
| 言語 | C  (gcc 13.3.0)  | 
                    
| 結果 | 
                             
                                AC
                                 
                             
                            
                         | 
                    
| 実行時間 | 4 ms / 2,000 ms | 
| コード長 | 1,191 bytes | 
| コンパイル時間 | 199 ms | 
| コンパイル使用メモリ | 32,180 KB | 
| 実行使用メモリ | 6,944 KB | 
| 最終ジャッジ日時 | 2024-07-01 12:08:32 | 
| 合計ジャッジ時間 | 1,663 ms | 
| 
                            ジャッジサーバーID (参考情報)  | 
                        judge4 / judge5 | 
(要ログイン)
| ファイルパターン | 結果 | 
|---|---|
| sample | AC * 4 | 
| other | AC * 41 | 
コンパイルメッセージ
main.c:29:1: warning: return type defaults to 'int' [-Wimplicit-int]
   29 | main(){
      | ^~~~
main.c: In function 'main':
main.c:78:9: warning: implicit declaration of function 'write' [-Wimplicit-function-declaration]
   78 |         write(1,wp,wbuf+sizeof wbuf-wp);
      |         ^~~~~
main.c:79:9: warning: implicit declaration of function '_exit'; did you mean '_Exit'? [-Wimplicit-function-declaration]
   79 |         _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;{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);
}
            
            
            
        
            
tails