LG P5655 基础数论函数练习题（GCD）

最新推荐文章于 2022-08-21 18:39:33 发布

Freopen

最新推荐文章于 2022-08-21 18:39:33 发布

阅读量264

点赞数

分类专栏：数论性质分析

本文链接：https://blog.csdn.net/qq_35950004/article/details/103064287

版权

数论同时被 2 个专栏收录

157 篇文章 0 订阅

订阅专栏

性质分析

86 篇文章 0 订阅

订阅专栏

题目

题解：
考虑如何求 $n$ 个数的 $l c m$ ，可以求出前 $n - 1$ 个数的 $l c m = x$ ，再和 $a_n =y$ 求 $\frac {xy}{gcd(x,y)}$
~~然后就可以python骗到50分了~~
$\frac {xy}{\gcd(x,y)} = \frac {a_n}{\gcd(x,y)} x$
$x$ 可以模。
考虑怎么在不用高精度的情况下计算 $\frac {y}{\gcd(x,y)} \mod 1000000007$
$\gcd(x,y) = \gcd(x \mod y , y)$
所以我们只需要再算 $x\mod y$ 就行了。
可以 $O(n^2+n\log n)$ 求出 $n$ 个数的 $l c m$

考虑多组询问。
我们考虑预处理出所有答案。
假设我们安排了在求 $l c m$ 中一个一个加入 $a_i$ 的顺序，
那么放 $a_i$ 进来的时候，可以直接令 $b_i = \frac {a_i}{\gcd(a_i,x)}$
那么答案就相当于是 $\prod b_i$
对于所有的区间，我们都安排一个统一的顺序（从大编号到小编号加入），那么我们可以寻求到一些可以重复利用的计算信息。
从小标号到大插入 $a_i$ ，然后对于每个 $b_j,j<i$ 更新为 $\frac {b_j}{gcd(b_j,\prod _{k=j+1}^{i} b_k)}$ ，其实只需要看 $a_i$ 就可以了。
然后（对于右端点是 $i$ 的）区间询问 $[l, i]$ 就是 $\prod_{j=l}^ib_j$
这个方法每次都要求 $\gcd$ ，是 $O(n^2\log n)$ ,总复杂度是 $O(Tn^2\log n)$ 难以通过。
~~O(1)gcd需要O(值域)的预处理~~

考虑优化每次求 $\gcd$ 的过程。
一个数 $a_i$ 和前面 $i - 1$ 段 $[l, i - 1]$ 分别的 $\gcd$ 只会有 $O(\log a)$ 种。
所以从小编号到大编号更新 $b_j$ ,用模运算来快速判断 $\gcd$ 是否会变小（一段变短乘积变小）以此计算。
复杂度是 $O(n^2 + n\log ^2n)$ 可以通过本题。

#include<bits/stdc++.h>
#define maxn 305
#define mod 1000000007
#define LL long long
using namespace std;

char cb[1<<15],*cs=cb,*ct=cb;
#define getc() (cs==ct&&(ct=(cs=cb)+fread(cb,1,1<<15,stdin),cs==ct)?0:*cs++)
template<class T>void read(T &res){
	char ch;
	for(;!isdigit(ch=getc()););
	for(res=ch-'0';isdigit(ch=getc());res=res*10+ch-'0');
}

int n,q,ans[maxn][maxn];
LL a[maxn],b[maxn];
inline LL mul(LL a,LL b,LL p){
	LL r = a*b-(LL)((long double)a/p*b+0.5)*p;
	return r < 0 ? r + p : r;
}
inline LL gcd(LL a, LL b){
    if(!a||!b)return a|b;
    register int t=__builtin_ctzll(a|b);
    a>>=__builtin_ctzll(a);
    do{
        b>>=__builtin_ctzll(b) ;
        if(a>b){LL t=b;b=a,a=t;}
        b-=a;
    }while(b);
    return a<<t;
}

int main(){
	int T;
	for(read(T);T--;){
		read(n),read(q);
		for(int i=1;i<=n;i++){
			read(a[i]);
			b[i] = 1;
			for(int j=i-1;j>=1;j--) b[j] = mul(a[j],b[j+1],a[i]);
			LL G = gcd(b[1] , a[i]) , t;
			for(int j=1;j<i;j++) 
				if(t=(b[j+1]%G)){
					t=gcd(t,G);
					a[j] /= G / t , G = t;
				}
			ans[i][i] = a[i] % mod;
			for(int j=i-1;j>=1;j--) ans[i][j] = 1ll * ans[i][j+1] * (a[j] % mod) % mod;
		}
		for(int l,r;q--;){
			read(l),read(r);
			printf("%d\n",(ans[r][l]+mod)%mod);
		}
	}
}