【2019/07/20测试 T1】最大公约数

最新推荐文章于 2023-03-30 17:48:53 发布

forever_dreams

最新推荐文章于 2023-03-30 17:48:53 发布

阅读量253

点赞数

分类专栏： --------数学问题-------- # 快速幂文章标签：数学问题

本文链接：https://blog.csdn.net/forever_dreams/article/details/96631266

版权

--------数学问题-------- 同时被 2 个专栏收录

141 篇文章 1 订阅

订阅专栏

快速幂

15 篇文章 0 订阅

订阅专栏

传送门

Problem

给出 $n$ 个数，计算其所有子集的 $G C D$ （即子集所有数的 $g c d$ ）的乘积对 $10^9+7$ 取模的结果。

其中空子集的 $G C D$ 为 $1$ 。

数据范围： $2≤n,a_i≤10^5$ 。

Solution

考场上一个细节没注意，硬生生把 $100$ 分变成了 $20$ 分。。。

先用桶排，算出每个数的出现次数。然后考虑枚举 $g c d$ ，假设当前枚举的数为 $x$ ，可以求出 $x$ 的倍数个数 $sum_x$ ，那么初步认为以 $x$ 为 $g c d$ 的集合数 $set_x=(2^{sum_x}-1)$ 个。

但是这样显然是不对的，举个例子，在枚举 $2$ 的时候，这样会计算到 ${6\},\{4,8\}$ 等集合。要想办法除去这些集合。

由于多算的集合的 $g c d$ 肯定是 $x$ 的倍数，因此我们枚举 $x$ 的倍数然后减掉相应的 $s e t$ 值即可，即以 $x$ 为 $g c d$ 的子集数量为 $set_x-\sum_{n|d}set_d$ 。

由于我们要先知道所有倍数的 $s e t$ 值才能推到 $x$ ，故应倒着循环。

最后的答案 $ans=\prod_{i=1}^{10^5}i^{set_i}$ ，时间复杂度 $O(n\log n)$ 。

有一个地方要注意，即 $set_i$ 可能比较大，我们要先对它进行降幂操作。根据费马小定律，当 $p$ 为质数时， $a^{p-1}\equiv 1(\,\mathrm{mod}\;p)$ ，因此要将 $set_i\%(p-1)$ 再用快速幂。

Code

#include<cstdio>
#include<cstring>
#include<algorithm>
#define N 100005
#define Mod 1000000007
using namespace std;
int cnt[N],Set[N];
int Power(int a,int b,int p){
	int ans=1;
	for(;b;b>>=1,a=1ll*a*a%p)
		if(b&1)  ans=1ll*ans*a%p;
	return ans;
}
int main(){
	int n,x,ans=1;
	scanf("%d",&n);
	for(int i=1;i<=n;++i)
		scanf("%d",&x),cnt[x]++;
	for(int i=2;i<=1e5;++i){
		int sum=0;
		for(int j=i;j<=1e5;j+=i)  sum+=cnt[j];
		Set[i]=Power(2,sum,Mod-1)-1;
	}
	for(int i=1e5;i>=2;--i){
		if(!Set[i])  continue;
		for(int j=2*i;j<=1e5;j+=i)  Set[i]=(Set[i]-Set[j]+Mod-1)%(Mod-1);
		ans=1ll*ans*Power(i,Set[i],Mod)%Mod;
	}
	printf("%d",ans);
	return 0;
}