week8实验C 咕咕东的奇妙序列

题目：
咕咕东 正在上可怕的复变函数，但对于稳拿A Plus的 咕咕东来说，她早已不再听课，此时她在睡梦中 突然想到了一个奇怪的无限序列：112123123412345 …这个序列由连续正整数组成的若干部分构成，其中第一部分包含1至1之间的所有数字，第二部分包含1至2之间的所有数字，第三部分包含1至3之间的所有数字，第i部分总是包含1至i之间的所有数字。

所以，这个序列的前56项会是11212312341234512345612345671234567812345678912345678910，其中第1项是1，第3项是2，第20项是 5，第38项是2，第56项是0。咕咕东 现在想知道第 k 项数字是多少！但是她睡醒之后发现老师讲的东西 已经听不懂了，因此她把这个任务交给了你。
输入：

输出：

样例输入：
5
1
3
20 38 56
样例输出：
1
2
5
2
0
这个序列要注意11212312341234512345612345671234567812345678912345678910是56项，而不是55项，10被分成了1和0两项。
首先，计算出i位数的总序列长度，就好像对于2位数，计算出112123…12345678910开始，到112123…1234…9899截止，这其中所有的序列的长度之和。
令a[i]记录i位数的总序列长度，b[i]记录i位数最长的序列长度，例如112123…1234…99的长度是b[2]的值，可以用等差数列的求和公式计算。

long long a[30],b[30];
long long fun1()
{
 a[0]=b[0]=0;
 long long count=9;
 long long i=1;
 while(1)
 {
  a[i]=(b[i-1]+i + b[i-1]+i*count)*count/2;
  b[i]=b[i-1]+i*count;
  count*=10;
  if(a[i]>(long long)pow(10,18))//计算到1e18为止
  {
   break;
  }
  i++;
 }
 return i;
}

然后判断输入的k是几位数：

 long long x;
 for(long long i=1;i<=max_len;i++)//确定k是几位数 
 {
  if(a[i]>=k)
  {
   x=i;
   break;
  }
  k=k-a[i];
 }

通过二分查找，找到k是x位数序列的第几个，然后减去前面的序列长度。

 //二分查找k是x位数序列的第几个 
 long long l=1,r=(long long)(9*pow(10,x-1));
 while(l<=r)
 {
  long long mid=(l+r)/2;
  long long temp=(b[x-1]+x + b[x-1]+x*mid)*mid/2;
  if(temp>=k) 
   r=mid-1;
  else 
   l=mid+1;
 }
 k=k-(b[x-1]+x + b[x-1]+x*r)*r/2;//去掉序列的前r个序列的长度

然后再确定是第r个序列中的第几部分

 long long y=1;
 while(1)
 {
  long long sum=9*pow(10,y-1)*y;//从序列首部开始，每次去掉序列的1~9，10~99以此类推 
  if(sum>=k)
  {
   break; 
  } 
  k=k-sum;
  y++;
 }
 long long cnt=(k+y-1)/y;//找出k是y位数序列中的第几个 
 k=k-(cnt-1)*y;//去掉前cnt-1个数

最后复原该数得到k所对应的那一位。

 long long number=pow(10,y-1)+cnt-1;//复原该数 
 long long t=y-k+1;//第y个数的第t位
 for(long long i=1;i<t;i++)
 {
  number/=10;//去掉后t-1位 
 }
 //取个位
 printf("%lld\n",number%10);

以下是完整代码：

#include<iostream>
#include<cmath>
using namespace std;
long long a[30],b[30];
long long fun1()
{
 a[0]=b[0]=0;
 long long count=9;
 long long i=1;
 while(1)
 {
  a[i]=(b[i-1]+i + b[i-1]+i*count)*count/2;
  b[i]=b[i-1]+i*count;
  count*=10;
  if(a[i]>(long long)pow(10,18))
  {
   break;
  }
  i++;
 }
 return i;
}
void fun2(long long max_len,long long k)
{
 long long x;
 for(long long i=1;i<=max_len;i++)//确定k是几位数 
 {
  if(a[i]>=k)
  {
   x=i;
   break;
  }
  k=k-a[i];
 }
 //二分查找k是x位数序列的第几个 
 long long l=1,r=(long long)(9*pow(10,x-1));
 while(l<=r)
 {
  long long mid=(l+r)/2;
  long long temp=(b[x-1]+x + b[x-1]+x*mid)*mid/2;
  if(temp>=k) 
   r=mid-1;
  else 
   l=mid+1;
 }
 k=k-(b[x-1]+x + b[x-1]+x*r)*r/2;//去掉序列的前r个
 long long y=1;
 while(1)
 {
  long long sum=9*pow(10,y-1)*y;//从序列首部开始，每次去掉序列的1~9，10~99以此类推 
  if(sum>=k)
  {
   break; 
  } 
  k=k-sum;
  y++;
 }
 long long cnt=(k+y-1)/y;//找出k是y位数序列中的第几个 
 k=k-(cnt-1)*y;//去掉前cnt-1个数
 long long number=pow(10,y-1)+cnt-1;//复原该数 
 long long t=y-k+1;//第y个数的第t位
 for(long long i=1;i<t;i++)
 {
  number/=10;//去掉后t-1位 
 }
 //取个位
 printf("%lld\n",number%10);
}
int main()
{
 long long len=fun1();
 int q;
 cin>>q;
 for(int i=0;i<q;i++)
 {
  long long kk;
  cin>>kk;
  fun2(len,kk);
 }
}