POJ 2886 【线段树&&反素数&&约瑟夫循环】

POJ 2886 【线段树&&反素数&&约瑟夫循环】
题目链接
题意：给你n,k<=5e5,n是人数下面还会给你n个人的信息，k是初始位置，每次都会根据当前的位置大于0决定顺时针走a[i]，否则就走逆时针a[i],循环的次数就是小于n的因子最大的那个。题意是这样
思路：
1.首先要确定循环次数，暴力找5e5以内因子数最大的那个n*sqrt(n),加上多组肯定超时，明眼人一看就是反素数，看到很多博客写的都是直接给出5e5以内的反素数及因子个数，很草率，于是博主写了个反素数筛，挺容易理解的。
2.每次根据当前位置及信息能不能推出下一次位置呢？
分两种情况，顺时针的话就它的后面的下标就会-1，以为该位置被删除了，循环节肯定是res,现在公式就是（k+m-1)%res,但是考虑到%res可能会出现0,因此可以在里面-1,外边+1消除这种影响（k+m-2）%res+1;逆时针同理：但是它之前的下标不变公式就是（k+m-1）%res+1;
3.但是这个过程怎么实现呢？每次都会删除相当于一个动态的过程，看过POJ2828你肯定就知道了:插空法，每次删除一个人就等于那个空被占用了，pos就是第pos个空位，直接上线段树或者树状数组就ok。

#include<iostream>
#include<cstring>
#include<algorithm>
#include<cstdio>
using namespace std;
const int maxn=5e5+5;
const int p[]= {0,2,3,5,7,11,13,17,19,23,29,31};
int ma=0,ans=0,n,k,pos;
struct node
{
    int l,r,num;
} a[maxn<<2];
struct node1
{
    int num;
    char name[10];
} b[maxn];
void dfs(int now,int cnt,int step)
{
    if(cnt>ma)
    {
        ma=cnt;
        ans=now;
    }
    else if(cnt==ma)
        ans=min(ans,now);
    if(step>11) return ;
    for(int i=1; i<=31; i++)
    {
        if(now>n/p[step])
            return;
        dfs(now*p[step],cnt*(i+1),step+1);
        now=now*p[step];
    }
}
void pushup(int i)
{
    a[i].num=a[i*2].num+a[i*2+1].num;
}
void build(int i,int l,int r)
{
    a[i].l=l,a[i].r=r;
    a[i].num=r-l+1;
    if(l==r)
        return ;
    int mid=(l+r)/2;
    build(i*2,l,mid);
    build(i*2+1,mid+1,r);
}
int update(int i,int l,int r,int pos)
{
    if(l==r)
    {
        a[i].num=0;
        return l;
    }
    int mid=(l+r)/2,ans;
    if(pos<=a[i*2].num)
        ans=update(i*2,l,mid,pos);
    else
        ans=update(i*2+1,mid+1,r,pos-a[i*2].num);
    pushup(i);
    return ans;
}
int main()
{

    while(~scanf("%d%d",&n,&k))
    {
        build(1,1,n);
        dfs(1,1,1);
        for(int i=1; i<=n; i++)
            scanf("%s %d",&b[i].name,&b[i].num);
        int res=n;
        while(ans--)
        {
            res--;
            pos= update(1,1,n,k);
            if(!res)break;
            if(b[pos].num>0)
                k=(k-2+b[pos].num)%res+1;
            else
                k=((k-1+b[pos].num)%res+res)%res+1;
        }
        printf("%s %d\n",b[pos].name,ma);
    }
}