【模拟】Vijos P1005 超长数字串

题目链接：

　　https://vijos.org/p/1005

题目大意：

　　无限的正整数按顺序拼接成字符串S(S=12345678910111213...)，给你一个字符串A(len<=200)求这个字符串在S中最早出现的位置。

　　(答案超过long long ，无法用KMP，不要相信标签)

题目思路：

　　【模拟】

　　这题简直了！！！！！！大模拟啊。细节超级多。疯狂TLE+WA+RE了17次才AC。

　　第一次写了一整天没写过，放了好久，昨天又写了一整天(哎效率低下。)

　　首先分两种普遍的情况和两种特殊情况考虑。

　　普遍情况一是S串是由三段组成，头，尾和中间数字串，头尾可能不完全。这时候枚举中间数字串的长度和起始位置，就可以判断出最早出现的数字

　　(如：123[98 12399 124]00  99[8 999 1000 10]01 等)

　　(一开始我是想把这个再分情况考虑，发现有进位什么的太麻烦了后来直接把这个数压到高精度中，先-1判断和头是否匹配，再每次+1判断和后面的字符串匹不匹配，尾部同理)

　　普遍情况二是S串是由两段组成，前面的数字的尾部和后面的数字的头部组成，这时要枚举重叠部分的长度和断开的位置。

　　(如：1[2 1]3  123[699 1237]00    60211602应该是 1[6021 1602]2 而不是[60211 602]12)

　　(一开始我忘记了这种情况WA了好久。。感觉自己好蠢。)

　　特殊情况一是全为0的情况，这是都不符合上面的情况，就需要在头部再添上一个1，肯定是最早出现的。

　　特殊情况二是全为9的情况，除了单个9是9以外，两位数以上的9都可以拆成 8[9999... 9]0000... 这个答案明显比[9999...] 10000...优。

　　（不知道自己还有没有漏情况，希望没有吧，望指教。）

　　

 //
 //by coolxxx
 //#include<bits/stdc++.h>
 #include<iostream>
 #include<algorithm>
 #include<string>
 #include<iomanip>
 #include<map>
 #include<memory.h>
 #include<time.h>
 #include<stdio.h>
 #include<stdlib.h>
 #include<string.h>
 //#include<stdbool.h>
 #include<math.h>
 #define min(a,b) ((a)<(b)?(a):(b))
 #define max(a,b) ((a)>(b)?(a):(b))
 #define abs(a) ((a)>0?(a):(-(a)))
 #define lowbit(a) (a&(-a))
 #define sqr(a) ((a)*(a))
 #define swap(a,b) ((a)^=(b),(b)^=(a),(a)^=(b))
 #define mem(a,b) memset(a,b,sizeof(a))
 #define eps (1e-8)
 #define J 10
 #define mod 1000000007
 #define MAX 0x7f7f7f7f
 #define PI 3.14159265358979323
 #define N 404
 using namespace std;
 typedef long long LL;
 int cas,cass;
 int n,m,lll,ans;
 double anss;
 char s[N];
 int a[N],b[N],t[N],c[N];
 void gjdchange(int a[],int l,int r)//l~r中间的字符转化为高精度数a
 {
     int i;
     mem(t,);
     t[]=r-l+;
     for(i=;i<=r-l+;i++)
         t[i]=s[r-i+]-'';
     memcpy(a,t,sizeof(t));
 }
 void plu(int a[])//高精度a+1
 {
     int i;
     for(i=,a[]++;a[i]>=J;i++)
         a[i]-=J,a[i+]++;
     a[]=max(a[],i);
 }
 void dec(int a[])//高精度a-1
 {
     int i;
     for(i=,a[]--;a[i]<;i++)
         a[i]+=J,a[i+]--;
     while(!a[a[]] && a[]>)a[]--;
 }
 int gjdbigger(int a[],int b[])//高精度a和高精度b比大小  1:a>b  0:a<b  -1:a=b
 {
     if(a[]!=b[])return a[]>b[];
     int i;
     for(i=a[];i;i--)if(a[i]!=b[i])return a[i]>b[i];
     return -;
 }
 bool tou(int i,int len)//判断头部是否满足
 {
     int ii;
     dec(a);
     for(ii=;ii<=a[] && i-ii;ii++)
         if(a[ii]!=s[i-ii]-'')return ;
     return ;
 }
 void gjdjia(int a[],int b[],int c[])//高精度a+高精度b，答案存在高精度c里
 {
     int i;
     mem(t,);
     t[]=max(a[],b[]);
     for(i=;i<=t[];i++)
         t[i]=a[i]+b[i];
     for(i=;i<=t[];i++)
         t[i+]+=t[i]/J,t[i]%=J;
     while(t[t[]+])t[]++;
     while(!t[t[]] && t[]>)t[]--;
     memcpy(c,t,sizeof(t));
 }
 void gjdjian(int a[],int b[],int c[])//高精度a-高精度b，答案存在高精度c里
 {
     int i;
     mem(t,);
     t[]=a[];
     for(i=;i<=t[];i++)
         t[i]=a[i]-b[i];
     for(i=;i<=t[];i++)
         if(t[i]<)t[i+]--,t[i]+=J;
     for(i=;i<=t[];i++)
         t[i+]+=t[i]/J,t[i]%=J;
     while(t[t[]+])t[]++;
     while(!t[t[]] && t[]>)t[]--;
     memcpy(c,t,sizeof(t));
 }
 void gjdchengdjd(int a[],int b,int c[])//高精度a*单精度b，答案存在高精度c里
 {
     int i;
     mem(t,);
     t[]=a[]+;
     for(i=;i<=a[];i++)
         t[i]=a[i]*b;
     for(i=;i<=t[];i++)
         t[i+]+=t[i]/J,t[i]%=J;
     while(t[t[]+])t[]++;
     while(!t[t[]] && t[]>)t[]--;
     memcpy(c,t,sizeof(t));
 }
 void cal(int len,int head)//已知a为第一个完整出现的数字，len为长度，head为出现的位置
 {
     int i,j;
     int tmp[N];
     mem(tmp,);
     //gjdchange(a,head,head+len-1);
     b[]=;b[]=;c[]=;c[]=;
     for(i=;i<a[];i++)
     {
         gjdchengdjd(b,i,tmp);
         gjdjia(c,tmp,c);
         b[b[]+]=;b[b[]++]=;
     }
     b[b[]]=;
     b[]=a[];b[b[]]=;
     for(i=;i<a[];i++)b[i]=;
     gjdjian(a,b,a);
     gjdchengdjd(a,len,b);
     gjdjia(c,b,c);
     mem(b,);
     b[]=;b[]=head-;
     while(b[b[]]>=J)b[b[]+]=b[b[]]/J,b[b[]++]%=J;
     gjdjian(c,b,c);
 }
 void work(int &len,int &head)//从中间断开分成两个数字的情况
 {
     int i,j,l;
     if(len<=n)
         gjdchange(a,head,head+len-);
     else a[]=MAX;
     b[]=;
     for(l=;l<=(n+)/+;l++)//最终长度n-l
     {
         for(j=;j<=l;j++)
             if(s[j]!=s[n-l+j])break;
         if(j<=l)continue;
         for(i=l+;i<=n-l;i++)
         {
             b[]=;
             for(j=i-;j;j--)b[++b[]]=s[j]-'';
             for(j=n-l;j>=i;j--)b[++b[]]=s[j]-'';
             plu(b);
             if(b[b[]]==)continue;
             for(j=;j<b[] && i+j<=n;j++)
                 if(b[b[]-j]!=s[i+j]-'')break;
             if(j<b[] && i+j<=n)continue;
             if(gjdbigger(a,b)){memcpy(a,b,sizeof(b));len=n-l;head=i;}
         }
     }
 }
 bool all9()//长度大于2且全是9的特殊情况 拆成8999999999...+9000000000...
 {
     int i;
     gjdchange(b,,n);
     for(i=;i<=b[];i++)if(b[i]!=)return ;
     if(b[]==)return ;
     b[b[]]=;plu(b);
     memcpy(a,b,sizeof(b));
     return ;
 }
 bool all0()//全是0的特殊情况 1+00000000...
 {
     int i;
     gjdchange(b,,n);
     for(i=;i<=b[];i++)if(b[i]!=)return ;
     b[++b[]]=;
     memcpy(a,b,sizeof(b));
     return ;
 }
 void gjdprint(int a[])//输出高精度a
 {
     int i;
     printf("%d",a[a[]]);
     for(i=a[]-;i;i--)
         printf("%d",a[i]);
     puts("");
 }
 int main()
 {
     #ifndef ONLINE_JUDGE
     freopen("1.txt","r",stdin);
 //    freopen("2.txt","w",stdout);
     #endif
     int i,j,k,l;
 //    for(scanf("%d",&cas);cas;cas--)
 //    for(scanf("%d",&cas),cass=1;cass<=cas;cass++)
     while(~scanf("%s",s+))
 //    while(~scanf("%d",&n))
     {
         s[]='.';
         n=strlen(s)-;
         for(l=;l<=n;l++)
         {
             for(i=;i<=l && i+l-<=n;i++)
             {
                 gjdchange(a,i,i+l-);
                 if(a[a[]]==)continue;//前导0不符合要求
                 if(i!=)
                 {
                     if(!tou(i,l))continue;
                     plu(a);
                 }
                 for(j=i+l;j<=n;j+=a[])
                 {
                     plu(a);
                     for(k=;j+k<=n && k<a[];k++)
                     {
                         if(a[a[]-k]!=s[j+k]-'')break;
                     }
                     if(j+k>n || k<a[])break;
                 }
                 if(j+k>n)break;
                 if(j>n)break;
             }
             if(i<=l && i+l-<=n)break;
         }
         //以上是将字符串拆成三段枚举，头，尾和中间的数字
         work(l,i);
         if(all9())cal(n,n);
         else if(all0())cal(n+,);
         else cal(l,i);
         gjdprint(c);
     }
     return ;
 }
 /*
 //

 //
 */