結果

問題 No.2183 LCA on Rational Tree
ユーザー tails
提出日時 2023-01-17 20:13:41
言語 cLay
(20241019-1)
結果
AC  
実行時間 64 ms / 2,000 ms
コード長 801 bytes
コンパイル時間 3,236 ms
コンパイル使用メモリ 196,152 KB
実行使用メモリ 6,820 KB
最終ジャッジ日時 2025-01-02 22:39:18
合計ジャッジ時間 4,008 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 6
権限があれば一括ダウンロードができます

ソースコード

diff # 

unsigned primes[3500];
ll nprimes;

ll factors(unsigned v,unsigned*a){
	ll n=0;
	rep[primes](p,nprimes){
		if(v%p==0){
			a[n++]=p;
			do v/=p;while(v%p==0);
		}
	}
	if(v!=1){
		a[n++]=v;
	}
	return n;
}

struct S{
	unsigned p,q,l,n,f[32],e[32];
	S(){
		rd(p,q);
		//n=Factor(q-=p,f);
		n=factors(q-=p,f);
		l=0;
	}
	void step(){
		p+=l;
		rep[f](x,n){
			while(p%x==0&&q%x==0){
				p/=x;
				q/=x;
			}
		}
		l=1<<30;
		rep[f](x,n){
			if(q%x==0){
				l<?=x-p%x;
			}
		}
	}
};

bool h(S&u,S&v){
	if(u.p<=v.p<u.p+u.l){
		wt(v.p,v.p+v.q);
		return true;
	}
	return false;
}

{
	nprimes=Prime(31622,(int*)primes);
	ll@q;
	rep(q){
		S u,v;
		while(1){
			int d=v.q-u.q;
			if(d==0){
				if(h(u,v)||h(v,u)){
					break;
				}
			}
			if(d<=0){
				u.step();
			}
			if(d>=0){
				v.step();
			}
		}
	}
}
0