CodeForces - 1025B Weakened Common Divisor
阿新 • • 發佈:2019-03-03
對數 its ++ pre 可能 至少 scanf ref max
http://codeforces.com/problemset/problem/1025/B
大意:n對數對(ai,bi),求任意一個數滿足是所有數對中至少一個數的因子(大於1)
分析:
- 首先求所有數對的lcm,把所有數的素因子提出來
- 求所有lcm的gcd,這樣做求出數對之間的公共素因子gcd
- 註意,公共素因子可能在某一組數對中狀態為某一部分是ai的素因子而剩下的一部分是bi的素因子,這種情況導致gcd既不是ai的因子又不是bi的因子,因此必須只保留下aibi公共的素因子,剔除非公共素因子,故只需執行gcd=__gcd(ai,gcd),gcd=__gcd(bi,gcd)(註意,需先判斷__gcd(xi,gcd)>1)
#include <bits/stdc++.h> using namespace std; typedef long long ll; const int MAXN=150000+100; ll a[MAXN],b[MAXN]; int main() { int n; scanf("%d",&n); for(int i=0;i<n;i++) scanf("%lld %lld",&a[i],&b[i]); ll gcd=a[0]*b[0]/__gcd(a[0],b[0]); for(int i=1;i<n;i++) gcd=__gcd(gcd,a[i]*b[i]/__gcd(a[i],b[i])); if(gcd==1) { printf("-1\n"); return 0; } ll temp; for(int i=0;i<n;i++) { temp=__gcd(gcd,a[i]); if(temp>1){ gcd=temp; continue; } temp=__gcd(gcd,b[i]); if(temp>1) gcd=temp; } printf("%lld\n",gcd); return 0; }
CodeForces - 1025B Weakened Common Divisor