六校联合周赛F题

最新推荐文章于 2018-11-22 21:22:35 发布

taozitaozi99

最新推荐文章于 2018-11-22 21:22:35 发布

阅读量181

点赞数 1

本文链接：https://blog.csdn.net/qq_42398530/article/details/83547991

版权

题目：

Problem Description

给定n个整数a1,a2,a3……an，且1<=ai<=n问有多少对i,j 满足i<j 且gcd(a[i],a[j])=1,2,3,……n。

Input

单组数据

第一行包含一个整数n。

接下一行有n个整数 a1,a2,a3……an

1<=ai<=n

1<=n<=105

Output

输出n行，每行包含一个整数，第k行代表gcd(a[i],a[j])=k的方案数

SampleInput
 
10 
1 1 1 1 1 1 1 1 1 1
SampleOutput
 
45 0 0 0 0 0 0 0 0 0

题目大意：

简而言之，输入n个数，输出在第k行有多少对 gcd(最大公约数)=k。

题目思路：

设F[i]表示 gcd=i的倍数的对数

设G[i]表示gcd=i 的对数，则结果等于F[i]-G[2*i]-G[3*i]-……

b[i]表示i的倍数有多少个，F[i]=b[i]*(b[i]-1)/2

代码：(整体容斥)

#include<cstdio>
#include<cstring>
#include<stdlib.h>
#include<algorithm>
#include<time.h>
#include<vector>
using namespace std;
long long a[1000005];
int main(){
    int n,i,x,y;
    scanf("%d",&n);

    for(i=1;i<=n;i++){
        scanf("%d",&x);
        a[x]++;
    }
    register int j;
    for(i=1;i<=n;i++)
    {
        for(j=i+i;j<=n;j+=i)
        {
            a[i]+=a[j];
        }
        a[i]=a[i]*(a[i]-1)/2;
    }
    for(i=n;i>=1;i--)
    {
        for(j=i+i;j<=n;j+=i)
        {
            a[i]-=a[j];
        }
    }
    for(i=1;i<=n;i++)
    {
        printf("%lld\n",a[i]);
    }
    return 0;
}