HDU 5072 Coprime

题目链接：Coprime

解题思路：这个题目最先没什么想法，看了别人的说和单色三角形有关。

单色三角形问题：平面内有N个点，每两个点都有一条边，每条边不是红色就是黑色，让你求出全是黑色的和全是红色的三角形的个数。我们从N个点选出三角形的个数是C(n,3)。我们只要减去那些不是单色的就好。怎么减去？对于每一个点来说，有a条红色边，b条黑色边。我们只要对于每个点计算a*b再加起来/2就得到结果。/2是因为一个三角形如果不是单色的话他会有两个顶点的边的颜色不同。一个顶点的边的颜色相同。所以每一个都计算了两遍。答案就是C(n, 2) - a * b / 2

回到我们这个题目，我们要求的就是全部互素或者全部不互素。那么就相当于红黑边。问题是对于每一个数怎么计算其他的数有多少和它是不互素的，只要算出了这个就可以得到互素的值那就用N-1再减去前面的值。暴力方法O(n^2)，不行。对于每个数分解质因数（每种质因数只保留一个），然后我们对于分解的质因数再进行组合。我们可以得到10^5以内的数分解质因数的种类不超过7个。因为2*3*5*7*11*13*17>100000。所以进行再组合不超过127个。举例

1680可以分解为2 3 5 7，这四种质因数，再组合就会出现6 10 14 15 21 35 30 42 70 105 210这几个因数。那么我们对于这些数制作一个统计量表示N个数字里面能被这些数整除的数有多少个。

容斥原理：当我们计算24和12的时候

24:2 3 6

12:2 3 6

15:3 5 15

18:2 3 6

所以我们得到

因数2 3 5 6 15 

个数3 4 1 3 1

那么对于24来说，和他有2这个公共因数的一共有3-1个，和他有3这个公共因数的一共有4-1个，和他有6这个公共因数的一共有3-1个

那么我们将是2的和是3的加起来发现得到的是5，其中12， 18被计算两次，所以要再减去6的2个，得到3.

这就是熔池定理，公共因数为奇数就要加上，偶数就要减去。这样就求出来了。

我这个复杂度可能有点高。。。。。。。跑了700+ms最快几十毫秒。。。。。。

中间结果可能会爆掉。记得用longlong

#include<cstdio>
#include<iostream>
#include<cmath>
#define MAX 100010
#define MAX_N 150
#define ll __int64



int seq[MAX], di[MAX][MAX_N], sum[MAX], tot[MAX];
int n;

void dfs(int head, int now, int index, int sgn, int siz){
    if(index >= siz){
        return;
    }
    for(int i = index; i < siz; i++){
        if(now != 1){
            di[head][tot[head]++] = now * di[head][i] * sgn;
        }
        dfs(head, now * di[head][i], i + 1, sgn * -1, siz);
    }
}


void factor(int x){
    int temp, now, i;
    temp = (int)((double)sqrt(x) + 1);
    tot[x] = 0;
    now = x;
    for(i = 2; i <= temp; i++){
        if(now % i == 0){
            di[x][tot[x]++] = i;
            while(now % i == 0){
                now /= i;
            }
        }
    }
    if(now != 1){
        di[x][tot[x]++] = now;
    }
}

ll calc(){
    int i, j, k;
    ll ret = 0, eve;
    for(i = 0; i < n; i++){
        eve = 0;
        for(j = 0; j < tot[seq[i]]; j++){
            int tem = di[seq[i]][j];
            if(tem < 0){
                tem *= -1;
                eve -= sum[tem] - 1;
            }
            else{
                eve += sum[tem] - 1;
            }
        }
        ret += eve * (n - eve - 1);
    }
    return ret;
}

void init(){
    memset(sum, 0, sizeof(sum));
}

__int64 get(int x){
    return 1LL * (x - 2) * (x - 1) / 2 * x / 3;
}

int main(){
    int i, j, k;
    int t;
    scanf("%d", &t);
    while(t--){
        init();
        scanf("%d", &n);
        for(i = 0; i < n; i++){
            scanf("%d", &seq[i]);
            factor(seq[i]);
            dfs(seq[i], 1, 0, 1, tot[seq[i]]);
            for(j = 0; j < tot[seq[i]]; j++){
                sum[abs(di[seq[i]][j])]++;
            }
        }
        ll ans = calc();
        printf("%I64d\n", get(n) - ans / 2);
    }
    return 0;
}