51nod1989 竞赛表格
1989 竞赛表格
定义rev(i)为i十进制下各位翻转所得的数，例如rev(2345)=5432，rev(3210)=123。
l0nl1f3比较无聊，他打算用rev函数进行一个游戏。
他找了一个无穷大的表格，设第i行第j列的格子上的数为f(i,j)，那么ff满足
l0nl1f3想要知道这个表格中[L,R]之间的数出现了多少次，mod P输出。
由于他有时候会改变主意，所以他可能会询问q次，每次的L、R、P可能不同。
由于输入较慢，建议使用输入优化，这里提供一份C++的输入优化（标程就使用这个输入优化）：
char ch,B[1<<20],*S=B,*T=B;
#define getc() (S==T&&(T=(S=B)+fread(B,1,1<<20,stdin),S==T)?0:*S++)
#define isd(c) (c>='0'&&c<='9')
ll aa;ll F(){ //positive only.
while(ch=getc(),!isd(ch));aa=ch-'0';
while(ch=getc(),isd(ch))aa=aa*10+ch-'0';
return aa;
}由于输出较慢，你只需要将每个询问的答案异或后输出。
样例解释：
第一个询问，数字表格前13行、前4列如下：
1 2 4 8
2 4 8 16
3 6 12 33
4 8 16 77
5 10 11 22
6 12 33 66
7 14 55 110
8 16 77 154
9 18 99 198
10 11 22 44
11 22 44 88
12 33 66 132
13 44 88 176
可以发现其中有9个10~13的数，9 mod 7=2，所以该询问答案为2。
四个询问的答案分别为：2 159 147 1825，2 xor 159 xor 147 xor 1825=1839，所以输出1839。
输入
第一行一个整数q，表示询问次数。（ 1≤q≤10^5）
第二行开始q行每行三个整数L、R、P。（1≤L≤R≤10^10，1≤P≤10^9）输出
一行一个整数，为每个询问的答案的异或和。输入样例
4
10 13 7
2333 233333 233
6666 666666 666
1234 567890 2333输出样例
1839解析：
例如一个数x=abcde，rev(x)为edcba，那么x+rev(x)=10001(a+e)+1010(b+d)+200c。
那么可以发现a+e、b+d都是[0,18]的一个整数，c是[0,9]的整数。
那么我们可以考虑枚举x的位数和a+e、b+d、c这样的数，dfs一波，进而求出x+rev(x)可能取到的每个数，同时计算取到每个x+rev(x)的x的个数，全部记下来sort就好。
这样的复杂度大约是19^(len/2)*len的，可以求出250w左右个可能的数，记这个集合为S。
记f[i]为i在表格中的出现次数，g[i]为j+rev(j)=i的j个数，那么 
，所以 
。
显然只有g[i]>0的数f[i]才大于1，那么我们只要从小到大枚举S中的每个数p，每次让 
，
即可。
求出S集合中所有f值之后只要前缀和一下，每次询问求和一下，再加上不在集合S中的数的个数即可。
放代码：
#include
#define ll long long
using namespace std;
char ch,B[1<<20],*S=B,*T=B;
#define getc() (S==T&&(T=(S=B)+fread(B,1,1<<20,stdin),S==T)?0:*S++)
#define isd(c) (c>='0'&&c<='9')
ll aa;ll F(){
while(ch=getc(),!isd(ch));aa=ch-'0';
while(ch=getc(),isd(ch))aa=aa*10+ch-'0';
return aa;
}
ll psb[3200000],sb[3200000],pw[20],psum[3200000],sum[3200000];
int c[3200000];
void dfs(int now,int lim,ll had,ll ways){
if(had > pw[10]) return;
if(now == (lim+1)/2){
psb[++psb[0]] = had , c[psb[0]] = psb[0] , psum[psb[0]] = ways;
return;
}
if(!now){
if(pw[now] == pw[lim-now-1]){
for(int i=1;i<=9;i++)
dfs(now+1,lim,had+2*i*pw[now],ways);
}
else{
for(int i=1;i<=18;i++)
dfs(now+1,lim,had+i*(pw[now]+pw[lim-1-now]),ways*(i<=9?i:9-(i-10)));
}
}
else{
if(pw[now] == pw[lim-now-1]){
for(int i=0;i<=9;i++)
dfs(now+1,lim,had+2*i*pw[now],ways);
}
else{
for(int i=0;i<=18;i++)
dfs(now+1,lim,had+i*(pw[now]+pw[lim-1-now]),ways*(i<=9?i+1:9-(i-10)));
}
}
}
ll rev(ll x){
ll ret=0;
for(;x;) ret=ret*10+x%10,x/=10;
return ret;
}
bool cmp(const int &u,const int &v){ return psb[u] < psb[v]; }
int main(){
pw[0] = 1;
for(int i=1;i<=10;i++) pw[i] = pw[i-1] * 10;
for(int i=1;i<=10;i++) dfs(0,i,0,1);
sort(c+1,c+1+psb[0],cmp);
for(int i=1;i<=psb[0]+1;i++)
if(psb[c[i]] == psb[c[i-1]])
psum[c[i]] += psum[c[i-1]];
else if(i>1)
sb[++sb[0]] = psb[c[i-1]] , sum[sb[0]] = psum[c[i-1]];
for(int i=1;i<=sb[0];i++){
sum[lower_bound(sb+1,sb+1+sb[0],sb[i]+rev(sb[i]))-sb] += sum[i];
sum[i] += sum[i-1];
}
int q=F(),ans=0;
ll l,r,p;
for(;q--;){
l=F(),r=F(),p=F();
int ret = (sum[upper_bound(sb+1,sb+1+sb[0],r)-sb-1] - sum[upper_bound(sb+1,sb+1+sb[0],l-1)-sb-1] + r - l + 1) % p;
ans ^= ret;
}
printf("%d\n",ans);
return 0;
}