SuperPrime Rib
题目(简述):
找出某一长度的超级素数。所谓的超级素数,例如7331:其中7331是素数,733是素数,73是素数,7是素数。即每次从该数中去掉最后一位数后,依然是素数的素数,被称为超级素数(superprime)。
分析:
判断一个数是不是超级素数,我们需要判断n次,n为这个数的位数。为了使程序更快,我们应该首先判断最高位是不是素数,然后判断最高两位是不是素数……直到判断到这个数本身是不是素数。其次,我们并不需要遍历所有的n位数,因为很多数直接就可以排除掉:除了最高位,其他任何位上有偶数的排除,1开头的数排除,至于为什么,原因很简单,大家可以自己思考。这样之后,我们其实可以生成这样的数,而不是遍历。生成这样的数的时候,注意要按照从小到大的顺序。下面是代码实现。
代码:
Executing...
Test 1: TEST OK [0.000 secs, 3048 KB]
Test 2: TEST OK [0.000 secs, 3048 KB]
Test 3: TEST OK [0.000 secs, 3048 KB]
Test 4: TEST OK [0.000 secs, 3048 KB]
Test 5: TEST OK [0.000 secs, 3048 KB]
/*
ID:haxc_xc1
PROG:sprime
LANG:C++
*/
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
int n;
bool prime(int num)
{
bool ret=true;
int i;
for(i=2;i<=sqrt(num);i++)
{
if(num%i==0)
{
ret=false;
break;
}
}
return ret;
}
void deal(int num,int k)
{
if(k==n)
{
if(prime(num)) //搜索到k位的时候,如果满足条件,则直接输出
{
printf("%d\n",num);
}
}
else
{
if(prime(num*10+1)) //如果这个数满足,则向下搜索
{
deal(num*10+1,k+1);
}
if(prime(num*10+3))
{
deal(num*10+3,k+1);
}
if(prime(num*10+5))
{
deal(num*10+5,k+1);
}
if(prime(num*10+7))
{
deal(num*10+7,k+1);
}
if(prime(num*10+9))
{
deal(num*10+9,k+1);
}
}
}
int main(void)
{
freopen("sprime.in","r",stdin);
freopen("sprime.out","w",stdout);
scanf("%d",&n);
deal(2,1);
deal(3,1);
deal(5,1);
deal(7,1);
return 0;
}
Checker Challenge
题目:
这个问题是8皇后问题的扩展,读入一个数n,代表n*n的格子,求出所有的方案数,并且输出前3种方案的具体内容。
分析:
笔者写的第一个程序,是基于深度优先搜索,没有任何优化,n=13时运算超时。后来看了若干别人的分析之后,了解到matrix67大牛的一个基于位运算的算法,深表佩服:具体链接:http://www.matrix67.com/blog/archives/266。具体内容里面已经讲述地很清楚了,也就不再赘述。代码部分附上本人所写的基于此思想的C语言的实现。另外,这个题还可以采用对称的思想对其进行优化。
代码(只输出了所有方案的数量,没有输出前3个具体的方案):
#include <stdio.h>
#include <stdlib.h>
int n;
unsigned upperlim;
int result;
void deal(unsigned row,unsigned ld,unsigned rd,int deep)
//row代表列有冲突的,ld代表左斜对角线有冲突的,rd代表右斜对角线有冲突的。deep代表层数或者深度,都可以。
{
if(row!=upperlim)
//还没有放完
{
unsigned pos=upperlim & ~(row|ld|rd),p;
while(pos)
{
p=pos&(-pos); //p为右起第一个安全位置
pos=pos^p; //设置禁止位
deal(row|p,(ld|p)>>1,(rd|p)<<1,deep+1);
}
}
else
{
result++;
}
}
int main(void)
{
freopen("checker.in","r",stdin);
freopen("checker.out","w",stdout);
scanf("%d",&n);
upperlim=(1<<n)-1; //upperlim的每一位都表示了该位置是否能够放下
deal(0,0,0,1);
printf("%d\n",result);
return 0;
}