1. 程式人生 > >CodeForces - 1025B Weakened Common Divisor

CodeForces - 1025B Weakened Common Divisor

對數 its ++ pre 可能 至少 scanf ref max

http://codeforces.com/problemset/problem/1025/B

大意:n對數對(ai,bi),求任意一個數滿足是所有數對中至少一個數的因子(大於1)

分析:

  1. 首先求所有數對的lcm,把所有數的素因子提出來
  2. 求所有lcm的gcd,這樣做求出數對之間的公共素因子gcd
  3. 註意,公共素因子可能在某一組數對中狀態為某一部分是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