HOJ 1016 Joseph's problem I 约瑟夫环问题

常规的约瑟夫环问题描述是：n个数字，每次删除第m个，最后剩下的是哪个？

它有如下的递推关系：

f(n,m)=0 n=1

f(n,m) =(f(n-1,m)+m)%n,  n>1

所以这里这道题m是素数，因此多一个步骤就是求出指定数量的素数，然后用递推关系计算即可。代码如下：

#include <iostream>
#include <stdio.h>
#define N 32620
#define Size (N/6*2 + (N%6 == 5? 2: (N%6>0)))
using namespace std;


int prime[3501];//提前算出32620以内有3501个素数
void all_prime_number();//寻找所有素数，将其存到prime数组中
int main()
{
    int i,a,fx;


    all_prime_number();
    while(cin>>a && a)
    {


        fx=0;
        for(i=0;i<a-1;i++)//逆向思考，从最初只有一个人计算到N个人；正向考虑序号会变，所以处理起来很麻烦
        {
            fx=(fx+prime[a-2-i])%(i+2);
        }
        cout<<fx+1<<endl;;
    }




    return 0;
}


void all_prime_number()
{
    int m,n,j,len,stp,k;//int型每一位表示是否为素数
    bool P[Size]={false};
    int gab = 4;


    for (k = 5; k * k <= N; k += gab ^= 6)
    {
        m=k/6*2 + ((k%6)>>2);
        n=m>>1;
        if(m%2)//6n+5
        {
            for(j=12*n*n+20*n+8,len=4*n+3,stp=12*n+10;(j+len)<Size ;j+=stp)
            {
                P[j] = 1;
                P[j+len] = 1;
            }
            if(j<Size)
                P[j] = 1;
        }
        else//6n+1
        {


            for(j=12*n*n+4*n,len=8*n+1,stp=12*n+2;(j+len)<Size;j+=stp)
            {
                P[j] = 1;
                P[j+len] = 1;
            }
            if(j<Size)
                P[j] = 1;
        }
        //cout<<k<<" ";
    }
    //cout<<(P[0] & 4095)<<endl;
    n=1;
    for(j=0;j<Size;j++)
    {


        if(!P[j])
        {
            prime[n]=j/2*6+((j%2)<<2)+1;
            //cout<<prime[n]<<endl;
            n++;
        }
    }
    prime[0]=2;
    prime[1]=3;
}