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
,a
;
LL mul
,p
,sumt
;

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
<<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;
}