hust2015暑假集训 0715 c a coprime

http://acm.hust.edu.cn/vjudge/contest/view.action?cid=82839#problem/C

14年鞍山区域赛的题...

记得是当时全场过得第四多的题...

讲题的时候丽洁姐说“莫比乌斯反演blablabla”

当时完全不懂==

嗯，第一次接触了莫比乌斯反演这个东西。。。好像这个东西在数论中挺重要的。

这道题需要注意的是，求和x的互质的个数的时候，1和任何正整数都是互质的，但是1和1互质这一组是不符合条件的....记得判断下。

求Mobius函数好像有一个筛法求（见贾志鹏线性筛.pdf），然后红书上有一个递推求1~n的mobius函数值的方法（见红书23）；

然后被强制转化类型坑了2发WA。。。要开long long ,不写强制转换就会WA。。。还没搞明白为什么。看了别人写的代码上有，加上，结果就过了==

知道了。感谢zj学长。 

 等号右边如果是int，左边是LL，超int了也会越界  

就是说c=a*b  如果a 和b 是int  C是long long  如果a*b超过了Int,就会越界...

需要先把两个int类型转化成long long  写成c=long long (a)*long long (b);

 

/*******************************************************************

Author :111qqz

*******************************************************************/



#include <algorithm>

#include <cstdio>

#include <iostream>

#include <cstring>

#include <string>

#include <cmath>

#include <map>

#include <stack>

#include <queue>

#include <set>



using namespace std;

typedef long long LL;

typedef unsigned long long ULL;

const int inf = 0x7fffffff;

const int N=1E5+3;

int T,n;

int mu[N],a[N];

LL mul[N],p[N],sumt[N];





void getmu()

{

    for (int i = 1 ; i <N ; i++)

    {

        int target = i ==1?1:0;

        int delta = target - mu[i];

        mu[i] = delta;

        for ( int j = i+i; j <N ; j=j+i)

            mu[j] = mu[j] + delta;

    }



}

int main()

{

    cin>>T;

    getmu();

    while (T--)

    {

         cin>>n;

       // cout<<mu[n]<<endl;

        memset(p,0,sizeof(p));

        memset(sumt,0,sizeof(sumt));

        memset(mul,0,sizeof(mul));

        memset(a,0,sizeof(a));



        for ( int i = 0 ;  i < n ; i++ )

        {

            scanf("%d",&a[i]);

            sumt[a[i]]++;

        }

        for ( int i = 1  ; i < N ; i++)

            for ( int j = i ;  j < N ; j=j+i)

            {

                mul[i] = mul[i] + sumt[j];

            }



        for ( int i = 1 ; i < N ; i++ )

            for ( int  j = i ;  j < N ; j=j+i)

            {

                p[j] = p[j] + mul[i]*mu[i];

            }

       // p[1]--;//减到1和1本身互质的那个...因为题目中说所有人的id不相同了

       // cout<<"p[1]:"<<p[1]<<endl;

      //  for ( int i = 0 ;  i < n ; i++)

        //    cout<<"p[i]："<<a[i]<<" "<<p[a[i]]<<endl;



        LL ans=(LL)n*(n-1)*(n-2)/6;

        LL sum=0;

        for (int i=0;i<n;i++)

        {

            if (a[i]==1) continue;

            sum+=(LL)p[a[i]]*(LL)(n-1-p[a[i]]);



        }

        sum =(LL)sum/2;

        ans =(LL)(ans-sum);

        cout<<ans<<endl;



    }









    return 0;

}