poj 2886 线段树 （用反素数）

题意：有n 个人，每个人有一个卡，上面有一个非零的数。第一次给一个数k，第k个人跳出，然后看第k的卡上数字是多少，如果大于0 则顺时针数那人人跳出。小于零，逆时针数到的那个人跳出，一次类推，第x次跳出的人会得到f（x），f（x）等于x的因子个数，求最大的f（x）

f（x）满足反素数。 线段树的区间存的是有多少人。

#include<cstdio>
#include<cstring>
#include<iostream>
using namespace std;
#define FF freopen("Input.txt","r",stdin)
const int N=5e5+2;
int RPrime[]=
{//反素数
1,2,4,6,12,24,36,48,60,120,180,240,360,720,840,1260,1680,2520,5040,7560,10080,15120,
20160,25200,27720,45360,50400,55440,83160,110880,166320,221760,277200,332640,498960,
554400
};
int fact[]=
{//反素数约数个数
1,2,3,4,6,8,9,10,12,16,18,20,24,30,32,36,40,48,60,64,72,80,84,90,96,100,108,120,128,
144,160,168,180,192,200,216
};
struct Node
{
char name[20];
int val;
}A
,res;
struct TREE
{
int l,r;
int sum;
int mid()
{
return (l+r)>>1;
}
}Tree[N*4];
void Build(int rt,int l,int r)
{
Tree[rt].l=l;
Tree[rt].r=r;
Tree[rt].sum=r-l+1;
if(l==r) return ;
int mid=Tree[rt].mid();
Build(rt<<1,l,mid);
Build(rt<<1|1,mid+1,r);
}
int Query(int rt,int x)
{
Tree[rt].sum--;
if(Tree[rt].l==Tree[rt].r) return Tree[rt].l;
if(x<=Tree[rt<<1].sum) return Query(rt<<1,x);
else return Query(rt<<1|1,x-Tree[rt<<1].sum);
}
int main()
{
//FF;
int n,k,i;
while(~scanf("%d%d",&n,&k))
{
for(i=1;i<=n;i++)
scanf("%s%d",A[i].name,&A[i].val);
int p,Max,id;
for(i=0;RPrime[i]<=n;i++);
p=RPrime[i-1];
Max=fact[i-1];
Build(1,1,n);
for(i=1;i<=p;i++)
{
id=Query(1,k);//删掉第k个.得到第k的位置id.
if(i==p) break;
n--;
if(A[id].val>0) k=(k-1+A[id].val-1)%n+1;
else k=((k-1+A[id].val)%n+n)%n+1;
}
printf("%s %d\n",A[id].name,Max);
}
return 0;
}

如果某个正整数x满足:对于任意i(0<i<x),都有g(i)<g(x),则称x为反素数. g(x) 表示x的约数的个数

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
const int N=500001;
const int prime[13]={0,2,3,5,7,11,13,17,19,23,29,31,37};
int ct[12],visit
;
int Rprime[100],fact[100];
void dfs(int x,int num,int sum)
{
visit[num]=sum;
for(int i=1,tmp=1;i<=min(19,ct[x-1]);i++)
{
tmp*=prime[x];
ct[x]=i;
if(num*tmp<N)
dfs(x+1, num*tmp, sum*(i+1));
else  break;
}
}
int main()
{
int i, mm = -1;
memset(visit,-1,sizeof(visit));
ct[0] = 19;
dfs(1,1,1);
int tmp=0;
for(i=1;i<N;i++)
{
if(visit[i]!=-1)
{
if(mm>=visit[i]) visit[i]=-1;
else
{
Rprime[tmp]=i;
fact[tmp++]=visit[i];
mm=visit[i];
}
}
}
return 0;
}