这道题的python版本居然不超时,真是惊喜了,毕竟在跑uDebug的测试数据时,感觉都过了好几秒才出的答案。。。
回到本题,题意是说给定一个全部由数字组成的字符串,字符串长度最多为256。给定的素数范围是100,000以内。要求从字符串中找出最大的素数并输出。输入时,仅有一个0组成的字符串"0"表示输入结束,输入的字符串不会超过1000个行(1行一个字符串)
题目其实还是比较简单的,只要根据题意,先构造出100,000以内的素数表(埃拉托斯特尼筛法或者欧拉筛法都可以),然后从5到1倒序取子串进行遍历,判断是否素数,是否即可。为了能够让python版本跑的更快,我增加了一下的代码,用来判断取到的子串素数是否是最大的。
if len(str(maxprime)) == i:
break
因为之前简单粗暴的版本有bug(如下,原意是想如果长子串构成的素数肯定比以后的短子串构成的素数要大,但是我错了),在uDebug的测试数据中发现了,是字符串转整型时的前导0的问题引起的,比如下面的测试用例,就会在长度为5的子串中发现一个11的素数,但不幸的是,之后还能在长度为3的子串中找到素数101,比11还大。
if maxprime != 2:
break
1101100110001100001100000110000001100000110000110001100110
11 ##错误结果
101##正确结果
python版本AC代码
maxn = 100000
visit = [True for _ in range(maxn)]
prime = [0 for _ in range(maxn)]
def init_prime():
num = 0
for i in range(maxn):
if i < 2:
visit[i] = False
continue
if visit[i] == True:
num += 1
prime[num] = i
j = 1
while j <= num and i * prime[j] < maxn:
visit[i*prime[j]] = False
if i % prime[j] == 0:
break
j += 1
init_prime()
while True:
line = input()
if line == '0':
break
i = 5
maxprime = 2
while i > 0:
if i > len(line):
i -= 1
continue
maxj = len(line)-i
j = 0
for j in range(maxj):
no = int(line[j:j+i])
if visit[no] == True:
if no > maxprime:
maxprime = no
if len(str(maxprime)) == i:
break
else:
i -= 1
print(maxprime)
C++版AC代码
#include <iostream>
#include<cstdio>
#include<cstring>
using namespace std;
//#define ZANGFONG
const int maxn = 100000;
bool visit[maxn];
int prime[maxn];
char line[256];
void init_prime()
{
memset(visit,true,sizeof(visit));
memset(prime,0,sizeof(prime));
int num = 0;
int i,j;
visit[1] = false;
visit[0] = false;
for(i = 2; i < maxn; i++)
{
if(visit[i]==true)
{
num++;
prime[num] = i;
}
for(j = 1; (j <= num)&&(i * prime[j]<=maxn); j++)
{
visit[i*prime[j]] = false;
if(i% prime[j]==0) break;
}
}
}
int main()
{
#ifdef ZANGFONG
freopen("in.txt","r",stdin);
freopen("out.txt","w",stdout);
#endif // ZANGFONG
int i,j,k,len,no;
int maxprime,flag;
init_prime();
while(scanf("%s\n",line))
{
len = strlen(line);
if(len == 1 && line[0] == '0') break;
flag = 0; maxprime = 2;
for(i = 5; i >= 1; i--)
{
if(i > len) continue;
for(j = 0; j <= len - i; j++)
{
no = 0;
for(k = j; k < i+j; k++)
{
no *= 10;
no += line[k] - '0';
}
if(visit[no] == true)
{
if(no > maxprime) maxprime = no;
}
}
}
printf("%d\n",maxprime);
memset(line,0,sizeof(line));
}
return 0;
}