codeforces round 911 D - Small GCD

最新推荐文章于 2024-10-02 23:40:03 发布

Heredy

最新推荐文章于 2024-10-02 23:40:03 发布

阅读量978

点赞数 31

分类专栏： codeforces题解文章标签：算法数据结构 c++

本文链接：https://blog.csdn.net/wyzz_yz/article/details/134731120

版权

codeforces题解专栏收录该内容

32 篇文章 0 订阅

订阅专栏

分析

$\sum\limits_{i=1}^{n}\sum\limits_{j=i+1}^{n}\sum\limits_{k=j+1}^{n}=f(a_i,a_j,a_k)\quad a_i,a_j$ 是较小的两数
易得简单做法排序后通过枚举前两个数结合 $GC D$ 计算，但是复杂的为 $n^2logn$ 会 $TE L$ 。

思路

因数分解：预处理出可能会用到的所有数的所有因数。因数分解比较耗时，所以在每个 $c a se$ 里再去处理是不现实的。
试想在排序后我们枚举 $a_i\;a_j$ 计算两者 $GC D$ 的贡献，计算的就是因数对答案的贡献~~不过是最大公约数~~ 。我们可以用一个 $f or$ 遍历数组，每到一个数就记录该数所有出现过的因数对答案的贡献~~最大公约数对答案才有贡献，但我们不计算最大公约数，我们记录每一个约数的贡献，其中一定包括最大公约数，而出现过就保证了是公约数~~ $sum[x]{+}{=}cnt[x]*(n-i)\;n-i$ 是距离，应题目。记录之后需 ${++}cnt[x]$ 记录出现过的次数
容斥：考虑 $x$ 就是最大公约数，那所有 $x$ 的倍数作为最大公约数时的贡献都应该减去。因为 $x$ 与 $x$ 的倍数是同在一个因数分解的桶中的，而我们 $f or$ 遍历时 ${++}cnt[x]$ 每一个因数。
Think Twice, Code Once

#include<bits/stdc++.h>
#define il inline
#define get getchar
#define put putchar
#define is isdigit
#define re register
#define int long long
#define dfor(i,a,b) for(re int i=a;i<=b;++i)
#define dforr(i,a,b) for(re int i=a;i>=b;--i)
#define dforn(i,a,b) for(re int i=a;i<=b;++i,put(10))
#define mem(a,b) memset(a,b,sizeof a)
#define memc(a,b) memcpy(a,b,sizeof a)
#define pr 114514191981
#define gg(a) cout<<a,put(32)
#define INF 0x7fffffff
#define tt(x) cout<<x<<'\n'
#define ls i<<1
#define rs i<<1|1
#define la(r) tr[r].ch[0]
#define ra(r) tr[r].ch[1]
#define lowbit(x) (x&-x)
using namespace std;
typedef unsigned int ull;
int read(void)
{
    re int x=0,f=1;re char c=get();
    while(!is(c)) (f=c==45?-1:1),c=get();
    while(is(c)) x=(x<<1)+(x<<3)+(c^48),c=get();
    return x*f;
}
void write(int x)
{
    if(x<0) x=-x,put(45);
    if(x>9) write(x/10);
    put((x%10)^48);
}
#define writeln(a) write(a),put(10)
#define writesp(a) write(a),put(32)
#define writessp(a) put(32),write(a)
const int N=1e5+10,M=3e4+10,SN=1e4+10,mod=19930726;
int n,ans,a[N],cnt[N],sum[N];
vector<int > ve[N];
signed main()
{
    int T=read();
    dfor(i,1,N)
        for(re int j=i;j<N;j+=i) ve[j].push_back(i);
    while(T--)
    {
        mem(cnt,0),mem(sum,0),ans=0;
        n=read();
        dfor(i,1,n) a[i]=read();
        sort(a+1,a+n+1);
        dfor(i,1,n)
        {
            int len=ve[a[i]].size()-1;
            dfor(j,0,len)
            {
                sum[ve[a[i]][j]]+=cnt[ve[a[i]][j]]*(n-i);
                ++cnt[ve[a[i]][j]];
            }
        }
        dforr(i,a[n],1)
        {
            for(re int j=i<<1;j<N;j+=i) sum[i]-=sum[j];
            ans+=sum[i]*i;
        }
        writeln(ans);
    }
    return 0;
}