数论中的Josehus环有关问题
数论中的Josehus环问题
这个问题吧,其实从我一开始接触到算法的时候,就开始纠结了。但是至今也没真正懂得其本质含义(冏)。所以,只能转一些别人的博客,记录一下,慢慢研究。不过坑爹的是居然很多人的博客的程序是有BUG的,而且是一个BUG的程序被一堆人转载着。其实,我也是其中转载BUG程序的一个。后来,做题的时候才发现了漏洞的。现在就给出一个正确的程序和解释吧。
下面利用数学推导,如果能得出一个通式,就可以利用递归、循环等手段解决。下面给出推导的过程:
(1)第一个被删除的数为 (m - 1) % n。
(2)假设第二轮的开始数字为k,那么这n - 1个数构成的约瑟夫环为k, k + 1, k + 2, k +3, .....,k - 3, k - 2。做一个简单的映射。
k -----> 0
k+1 ------> 1
k+2 ------> 2
...
...
k-2 ------> n-2
这是一个n -1个人的问题,如果能从n - 1个人问题的解推出 n 个人问题的解,从而得到一个递推公式,那么问题就解决了。假如我们已经知道了n -1个人时,最后胜利者的编号为x,利用映射关系逆推,就可以得出n个人时,胜利者的编号为 (x + k) % n。其中k等于m % n。代入(x + k) % n <=> (x + (m % n))%n <=> (x%n + (m%n)%n)%n <=> (x%n+m%n)%n <=> (x+m)%n
(3)第二个被删除的数为(m - 1) % (n - 1)。
(4)假设第三轮的开始数字为o,那么这n - 2个数构成的约瑟夫环为o, o + 1, o + 2,......o - 3, o - 2.。继续做映射。
o -----> 0
o+1 ------> 1
o+2 ------> 2
...
...
o-2 ------> n-3
这是一个n - 2个人的问题。假设最后的胜利者为y,那么n -1个人时,胜利者为 (y + o) % (n -1 ),其中o等于m % (n -1 )。代入可得 (y+m) % (n-1)
要得到n - 1个人问题的解,只需得到n - 2个人问题的解,倒推下去。只有一个人时,胜利者就是编号0。下面给出递推式:
f [1] = 0;
f [ i ] = ( f [i -1] + m) % i; (i>1)
int JosephusProblem_Solution4(int n, int m) { if(n < 1 || m < 1) return -1; vector<int> f(n+1,0); for(unsigned i = 2; i <= n; i++) f[i] = (f[i-1] + m) % i; return f[n]; }
还有一个高深的链表法:
struct ListNode { int num; //编号 ListNode *next; //下一个 ListNode(int n = 0, ListNode *p = NULL) { num = n; next = p;} }; //自定义链表实现 int JosephusProblem_Solution1(int n, int m) { if(n < 1 || m < 1) return -1; ListNode *pHead = new ListNode(); //头结点 ListNode *pCurrentNode = pHead; //当前结点 ListNode *pLastNode = NULL; //前一个结点 unsigned i; //构造环链表 for(i = 1; i < n; i++) { pCurrentNode->next = new ListNode(i); pCurrentNode = pCurrentNode->next; } pCurrentNode->next = pHead; //循环遍历 pLastNode = pCurrentNode; pCurrentNode = pHead; while(pCurrentNode->next != pCurrentNode) { //前进m - 1步 for(i = 0; i < m-1; i++) { pLastNode = pCurrentNode; pCurrentNode = pCurrentNode->next; } //删除报到m - 1的数 pLastNode->next = pCurrentNode->next; delete pCurrentNode; pCurrentNode = pLastNode->next; } //释放空间 int result = pCurrentNode->num; delete pCurrentNode; return result; }
标准库法:
//使用标准库 int JosephusProblem_Solution2(int n, int m) { if(n < 1 || m < 1) return -1; list<int> listInt; unsigned i; //初始化链表 for(i = 0; i < n; i++) listInt.push_back(i); list<int>::iterator iterCurrent = listInt.begin(); while(listInt.size() > 1) { //前进m - 1步 for(i = 0; i < m-1; i++) { if(++iterCurrent == listInt.end()) iterCurrent = listInt.begin(); } //临时保存删除的结点 list<int>::iterator iterDel = iterCurrent; if(++iterCurrent == listInt.end()) iterCurrent = listInt.begin(); //删除结点 listInt.erase(iterDel); } return *iterCurrent; }
再次声明,此文代码和解释非原创。
给出一题说是模板题的题吧,纠结了好久,最后才知道。
题目链接:Click Here~
#include <iostream> #include <cstdio> #include <cstring> using namespace std; int Josehus(int n,int m) { int s = 0; //第一个人的编号 for(int i = 2;i <= n;++i) //遍历的是总人数,与编号没有任何关系!!! s = (s+m)%i; return s; //返回最后一个人的编号 } int main() { int n; while(scanf("%d",&n),n) { int m = 2; /*假设编号为1的已经出列,则此时从2开始重新编号为1。 所以,此时总人数应该是N-1.而且,我们知道我们在求解的时候 是假设第一个人的编号为0的。所以,现在如果要2是最后出列的则 需要在Josehus的函数中返回第一个值,即0*/ while(Josehus(n-1,m)!=0)m++; printf("%d\n",m); } return 0; }