BeginOJ 入门OJ : Problem 1259: 约数最多 - DFS 非暴搜打表 - AC代码 - c/c++

标签 : BeginOJ 入门OJ : Problem 1259: 约数最多 - DFS 非暴搜打表 - AC代码标签 : BeginOJ 入门OJ : Problem 1259: 约数最多 - DFS 非暴搜打表 - AC代码

---

最近看见有些人在找这些题,这里分享一下AC代码
--------呵呵呵被你发现了--------
不想看就点这里吧

原题目:

---

Problem X: 约数最多

Time Limit: 1 Sec   Memory Limit: 128 MB
Submit: 209   Solved: 44
[ Submit][ Status][ Web Board]

Description

数学家们喜欢各种类型的有奇怪特性的数。例如，他们认为945是一个有趣的数，因为它是第一个所有约数之和大于本身的奇数。为了帮助他们寻找有趣的数，你将写一个程序扫描一定范围内的数，并确定在次范围内约数格式最多的那个数。 
不幸的是，这个数和给定的范围的都比较大，用简单的方法寻找可能需要较多的运行时间。所以请确定你的算法能在1秒内完成最大范围内的扫描。 

Input

只有一行，给出扫描的范围，由下界L和上界U确定，满足2<=L<=U<=1000000000且0<=U-L<=1000000.

Output

对于给定的范围，输出该范围内约数的个数D最多的数P.。若有多个，则输出最小的那个。请输出“Between L and U,P has a maximum of D divisors.”，其中L,U,P和D的含义同前面所述。

Sample Input

输入1
1 10
输入2
1000 1000

Sample Output

输出1
Between 1 and 10, 6 has a maximum of 4 divisors.
输出2
Between 1000 and 1000, 1000 has a maximum of 16 divisors.

HINT

[ Submit][ Status]

---

DFS完美AC代码分享:

--------呵呵呵又被你发现了--------

c++:

#include<bits/stdc++.h>
using namespace std;
int num,ans,minn,anx,yes[100001],no[100001];
void dfs(int start,int s,int x,int L,int R) {
	if (x>=minn&&((s>ans) || (s==ans && x<anx))) ans=s,anx=x;
	if ((L==R) && (L>x)) dfs(start,s<<1,x*L,1,1);
	for (int i=start; i<num; i++) {
		if (yes[i]>R) return;
		int j=yes[i],ll=L-1,l=L,r=R,y=x,ss=s,m=1;
		while (1) {
			m++,ss+=s,ll/=j,l/=j,r/=j;
			if (ll==r) break;
			y*=j;
			dfs(i+1,ss,y,l,r);
		}
		m=1<<m;
		if (s<(ans/m)) return;
	}
}
int main() {
	int l,r;
	cin>>l>>r;
	no[0]=no[1]=1;
	for (int i=2; i<100001; i++) {
		if (!no[i]) yes[num++]=i;
		for (int j=0; j<num && i*yes[j]<100001; j++) {
			no[i*yes[j]]=1;
			if (!(i%yes[j])) break;
		}
	}
	ans=anx=1;
	if (l!=1 || r!=1){
		minn=l,ans=2,anx=l;
		dfs(0,1,1,l,r);
	}
	printf("Between %d and %d, %d has a maximum of %d divisors.\n",l,r,anx,ans);
	return 0;
}

c:

#include<bits/stdc++.h>
using namespace std;
int num,ans,minn,anx,yes[100001],no[100001];
void dfs(int start,int s,int x,int L,int R) {
	if (x>=minn&&((s>ans) || (s==ans && x<anx))) ans=s,anx=x;
	if ((L==R) && (L>x)) dfs(start,s<<1,x*L,1,1);
	for (int i=start; i<num; i++) {
		if (yes[i]>R) return;
		int j=yes[i],ll=L-1,l=L,r=R,y=x,ss=s,m=1;
		while (1) {
			m++,ss+=s,ll/=j,l/=j,r/=j;
			if (ll==r) break;
			y*=j;
			dfs(i+1,ss,y,l,r);
		}
		m=1<<m;
		if (s<(ans/m)) return;
	}
}
int main() {
	int l,r;
	scanf("%d%d",&l,&r);
	no[0]=no[1]=1;
	for (int i=2; i<100001; i++) {
		if (!no[i]) yes[num++]=i;
		for (int j=0; j<num && i*yes[j]<100001; j++) {
			no[i*yes[j]]=1;
			if (!(i%yes[j])) break;
		}
	}
	ans=anx=1;
	if (l!=1 || r!=1){
		minn=l,ans=2,anx=l;
		dfs(0,1,1,l,r);
	}
	printf("Between %d and %d, %d has a maximum of %d divisors.\n",l,r,anx,ans);
	return 0;
}