poj2886 Who Gets the Most Candies?
这道题题意对于一般的ACMER应该不是什么问题吧！
首先是个约瑟夫问题的升级版
1）.告诉你第一次是顺时针第k个人出队
2）.某个人（编号为 i） 出队后，下个出队的是从这个开始查第num[i]个人，若num[i]<0,为逆时针数，否则顺时针数。
3）.重复步骤2）.，直到全部出队。
4）.输出出队序号中，那个序号的约数最多的那个人名及约数个数。
这里需要知道一个知识：反素数。
详见http://www.cnblogs.com/shen1000/archive/2012/07/29/2613877.html
有了反素数，我们就知道那个序号为最终答案了。
怎么求出那个序号是谁呢？对于这个，目前我只知道用模拟来找出。
思路:每次找出第几个需要要出队的人得位置，然后根据位置，找到下个需要出队的位置。
线段树功能:update:出队一个人，以及找到出队人的下标。 query:找出出对人的左边有多少未出队的人。
详见代码：
View Code
1 #include<iostream> 2 #include<string> 3 #include<queue> 4 #include<map> 5 #include<cmath> 6 #include<stack> 7 #include<algorithm> 8 using namespace std; 9 const int maxn = 500555; 10 char name[maxn][12]; 11 int num[maxn]; 12 int rprim[35][2] = { 13 498960,200,332640,192,277200,180,221760,168,166320,160, 14 110880,144,83160,128,55440,120,50400,108,45360,100, 15 27720,96,25200,90,20160,84,15120,80,10080,72, 16 7560,64,5040,60,2520,48,1680,40,1260,36, 17 840,32,720,30,360,24,240,20,180,18, 18 120,16,60,12,48,10,36,9,24,8, 19 12,6,6,4,4,3,2,2,1,1}; 20 21 //1.bulid(); 22 //2.query(a,b) 23 //3.update(a,b) 24 #define lson l , m , rt << 1 25 #define rson m + 1 , r , rt << 1 | 1 26 27 int sum[maxn<<2]; 28 int n,k; 29 //根据题意做相关修改，询问时的操作 30 int operate(int a,int b){ 31 return a+b; 32 } 33 void PushUp(int rt){ 34 sum[rt]=operate(sum[rt<<1],sum[rt<<1|1]); 35 } 36 void bulid(int l=1,int r=n,int rt=1){ 37 if(l==r){// 据题意做相关修改 38 sum[rt]=1;return ; 39 } 40 int m=(l+r)>>1; 41 bulid(lson); 42 bulid(rson); 43 PushUp(rt); 44 } 45 void update(int p,int add=0,int l=1,int r=n,int rt=1){ 46 if(l==r){// 据题意做相关修改 47 sum[rt]=add;k=l; 48 //printf("k=%d\n",k); 49 return ; 50 } 51 int m=(l+r)>>1; 52 if(p<=sum[rt<<1])update(p,add,lson); 53 else update(p-sum[rt<<1],add,rson); 54 PushUp(rt); 55 } 56 57 int query(int L,int R,int l=1,int r=n,int rt=1){ 58 if(L<=l && r<=R){ 59 return sum[rt]; 60 } 61 int m=(l+r)>>1; 62 int ret=0; 63 if(L<=m)ret=operate(ret,query(L,R,lson)); 64 if(R> m)ret=operate(ret,query(L,R,rson)); 65 return ret; 66 } 67 68 int main(){ 69 int p,x,now,leftnum,rightnum,m; 70 while(~scanf("%d%d",&n,&k)){ 71 p=0; 72 while(n<rprim[p][0])p++; 73 x=rprim[p][0]; 74 bulid(); 75 //printf("sum[1]=%d\n",sum[1]); 76 for(int i=1;i<=n;i++){ 77 scanf("%s%d",name[i],&num[i]); 78 } 79 m=n; 80 now=k; 81 //printf("x=%d\n",x); 82 for(int i=1;i<x;i++){ 83 update(now);//now出队 84 m--; 85 //printf("m=%d\n",m); 86 //转化为向左查多少人 87 if(num[k]%m==0){ 88 if(num[k]>0)num[k]=m; 89 else num[k]=1; 90 }else{ 91 num[k]%=m; 92 if(num[k]<0)num[k]+=m+1; 93 } 94 leftnum=query(1,k);//找出now左边有多少个人 95 rightnum=m-leftnum; 96 //printf("num[%d]=%d left=%d right=%d\n",k,num[k],leftnum,rightnum); 97 if(num[k]<=rightnum){ 98 now=leftnum+num[k]; 99 }else{ 100 now=num[k]-rightnum; 101 } 102 } 103 update(now); 104 printf("%s %d\n",name[k],rprim[p][1]); 105 106 } 107 108 return 0; 109 } 110 /* 111 4 1 112 Tom 2 113 Jack 4 114 Mary -1 115 Sam 1 116 */