結果

問題 No.1200 お菓子配り-3
コンテスト
ユーザー 👑 tails
提出日時 2021-05-19 19:10:23
言語 C
(gcc 13.3.0)
結果
AC  
実行時間 78 ms / 4,000 ms
コード長 1,695 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 424 ms
コンパイル使用メモリ 33,536 KB
実行使用メモリ 5,248 KB
最終ジャッジ日時 2024-10-09 23:46:39
合計ジャッジ時間 1,792 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 31
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.c:74:1: warning: return type defaults to 'int' [-Wimplicit-int]
   74 | main(){
      | ^~~~
main.c: In function 'main':
main.c:106:9: warning: implicit declaration of function 'write' [-Wimplicit-function-declaration]
  106 |         write(1,wbuf,wp-wbuf);
      |         ^~~~~
main.c:107:9: warning: implicit declaration of function '_exit'; did you mean '_Exit'? [-Wimplicit-function-declaration]
  107 |         _exit(0);
      |         ^~~~~
      |         _Exit

ソースコード

diff #
raw source code 

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

char*mmap();
#define rd_skip() while(*rp++>=48)
#define rd(v) int v=0;{int _c;while(_c=*rp++-48,_c>=0)v=v*10+_c;}

char wbuf[1<<28];

#define WTHI(v) {int _z=v,_n=0;long _d=0;while(++_n,_d=_d<<8|0x30|_z%10,_z/=10);*(long*)wp=_d;wp+=_n;}
#define WTLO(v) {int _z=v,_n=8;long _d=0;while(_d=_d<<8|0x30|_z%10,_z/=10,--_n);*(long*)wp=_d;wp+=8;}
#define wt(v) WTHI(v)

#define rep(v,e) for(int v=0;v<e;++v)

#define SIEVE_N 46340
int nprimes;
int primes[SIEVE_N];
char sieve[SIEVE_N];

void mkprimes(){
	primes[nprimes++]=2;
	for(int i=3;i<SIEVE_N;i+=2){
		if(sieve[i]==0){
			primes[nprimes++]=i;
			for(int j=i*i;j<SIEVE_N;j+=i*2){
				sieve[j]=1;
			}
		}
	}
}

int divisors(int a,int*divs){
	int n=0;
	divs[n++]=1;
	int i=0;
	int p;
	while(p=primes[i],a>=p*p){
		if(a%p==0){
			int k=1;
			while(a/=p,a%p==0){
				++k;
			}
			int n0=n*k;
			for(int j=0;j<n0;++j){
				divs[n++]=divs[j]*p;
			}
		}
		++i;
	}
	if(a>1){
		int n0=n;
#if 0
		for(int j=0;j<n0;++j){
			divs[n++]=divs[j]*a;
		}
#else
		int j=0;
		while(j<n0-1){
			*(long*)&divs[n]=*(long*)&divs[j]*a;
			n+=2;
			j+=2;
		}
		if(j<n0){
			divs[n++]=divs[j]*a;
		}
#endif
	}
	return n;
}

int d[100000];

main(){
	mkprimes();

	char*rp=mmap(0l,1l<<28,1,2,0,0ll);
	char*wp=wbuf;

	rd_skip();
	while(*rp){
		int res=0;
		rd(x);
		rd(y);
		if(x==y){
			int nd=divisors(x,d);
			res+=x+(x&1)+nd-3;
		}else{
			if(x>y){
				int t=x;x=y;y=t;
			}
			int nd=divisors(y-x,d);
			rep(i,nd){
				int a=d[i]+1;
				int c=(y-x)/(a-1);
				if(x>c){
					int b=(x-c)/(a+1);
					c+=b;
					res+=a*b+c==x&&a*c+b==y;
				}
			}
		}
		wt(res);
		*wp++=10;
	}
	write(1,wbuf,wp-wbuf);
	_exit(0);
}
0