UOJ 48 次最大公约数

次最大公约数 = gcd / 其中一个数质因数中最小的。

gcd(42,12) = 6;    div(42) = 2*3*7   div(12) = 2^2*3

sgcd(42,12) = 6 / 2 = 3;

之前素数筛选，分解质因数总是找模板，整理后就用红书上的模板了。

#include <bits/stdc++.h>

using namespace std;

typedef long long ll;

ll gcd(ll a,ll b)
{
    return b ==  ? a : gcd(b,a%b);
}

ll a[];

int tot;
ll primes[];
int num_primes[];
void divide(ll x)
{
    ll tmp = (int)(double(sqrt(x))+);
    tot = ;
    ll now = x;

    for(ll i=; i<=tmp; i++) {
        if(now%i==)
        {
            primes[++tot] = i;
            num_primes[tot] = ;
            while(now%i==)
            {
                ++num_primes[tot];
                now/=i;
            }
        }
    }

    if(now!=)
    {
        primes[++tot] = now;
        num_primes[tot] = ;
    }

}

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

    for(int i=; i<n; i++)
        scanf("%lld",&a[i]);

    divide(a[]);

    for(int i=; i<n; i++)
    {
        ll g = gcd(a[],a[i]);
        ll tmp = g;
        for(int j=; j<=tot; j++)
        {
            if(tmp%primes[j]==)
            {
                tmp = tmp / primes[j];
                break;
            }
        }
        if(tmp==g)
            printf("%d ",-);
        else printf("%lld ",tmp);
    }

    return ;
}