Dream Team（最小生成树）

Input file: dream.in

Output file: standard output

Balloon Color: Pink

Before the World Cup, we would like to form the ultimate dream team. The team that will win against all teams. We have a big group of the top N players in the world, where each player has a number representing his skill level. Two players i and j with skill levels ai and aj have a compatibility of passing the ball between each other; this compatibility level is measured by the greatest common divisor of their skill levels (i.e. compatibility(i, j) = gcd(ai , aj )). We would like to decide a strategy that connects all the players of the dream team with a tree of connections between the players. If two players are directly connected in the chosen strategy, then they will pass the ball between each others during the matches, with the compatibility level between them. The compatibility of the strategy is the sum of the compatibility levels of its connections. What is the maximum total compatibility of the chosen strategy that connects all players?

Input

The first line of the input contains a single integer 1 ≤ T ≤ 50 the number of test cases. Each test case consists of one line that contains N + 1 space separated integers; the first integer is N, followed by the skill levels a1, · · · , an; where 1 ≤ N ≤ 105 and 1 ≤ ai ≤ 105 for all i.

Output

For each test case, print a single line displaying the case number and the maximum compatibility of the dream team’s strategy.

Example

dream.in

2

2 3 6

5 5 6 7 10 21

standard output

Case 1: 3

Case 2: 17

Note

In the second test case we connect the players 3−5 with compatibility 7, players 1−4 with compatibility 5, players 2−5 with compatibility 3 and players 2−4 with compatibility 2; therefore the total compatibility of the dream team’s strategy is 17.

从样例分析可知，从大到小一直找最大的最大公约数，即可得到答案

#include <bits/stdc++.h>

using namespace std;

const int maxn = 1e6 + 50;
int pre[maxn];

vector<int> ve[maxn];

int Find(int x)
{
    if(pre[x] == x) return x;
    else return pre[x] = Find(pre[x]);
}

int main()
{
    freopen("dream.in", "r", stdin);
    int t, i, j, a, b, x, k;
    long long sum;
    cin >> t;
    k = 1;
    while(t--)
    {
        int n;
        cin >> n;
        for(i=1;i<=1e5;i++)
        {
            pre[i] = i;
            ve[i].clear();
        }
        for(i=1;i<=n;i++)
        {

            cin >> x;
            for(j=1;j<=x;j++)
            {
                if(x % j == 0)
                {
                    ve[j].push_back(i);
                }
            }
        }
        sum = 0;
        for(i=1e5; i>=1;i--)
        {
            if(ve[i].empty()) continue;
            for(j=0; j<ve[i].size()-1; j++)
            {
                a = Find(ve[i][j]);
                b = Find(ve[i][j+1]);
                if(a!=b)
                {
                    pre[b] = a;
                    sum += i;
                }
            }
        }
        cout << "Case " << k++ << ": ";
        cout << sum << endl;
    }
    return 0;
}