数论题目整理（持续补充）

数论错题整理

1、最大的不能被表示成素数幂的数+最小公倍数求法

题目链接

对于2，不止有因子2还有其他素数因子，能达到的最大数为$3\ast 2^k$（只有乘以3，k才有可能最大）

对于x>2的数，不止有因子x，能达到的最大数为$2\ast x^k$

题目要求我们求出小于n的，不能被表示成素数的整数次幂的数，的最小公倍数

由算数基本定理，$lcm=p{1}^{max(a1,b1)}p{2}^{max(a2,b2)}p{3}^{max(a3,a3)}....$

我们只需要求出最大的幂次即可

#include<bits/stdc++.h>
using namespace std;
const int N=8e7+5;
const int mod=1e9+7;
typedef long long ll;
int primes[N],cnt;
bool st[N];
void get_primes(int n){
	for(int i=2;i<=n;i++){
		if(!st[i]){
			primes[cnt++]=i;
		}
		for(int j=0;primes[j]<=n/i;j++){
			st[primes[j]*i]=true;
			if(i%primes[j]==0) break;
		}
	}
}
int main(){
	int n;
	int i;
	scanf("%d",&n);
	if(n<6){
		puts("empty");
		return 0;
	}
	ll res=1;
	get_primes(N);//先筛出8e7+5以内的素数
	for(i=1;primes[i]<=n/2;i++){
		ll p=primes[i];
		ll tmp=1;
		while(tmp*p<=n/2) tmp*=p;//看看2*这个素数最高能到最少次幂
		res=res*tmp%mod;
	}
	ll tmp2=1;
	while(tmp2*2<=n/3) tmp2*=2;//看看3*2最高能到最少次幂
	res=res*tmp2%mod;
	cout<<res;
}

2、构造gcd>1的数对

题目链接

题眼是k<n/2，因为从1到n中，偶数恰好为n/2个，所以

• 若k<n/2，直接摆放2到n/2中的偶数在前面，后面依次填写1,3,5，…
• 若k==n/2，因为所有的偶数放在前面不足以构成n/2个数对，所以拿出一个6放在最后，和奇数中的3凑成一对

注意特判n<=5的情况，n=5,k=2,n=4,k=2都是不存在解的

#include<bits/stdc++.h>
using namespace std;
int main(){
	int n,k;
	cin>>n>>k;
	if(n<6 && k==n/2){
		puts("-1");
		return 0;
	}
	if(k<n/2){
		int i=2;
		for(int j=1;i<=n&&j<=k+1;i+=2,j++){
			cout<<i<<" ";
		}
		i-=2;
		for(int j=i+1;j<=n;j++){
			cout<<j<<" ";
		}
		for(int j=1;j<=i;j+=2){
			cout<<j<<" ";
		}
	}else{
		int i=8;
		for(int j=1;j<=k-2&&i<=n;j++,i+=2){
			cout<<i<<" ";
		}
		for(int j=2;j<=6;j+=2){
			cout<<j<<" ";
		}
		cout<<3<<" "<<1<<" ";
		for(int j=5;j<=n;j+=2){
			cout<<j<<" ";
		}
	}
}

3、求整个数组、超大数的最大公约数

试题链接：2021年牛客寒假集训营第四场–J题

#include<cstdio>
#include<algorithm>
using namespace std;
typedef long long ll;
const int N=1e4+5;
const int mod=1e9+7;
int x[N],p[N];
int main(){
	int n;
	scanf("%d",&n);
	for(int i=0;i<n;i++) scanf("%d",&x[i]);
	for(int i=0;i<n;i++) scanf("%d",&p[i]);
	ll res=1;
	for(int i=2;i<=10000;i++){//ö¾ÙÒò×Ó 
		int cnt=0x3f3f3f3f;//ÿ¸öÊý¡¢Õâ¸öÒò×ÓµÄ×îС´ÎÃÝ 
		for(int j=0;j<n;j++){
			int s=0;
			while(x[j]%i==0){
				x[j]/=i;
				s++;
			}
			cnt=min(s*p[j],cnt);
		}
		for(int k=1;k<=cnt;k++){
			res*=i;
			res%=mod;
		}
	}
	printf("%d\n",res);
}