约瑟夫环的数学推导、数学方法求最后出圈的数字、循环单链表求所有出圈数字顺序

约瑟夫环是一个数学的应用问题:已知n个人(以编号1,2,3…n分别表示)围坐在一张圆桌周围;从编号为k的人开始报数,数到m的那个人出列;他的下一个人又从1开始报数,数到m的那个人又出列;依此规律重复下去,直到圆桌周围的人全部出列。

前几天,在一篇文章中得知了约瑟夫环的问题。然后,就涉及了解决办法。这个问题,在许多计算机或者关于数据结构的书中都有提及,而其中的解决办法便是使用循环链表——无论这个循环链表是使用指针还是数组实现,模拟约瑟夫环的进行,最后得到解决方案。具体方法肯定早有某人披露。但是,令人深感奇妙的还是这个问题的数学解决。在不统计解决过程,即不统计每次都需要出列哪个序号时,就可以应用数学递推公式,直接得出最后剩下的那个人的序号。

先总结一下约瑟夫环的递推公式:

  1. f[1]=0; f[i]=(f[i-1]+m)%i; (i>1)

公式得到的是以0~n-1标注的最终序号;下面我们就推导公式。


给出一个序列,从0~n-1编号。其中,k代表出列的序号的下一个,即k-1出列。

a 0, 1, …, k-1, k, k+1, …, n-1

那么,出列的序号是(m-1)%n,k=m%n。出列k-1后,序列变为

b 0, 1, …, k-2, k, k+1, …, n-1

然后,我们继续从n-1后延长这个序列,可以得到

c` 0, 1, …, k-2, k, k+1, …, n-1, n, n+1, …, n+k-2

我们取从k开始直到n+k-2这段序列。其实这段序列可以看作将序列b的0~k-2段移到了b序列的后面。这样,得到一个新的序列

c k, k+1, …, n-1, n, n+1, …, n+k-2

好了,整个序列c都减除一个k,得到

0, 1, …, n-2

c序列中的n-1, n, n+1都减除个k是什么?这个不需要关心,反正c序列是连续的,我们知道了头和尾,就能知道d序列是什么样的。

这样你看,从序列a到序列d,就是一个n序列到n-1序列的变化,约瑟夫环可以通过递推来获得最终结果。ok,继续向下。

剩下的就是根据n-1序列递推到n序列。假设在n-1序列中,也就是序列d中,我们知道了最终剩下的一个序号是x,那么如果知道了x转换到序列a中的编号x`,不就是知道了最终的结果了么?

下面我们就开始推导出序列a中x的序号是什么。

d->c,这个变换很容易,就是x+k;

c->b,这个变换是网上大家都一带而过的,也是令我郁闷的一个关键点。从b->c,其实就是0~k-2这段序列转换为n~n+k-2这段序列,那么再翻转回去,简单的就是%n,即(x+k)%n。%n以后,k~n-1这段序列值不会发生变化,而n~n+k-2这段序列则变成了0~k-2;这两段序列合起来,就是序列b。

于是乎,我们就知道了,x`=(x+k)%n。并且,k=m%n,所以x`=(x+m%n)%n=(x+m)%n。公式1就出来了:f[i]=(f[i-1]+m)%i。当然,i=1就是特殊情况了,f[1]=0。这里还有一个小问题。也许你会迷惑为什么x`=(x+m%n)%n=(x+m)%n中的%n变成公式中f[i]=(f[i-1]+m)%i中的%i?其实这个稍微想想就能明了。我们%n就是为了从序列c转换到序列b——这是在n-1序列转换成n序列时%n;那么从n-2转换到n-1呢?不是要%(n-1)了吗?所以这个值是变量,不是常量。

好了,这个最后需要注意的就是从一开始,我们将n序列从0~n-1编号,所以依据公式得出的序号是基于0开始的。下面给出代码:

    #include<iostream>  
    #include<cstdio>  
    using namespace std;  
      
    int main()  
    {  
        int ncase, n, m, ans;  
        scanf("%d", &ncase);  
        while(ncase--)  
        {  
            scanf("%d%d", &n, &m);  
            ans = 0;  
            for(int i = 2; i <= n; ++i)  
            {  
                ans = ( ans + m ) % i;  
            }  
            printf("%d\n", ans + 1); //易错点,编号从1开始,所以结果+1
        }  
        return 0;  
    }  

循环单链表解法:

#include<stdio.h>   
#include <stdlib.h>  

typedef struct node//节点存放一个数据和指向下一个节点的指针  
{  
    int data;  
    struct node* pnext;  
} Node;  

Node *link_create(int n)//创建n个节点的循环链表  
{  
    //先创建第1个节点  
    Node *p,*q,*head;  
    int i ;  
    p = (Node *)malloc(sizeof(Node));  
    head = p;  
    p->data = 1;  

    for(i = 2;i<=n;i++)//再创建剩余节点  
    {  
        q = (Node *)malloc(sizeof(Node));  
        q->data = i;  
        p->pnext = q;  
        p = q;  
    }  
    p->pnext = head;//最后一个节点指向头部,形成循环链表  
    return head;  
}  

void link_process(Node *head,int k,int m)//从编号为k(1<=k<=n)的人开始报数,数到m的那个人出列;  
{  
    int i;  
    Node *p = head,*tmp1;  
    while(p->data != k)//p指向的节点data为k 
    {
        tmp1 = p;
        p = p->pnext;  
    }

    while(p->pnext != p)  
    {  
        for(i=1;i<m;i++)  
        {  
            tmp1 = p;  
            p = p->pnext;  
        }  
        printf("%d ",p->data);  
        tmp1->pnext= p->pnext;  
        free(p);  
        p = tmp1->pnext;  

    }  
    printf("%d ",p->data);  
    free(p);  
}  

int main()  
{  
    Node *head = link_create(5);      
    link_process(head,3,1);  
    system("pause");  
    return 0;  
}  




  • 4
    点赞
  • 7
    收藏
    觉得还不错? 一键收藏
  • 1
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值