反素数

反素数

其实反素数是什么并不是那么重要，重要的是它的贪心思想可以在其他类似题中使用
反素数n就是【1，n】区间内，n是因数最多的数
做反素数题就是爆搜，dfs，但有两个重要的剪枝：

1. 反素数必由从2开始的连续素数相乘而得
   原因：贪心，25 的因数和 23的因数一样多，但23比25小
2. n = 2^a1 * 3^a2 * 5^a3 *7^a4 *…中 a1>=a2>=a3>=a4>=…
   原因：也是贪心，道理也比较浅显

还有最重要的一条原理就是，n的因字数等于(1 + a1)(1 + a2)(1 +a3)*…

例题

CF 27E
求因字数为n的最小数

#include<bits/stdc++.h>
using namespace std;
typedef long long LL;

LL pri[16] = {2,3,5,7,11,13,17,19,23,29,31,37,41,43,47,53};
    //为什么就算是1e18也只需要这些数就够了呢? 因为越往后,前面的次方就越多
LL ans,n;
void dfs(LL Y,LL dep,LL ci,LL num){
        //因子数 深度 上一个数的次数 当前的数
    if(dep==16)return;
    if(Y==n){
    	ans=min(ans,num);
    	return;
	}
    for(int i=1;i<=ci;i++){
        if(ans/pri[dep]<num)break;
        dfs(Y*(i+1),dep+1,i,num*=pri[dep]);
    }
}

int main(){
    while(scanf("%lld",&n)==1){
        ans=2e18;
        dfs(1,0,64,1);
        printf("%lld\n",ans);
    }
}

P1463
求范围内最小反素数

#include<bits\stdc++.h>
using namespace std;
typedef long long LL;

LL pri[16] = {2,3,5,7,11,13,17,19,23,29,31,37,41,43,47,53};
    //为什么就算是1e18也只需要这些数就够了呢? 因为越往后,前面的次方就越多
LL ans,n,ys;
void dfs(LL Y,LL dep,LL ci,LL num){
        //因子数 深度 上一个数的次数 当前的数
    if(num>n||(dep==14))return;
    if(Y>ys){
    	ans=num;
    	ys=Y;
	}	
	else{
		if(Y==ys&&num<ans){
			ans=num;
		}
	}
    for(int i=1;i<=ci;i++){
        if(n/pri[dep]<num)break;
        dfs(Y*(i+1),dep+1,i,num*=pri[dep]);
    }
}

int main(){
scanf("%lld",&n);
        ans=0;ys=0;
        dfs(1,0,34,1);
        printf("%lld\n",ans);
    
}

P1574
这题的问题在于数据过大，询问过多，如果直接对每个询问都求反素数会超时，所以应该先打表

#include<bits/stdc++.h>
using namespace std;
typedef long long LL;

LL pri[17] = {2,3,5,7,11,13,17,19,23,29,31,37,41,43,47,53,59};
    //为什么就算是1e18也只需要这些数就够了呢? 因为越往后,前面的次方就越多
LL ans,n,ys,MAX;
LL query[100005],answer[100005];
void dfs(LL Y,LL dep,LL ci,LL num){
        //因子数 深度 上一个数的次数 当前的数
    if(num>n||(dep==17))return;
    if(Y>ys){
    	ans=num;
    	ys=Y;
	}	
	else{
		if(Y==ys&&num<ans){
			ans=num;
		}
	}
    for(int i=1;i<=ci;i++){
        if(n/pri[dep]<num)break;
        dfs(Y*(i+1),dep+1,i,num*=pri[dep]);
    }
}

int main(){
	int T;
	scanf("%d",&T);
	for(int i=1;i<=T;i++){
		scanf("%lld",&query[i]);
		MAX=max(MAX,query[i]);
	}
	int j=-1;
	while(MAX>=2){
		n = MAX;ys=0;
		dfs(1,0,66,1);
		answer[++j]=ans;
		MAX=ans-1;
	}
	for(int i=1;i<=T;i++){
		LL q=query[i];
		int l=0,r=j,mid;
		LL ans1=0;
		while(l<=r){
			mid=(l+r)/2;
			//cout<<mid<<endl;
			if(answer[mid]>q){
				l=mid+1;
			}
			else{
				ans1=answer[mid];
				r=mid-1;
			}
		}
		printf("%lld\n",ans1);
	}
    //cout<<'a'<<endl;
}