給出一個reverse_prime，自身是一個7位數，反轉後是一個<=1e6的素數。
首先求出所有的這種數
兩種操作
q k :表示刪除數字K
http://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&page=show_problem&problem=2657
首先預處理出所有的reverse_prime，我的做法是先打出<=1e6的所有素數
然後將其反轉，可以看出所有素數都是6位數，但是題目要求是7位數，可見原數的最低位都為0，可以不考慮這個0，之後的效率也許會快點。
將所有素數反轉，然後湊成6位數。再算出每個數的素因子個數，不要忘記之前少算的末尾的0，也就是因子2和因子5.
對於統計部分，建立兩個樹狀數組，第一 個表示區間內還有多少個數，第二個表示區間內因子個數和。
對於D操作， 直接用map記錄某個reverse_prime的下標，然後更新兩個樹狀數組
對於q操作，二分位置，然後用第一個樹狀數組，可以知道區間內有多少個數。最後用第二個樹狀數組求和。
[cpp]
#include<iostream>
#include<cstdio>
#include<map>
#include<cstring>
#include<cmath>
#include<vector>
#include<algorithm>
#include<set>
#include<string>
#include<queue>
#define inf 1600005
#define M 40
#define N 1000000
#define maxn 300005
#define eps 1e-12
#define zero(a) fabs(a)<eps
#define Min(a,b) ((a)<(b)?(a):(b))
#define Max(a,b) ((a)>(b)?(a):(b))
#define pb(a) push_back(a)
#define mp(a,b) make_pair(a,b)
#define mem(a,b) memset(a,b,sizeof(a))
#define LL unsigned long long
#define MOD 1000000007
#define lson step<<1
#define rson step<<1|1
#define sqr(a) ((a)*(a))
#define Key_value ch[ch[root][1]][0]
#define test puts("OK");
#define pi acos(-1.0)
#define lowbit(x) ((-(x))&(x))
#pragma comment(linker, "/STACK:1024000000,1024000000")
using namespace std;
int flag[N]={0},prime[N],cnt=0;
int fac[N],a[N],tot=0,p[N];
LL s1[N],s2[N];
map<int,int>m;
int slove(int num){
int len=0,ret=0,bit[20];
while(num){
bit[len++]=num%10;
num/=10;
}
for(int i=0;i<len;i++)
ret=ret*10+bit[i];
while(ret<100000) ret*=10;
return ret;
}
void Init(){
for(int i=2;i<N;i++){
if(flag[i]) continue;
prime[++cnt]=i;
for(int j=2;j*i<N;j++)
flag[i*j]=1;
}
for(int i=1;i<=cnt;i++){
a[i]=slove(prime[i]);
}
sort(a+1,a+1+cnt);
for(int i=1;i<=cnt;i++){
m[a[i]]=i;
}
for(int i=1;i<=cnt;i++){
fac[i]=2; //去掉了最後一個0，肯定包括2，5兩個因子
int tmp=a[i];
for(int j=1;j<=cnt&&prime[j]*prime[j]<=tmp;j++)
while(tmp%prime[j]==0){
tmp/=prime[j];
fac[i]++;
}
if(tmp>1) fac[i]++;
}
}
void Update(LL *s,int x,int val){
for(int i=x;i<=cnt;i+=lowbit(i))
s[i]+=val;
}
LL sum(LL *s,int x){
LL ret=0;
for(int i=x;i>0;i-=lowbit(i))
ret+=s[i];
return ret;
}
int main(){
//freopen("in.txt","r",stdin);
Init();
char str[5];int k;
mem(s1,0);mem(s2,0);
for(int i=1;i<=cnt;i++) s1[i]=lowbit(i);
for(int i=1;i<=cnt;i++) Update(s2,i,fac[i]);
while(scanf("%s%d",str,&k)!=EOF){
if(str[0]=='d'){
int pos=m[k/10];
Update(s1,pos,-1);
Update(s2,pos,-fac[pos]);
}
else{
k++; www.2cto.com
int low=1,high=cnt,mid;
while(low<=high){
mid=(low+high)>>1;
LL tmp=sum(s1,mid);
if(tmp==k) break;
if(tmp<k) low=mid+1;
else high=mid-1;
}
printf("%lld\n",sum(s2,mid));
}
}
return 0;
}
