【模拟】Vijos P1005 超长数字串

时间:2023-03-10 07:38:27
【模拟】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 ;
}
/*
// //
*/