[SGU 108]Self-numbers II(筛法活用优化+滚动数组)

【题目大意】：

定义一个数为self-number x，当且仅当找不到一个n使得x=d(n)。定义d(n)为n+(n的各位数之和)。给出一个n，求小于n的self-number个数，并且给出k个询问，每个询问为这个范围内第Sk个self-number。

【题目分析】：

首先要想一想这个题的特性。

self-number是一串非self-number派生的根源，所以我们想到了素数的筛法。但是我们注意到这样一点，那就是这道题的空间限制卡得太严了~

10^7的布尔数组完全无法承受，所以我们就要考虑如何精简空间的问题。

我们发现，两个相派生的数之间的差距不会超过63(7*9)。所以可以使用滚动数组~数组只要开到64就足够了……

对于询问进行排序，然后在筛出self-number的过程中进行回答。

还有一个小优化，就是说一个数可以被断成两截，预处理出0~9999的各位数之和，这样得到d(n)只需要O(1)的复杂度了~

 

P.s. 最开始犯了一个极其缺心眼的错误，因为很晚了，神智不太清醒，快排的递归部分没写~结果WA on 6(这种情况居然~~还能过5个点…………)

 

【代码】：

program SGU_108; var ques,id,ans:array[1..5000] of longint; v:array[0..63] of boolean; sum:array[0..10000] of longint; i,m,n,tmp,now,now1,tot,now2:longint; procedure qsort(s,t:longint); var i,j,key,k:longint; begin i:=s; j:=t; key:=ques[(s+t) shr 1]; repeat while ques[i]<key do inc(i); while ques[j]>key do dec(j); if i<=j then begin k:=ques[i]; ques[i]:=ques[j]; ques[j]:=k; k:=id[i]; id[i]:=id[j]; id[j]:=k; inc(i); dec(j); end; until i>j; if s<j then qsort(s,j); //A tragedy~ WA on 6 if i<t then qsort(i,t); end; begin assign(input,'a.in'); assign(output,'a.out'); reset(input); rewrite(output); readln(n,m); for i:=1 to m do begin read(ques[i]); id[i]:=i; end; qsort(1,m); //make the questions up-leveled for i:=0 to 10000 do begin sum[i]:=0; tmp:=i; while tmp>0 do //a very good predo begin inc(sum[i],tmp mod 10); tmp:=tmp div 10; end; end; now:=1; now2:=1; tot:=0; for i:=1 to n do begin if not v[now] then begin inc(tot); while (ques[now2]=tot) and (now2<=m) do //Caution, there may be same questions begin //or WA on 5 ans[id[now2]]:=i; inc(now2); end; end; tmp:=(sum[i div 10000]+sum[i mod 10000]+i) mod 64; //here can be O(1) v[tmp]:=true; v[now]:=false; //Caution~ This is a rolling array... This is very important now:=(now+1) mod 64; end; writeln(tot); write(ans[1]); for i:=2 to m do write(' ',ans[i]); writeln; close(input); close(output); end.