poj 2886 Who Gets the Most Candies?

写代码的时间长达3个半小时，当然没刨去中间吃饭什么的........

思路：先直接dfs求出最大的 p 和 f(p) 这里用到了反素数，详解请看zoj 2562那篇文章，里面有解释

然后就是用线段树模拟约瑟夫，但是悲剧的是我没太搞懂出圈顺序，此处系参考而来....................

code

#include <iostream>
#include <fstream>
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <string>
#include <string.h>
#include <vector>
#include <bitset>
#include <cmath>
#include <queue>
#include <stack>
#include <set>
#include <ctime>
#include <map>
#include <limits>
#define LL long long
#define Vi vector<int>
#define Si set<int>
#define readf freopen("input.txt","r",stdin)
#define writef freopen("output.txt","w",stdout)
#define FF(i,a) for(int i(0); i < (a); i++)
#define FD(i,a) for(int i(a); i >= (1); i--)
#define FOR(i,a,b) for(int i(a);i <= (b); i++)
#define FOD(i,a,b) for(int i(a);i >= (b); i--)
#define PD(a) printf("%d",a)
#define SET(a,b) memset(a,b,sizeof(a))
#define SD(a) scanf("%d",&(a))
#define LN printf("\n")
#define PS printf(" ")
#define pb push_back
#define lson l , m , rt << 1
#define rson m + 1 , r , rt << 1 | 1
const double pi = acos(-1.0);
const int maxn = 500001;
const int INF = 99999999;
const int dx[]={0,1,0,-1};
const int dy[]={1,0,-1,0};
using namespace std;
struct stu{
    char name[20];
    int val;
}p[maxn];
int prime[7]={2,3,5,7,11,13,17};
int N,K;
int ans,mcnt;   //ans为最大的反质数，mcnt为最大反质数的因子数
void dfs(int num,int cnt,int k,int lim){
    if(k>7) return;
    if(cnt>mcnt){
        mcnt=cnt;
        ans=num;
    }
    if(cnt==mcnt && num<ans) ans=num;
    int tmp=num;
    FOR(i,1,lim){
        if(tmp*prime[k]>N) return ;
        tmp*=prime[k];
        dfs(tmp,cnt*(i+1),k+1,i);
    }
}
int len[maxn<<2];
void PushUP(int rt){
    len[rt]=len[rt<<1]+len[rt<<1|1];
}
void build(int l,int r,int rt){
    if(l==r){
        len[rt]=1;
        return ;
    }
    int m=(l+r)>>1;
    build(lson);
    build(rson);
    PushUP(rt);
}
int query(int p,int l,int r,int rt){
    if(l==r){
        len[rt]=0;
        return l;
    }
    int m=(l+r)>>1;
    int ret;
    if(p<=len[rt<<1]) ret=query(p,lson);
    else ret=query(p-len[rt<<1],rson);
    PushUP(rt);
    return ret;
}
int main()
{
    int &mod=len[1];
    while(~scanf("%d%d",&N,&K)){
        ans=N;mcnt=0;
        dfs(1,1,0,20);
        FOR(i,1,N) scanf("%s%d",p[i].name,&p[i].val);
        build(1,N,1);
        int pos;
        FOR(i,0,ans-1)
        {
            if(p[pos].val>0) K=((K+p[pos].val-2)%mod+mod)%mod+1;   //
            else    K=((K+p[pos].val-1)%mod+mod)%mod+1;         //说实话这个真心不懂
            pos=query(K,1,N,1);
        }
        printf("%s %d\n",p[pos].name,mcnt);
    }
    return 0;
}