D e s c r i p t i o n Description Description
给定一个长度为 n n n的序列 A A A,你需要将其分成两个部分,使得两部分的积互质
数据范围: n ≤ 1 0 5 n\leq 10^5 n≤105
S o l u t i o n Solution Solution
首先如果两个数的 g c d ≠ 1 gcd\neq 1 gcd=1,则它们一定不能在同一部分,将它们连边,最后求联通快数量。
这样做的复杂度是 O ( n 2 ) O(n^2) O(n2)的,我们考虑优化
对于一个数,如果它们含有相同的质因子,则它们一定不互质,所以我们只需要处理最小质因子即可
于是我们在筛素数的过程中顺便更新即可,时间复杂度近似于 O ( n + 1 0 6 ) O(n+10^6) O(n+106)
C o d e Code Code
#include<cctype>
#include<cstdio>
#include<cstring>
#define LL long long
#define mod 1000000007
using namespace std;int n,T,x,f[1000010],vis[1000010],last,tot,m;
inline int find(int x){return x==f[x]?x:f[x]=find(f[x]);}
bool v[1000010];
inline LL ksm(LL x,LL y)
{
LL res=1;
for(;y;y>>=1,(x*=x)%=mod) if(y&1) (res*=x)%=mod;
return res;
}
inline LL read()
{
char c;LL d=1,f=0;
while(c=getchar(),!isdigit(c)) if(c=='-') d=-1;f=(f<<3)+(f<<1)+c-48;
while(c=getchar(),isdigit(c)) f=(f<<3)+(f<<1)+c-48;
return d*f;
}
signed main()
{
T=read();
while(T--)
{
n=read();memset(vis,0,sizeof(vis));memset(v,0,sizeof(v));tot=0;m=0;
for(register int i=1;i<=n;i++)
{
x=read();
if(x==1) tot++;
if(x==1||vis[x]) continue;
vis[x]=++m;//存下来
f[m]=m;//并查集初始化
}
for(register int i=2;i<=1e6;i++)
{
if(v[i]) continue;last=0;
for(register int j=i;j<=1e6;j+=i)
{
v[j]=true;
if(vis[j])
{
if(vis[last]) f[find(vis[last])]=find(vis[j]);//相连
last=j;
}
}
}
for(register int i=1;i<=m;i++) if(find(i)==i) tot++;
printf("%lld\n",(ksm(2,tot)-2+mod)%mod);
}
}