UR #3 核聚变反应强度( gcd )

tags: -分解质因数 , gcd

题目大意

给定\(n\)个数，求\(a_1\)与\(a_i\)次小公约数

分析

易知次小公约数是\(\gcd\)的因数，于是用\(\gcd\)除去它的最小质因子即可。

没有次小公约数的情况是\(\gcd = 1\)，特判一下即可

直接枚举的时间复杂度为\(O(n \sqrt a)\)

由于数据规模较大考虑优化

由于是求\(sgcd(a_1,a_i)\)于是结果一定是\(a_1\)的质因数组成，于是预处理\(a_1\)的质因数，然后每次处理时除去最小的即可，\(10^{12}< 2^{38}\)于是可以知道得到质因数的个数小于\(38\)个，于是时间复杂度就变为了\(O(50n)\)啦！

代码

#include<iostream>
#include<cmath>
#include<cstring>
#include<cstdio>
#include<algorithm>
#define LL long long
#define maxn 100010
#define rep(i,a,b) for(int i = a; i <= b ;i++)

using namespace std;

int n;
LL ys[40];

LL gcd(LL a, LL b) {
    while (b ^= a ^= b ^= a %= b);
    return a;
}

int main(){
    scanf("%d",&n);

    LL l,tmp;   int tot = 0;
    scanf("%lld", &l);
    tmp = l;
    for (int i = 2; 1ll * i * i <= l ;i ++)
    if (l % i == 0)
    {
        ys[++tot] = i;
        while (l % i == 0) l /= i;
    }
    if (l != 1) ys[++tot] = l;
    l = tmp ;
    if (l != 1) printf("%lld ",l / ys[1]);
      else printf("-1 ");

    rep(i,2,n){
      LL aa;
      scanf("%lld",&aa);

      LL g = gcd(aa,l);
      if (g == 1) printf("-1 ");
      else rep (j,1,tot)
             if (g % ys[j] == 0 )
              {
                printf("%lld ",g / ys[j]);
                break;
              }
    }

    return 0;
}