首页 > 代码库 > POJ1019 NumberSequence 【数学公式转化题】

POJ1019 NumberSequence 【数学公式转化题】

这题也是一看感觉自己能做,就拿出笔和纸开始写写画画了,写的也蛮爽(下午1点到晚上8点,,,),其实,,,也就俩方程的事情。。。但是中间出了个逗比错误,10+11+12+...+99算成了(1+90)*90/2(找出这个错误真心不容易,现学了点linuxshell知识,也正好写了写相关的脚本,比如对拍脚本,运行脚本,挺好的)

这个错误也是奇葩,,前7000个竟然都没错,,真是 奇葩~

感觉数学题还是要做慢点,尤其是写完一个重要的东西,检验一下,,,错了就毙掉了,很花时间,在代码量还算少的时候就调试好

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <cmath>
using namespace std;
long long Fsave[99999];
long long fsave[99999];
char aint[59999];
inline long long  g(long long  n)//f函数的一部分,f函数是1234567891011...到每个数有几个数位f=g(n)+xn
{
    long long  tens=1;
    for(long long  i=0;i<n-1;i++)
    {
        tens*=10;
    }
    return (1-tens)/9-tens+n;
}
inline long long  getn(long long  x)//获得一个10进制数是几位数
{
    long long  co=0;
    while(x)
    {
        x/=10;
        co++;
    }
    return co;
}
inline long long  F(long long  x,long long  n)//1121231234最终序列的,到每个数而言有几个数位
{
    long long  sum=0;
    long long  tenPow=1;
    for(long long  i=1;i<=n-1;i++)
    {
        sum+=9*tenPow*g(i)+(tenPow+tenPow*10-1)*9*tenPow/2*i;
        tenPow*=10;
    }
    //cout<<"sum:"<<sum<<endl;
    n>=2?tenPow:tenPow=1;
    sum+=(x-tenPow+1)*g(n)+n*(tenPow+x)*(x-tenPow+1)/2;
    //cout<<tenPow<<endl;
    //cout<<n*(tenPow+x)*(x-tenPow+1)/2<<endl;
    return sum;
}
inline int binsearch(long long x,long long* A)从A中找出x的位置,夹在中间则取左边
{
    int L=0,R=104000;
    int mid=(L+R)>>1;
    while(L<R)
    {
        if(A[mid]>x)
            R=mid;
        if(A[mid]<x)
            L=mid+1;
        if(A[mid]==x)
            return mid;
        mid=(L+R)>>1;
    }
    return mid-1;
}
int main()
{
    //freopen("/home/test/in.txt","r",stdin);
    //freopen("/home/test/list.txt","w",stdout);
    for(long long  x=1;x<=52000;x++)//俩函数打表,用于binarysearch二分搜索
    {
        long long  n=getn(x);
        Fsave[x]=F(x,n);
    }
    for(long long  x=1;x<=52000;x++)
    {
        long long  n=getn(x);
        fsave[x]=g(n)+n*x;
    }
    long long t;scanf("%lld",&t);
    while(t--)
    {
        long long x;scanf("%lld",&x);
		    long long index_F=binsearch(x,Fsave);
		    x=x-Fsave[index_F];//x对应F序列中123...index_F的第几位
		    if(x==0)//如果恰好是结尾
		    {
		        printf("%lld\n",index_F%10);
		        continue;
		    }
		    long long index_f=binsearch(x,fsave);//x对应f中哪个数,f是123456789101112...
		    x=x-fsave[index_f];//x为该数字的第几位,比如122354第四位就是3
		    if(x==0)//刚好结尾
		    {
		        printf("%lld\n",index_f%10);
		        continue;
		    }
		    index_f++;
		    sprintf(aint,"%lld",index_f);
		    printf("%c\n",aint[x-1]);
        
    }
	return 0;
}