呃 等重现放出来再说吧。。。
#include <bits/stdc++.h>
inline int gcd(int a,int b)
{
return a==0?b:gcd(b%a,a);
}
#define N 1e8
using namespace std;
int n,m;
long long xx;
long long a[1010];
int fai[10000];
int check[10000];
int p[10000];
vector<long long> vec;
void getfai()
{
fai[1]=1;
int tot=0;
for(int i=2; i<=10000; i++)
{
if(!check[i])
{
p[tot++]=i;
fai[i]=i-1;
}
for(int j=0; j<tot; j++)
{
if(i*p[j]>N) break;
check[i*p[j]]=1;
if(i%p[j]==0)
{
fai[i*p[j]]=fai[i]*p[j];
break;
}
else
{
fai[i*p[j]]=fai[i]*(p[j]-1);
}
}
}
}
int main()
{
//freopen("in.txt","r",stdin);
int t;
scanf("%d",&t);
getfai();
while(t--)
{
vec.clear();
scanf("%d %lld",&n,&m);
for(int i=1; i<=n; i++)
{
scanf("%d",&xx);
a[i]=(long long)gcd(xx,m);
}
sort(a+1,a+n+1);
n=unique(a+1,a+n+1)-a-1;
for(int i=2; i*i<=m; i++)
{
if(m%i==0)
{
vec.push_back(i);
if(i*i!=m) vec.push_back(m/i);
}
}
long long re=0;
vector<long long>::iterator it;
for(it=vec.begin(); it!=vec.end(); it++)
for(int i=1; i<=n; i++)
{
long long x=*it;
if(x%a[i]==0)
{
re+=m*fai[m/x]/2;
break;
}
}
printf("%lld\n",re);
}
return 0;
}