约瑟夫问题(多解)——POJ 3750
对应POJ题目：点击打开链接
| Time Limit: 1000MS | Memory Limit: 65536KB | 64bit IO Format: %I64d & %I64u |
Description
Input
接下来每行输入一个小孩的名字(人名不超过15个字符)
最后一行输入W,S (W < N)，用逗号","间隔
Output
Sample Input
5 Xiaoming Xiaohua Xiaowang Zhangsan Lisi 2,3
Sample Output
Zhangsan Xiaohua Xiaoming Xiaowang Lisi
Source
思路：
1，循环链表，直接模拟。O(n*s)
2，线段树，O(nlg(n))
例子如下：n = 5, w = 1, s = 3 (括号内为原始编号)
1(1)，2(2)，3(3)，4(4)，5(5) 第3个数3出列
1(1)，2(2)，3(4)，4(5) 第1个数1出列
1(2)，2(4)，3(5) 第3个数5出列
1(2)，2(4) 第1个数2出列
1(4) 第1个数4出列
先引入Joseph递推公式，设有n个人（0,...,n-1），数到m退出，则第i轮出局的人为f(i)=(f(i-1)+m-1)%(n-i+1),f(0)=0; f(i) 表示当前子序列中要退出的那个人（当前序列编号为0~(n-i));
f(0) = 0;
第1轮：f(1) = (0 + 3 - 1) % 5 = 2 （下标为2即为第3个数）
第2轮：f(2) = (2 + 3 - 1) % 4 = 0
第3轮：f(3) = (0 + 3 - 1) % 3 = 2
第4轮：f(4) = (2 + 3 - 1) % 2 = 0
第5轮：f(5) = (0 + 3 - 1) % 1 = 0
所以使用这个公式每一轮中每个人的编号的，接着我们需要快速地找到这一轮中要退出的人的编号对应的原始编号。一个人在当前剩余队列中编号为i，则说明他是从左到右数第i个人，这启发我们可以用线段树来解决问题。用线段树维护原编号i..j内还有多少人没 有被淘汰，这样每次选出被淘汰者后，在当前线段树中查找位置就可以了。
例如我们有5个原编号，当前淘汰者在剩余队列中编号为3，先看左子树，即原编号1..3区间内，如果剩下的人不足3个，则说明当前剩余编号为3的 这个人原编号只能是在4..5区间内，继续在4..5上搜索；如果1..3内剩下的人大于等于3个，则说明就在1..3内，也继续缩小范围查找，这样既可 在logn时间内完成对应。问题得到圆满的解决。
3，网上借鉴的令我无法理解的方法：O(n)
http://www.cnblogs.com/void/archive/2011/04/21/2024377.html
方法1：
#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<map>
#include<queue>
#include<stack>
#include<vector>
#include<algorithm>
#include<cstring>
#include<string>
#include<iostream>
#define ms(x,y) memset(x,y,sizeof(x))
const int MAXN=1000+10;
const int INF=1<<30;
using namespace std;
struct Node
{
char name[20];
Node *prior;
Node *next;
};
void CreateList(Node *&head, int n)
{
Node *p,*s;
for(int i=0; i<n; i++){
s=(Node *)malloc(sizeof(Node));
scanf("%s", s->name);
if(head == NULL) head=s;
else{
p->next = s;
s->prior = p;
}
p=s;
}
p->next = head;
head->prior = p;
}
Node *Begin(Node *node, int s)
{
Node *p=node;
for(int i=0; i<s-1; i++)
p=p->next;
return p;
}
void Delete(Node *&node, int s, int n)
{
Node *q,*p=node;
for(int i=0; i<n; i++){
for(int i=0; i<s-1; i++)
p=p->next;
q=p->next;
printf("%s\n", p->name);
p->prior->next = p->next;
p->next->prior = p->prior;
free(p);
p=q;
}
}
int main()
{
//freopen("in.txt","r",stdin);
int n,w,s;
scanf("%d", &n);
Node *head=NULL;
CreateList(head, n);
scanf("%d,%d", &w,&s);
Node *p=Begin(head, w);
Delete(p, s, n);
return 0;
}
方法2
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#define N 70
#define L 20
char name[N][L];
int sum[N<<4];
void build(int root, int l, int r)
{
if(1 == r - l){
sum[root] = 1;
return;
}
int mid = ((l + r)>>1);
build(root<<1, l, mid);
build(root<<1|1, mid, r);
sum[root] = sum[root<<1] + sum[root<<1|1];
}
int query(int root, int l, int r, int num)
{
if(1 == r - l){
sum[root] = 0;
return l;
}
int mid = ((l + r)>>1);
int ans;
if(sum[root<<1] > num) ans = query(root<<1, l, mid, num);
else ans = query(root<<1|1, mid, r, num - sum[root<<1]);
sum[root] = sum[root<<1] + sum[root<<1|1];
return ans;
}
int main()
{
//freopen("in.txt", "r", stdin);
int n, w, s, i, out, num;
while(~scanf("%d", &n))
{
for(i=0; i<n; i++)
scanf("%s", name[i]);
scanf("%d,%d", &w, &s);
build(1, 0, n);
out = w - 1;
for(i=0; i<n; i++){
out = (out + s - 1) % (n - i); //第i+1轮退出的相对编号
num = query(1, 0, n, out); //根据相对编号找原始编号
//printf("%s\n", name[(num + w - 1) % n]);
printf("%s\n", name[num]);
}
}
return 0;
}
方法3
#include <stdio.h>
#define N 70
#define L 20
char name[N][L];
int main()
{
//freopen("in.txt", "r", stdin);
int n, i, m, p, w, s;
while(~scanf("%d", &n))
{
for(i=0; i<n; i++)
scanf("%s", name[i]);
scanf("%d,%d", &w, &s);
i = 0;
while( ++i <= n )
{
p = i * s;
while (p > n)
p = p - n + (p - n - 1)/(s - 1);
printf("%s\n",name[(p - 2 + w) % n]);
}
}
return 0;
}