题目链接：
题目大意：
无限的正整数按顺序拼接成字符串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...优。
（不知道自己还有没有漏情况，希望没有吧，望指教。）
1 //
2 //by coolxxx
3 //#include<bits/stdc++.h>
4 #include<iostream>
5 #include<algorithm>
6 #include<string>
7 #include<iomanip>
8 #include<map>
9 #include<memory.h>
10 #include<time.h>
11 #include<stdio.h>
12 #include<stdlib.h>
13 #include<string.h>
14 //#include<stdbool.h>
15 #include<math.h>
16 #define min(a,b) ((a)<(b)?(a):(b))
17 #define max(a,b) ((a)>(b)?(a):(b))
18 #define abs(a) ((a)>0?(a):(-(a)))
19 #define lowbit(a) (a&(-a))
20 #define sqr(a) ((a)*(a))
21 #define swap(a,b) ((a)^=(b),(b)^=(a),(a)^=(b))
22 #define mem(a,b) memset(a,b,sizeof(a))
23 #define eps (1e-8)
24 #define J 10
25 #define mod 1000000007
26 #define MAX 0x7f7f7f7f
27 #define PI 3.14159265358979323
28 #define N 404
29 using namespace std;
30 typedef long long LL;
31 int cas,cass;
32 int n,m,lll,ans;
33 double anss;
34 char s[N];
35 int a[N],b[N],t[N],c[N];
36 void gjdchange(int a[],int l,int r)//l~r中间的字符转化为高精度数a
37 {
38 int i;
39 mem(t,0);
40 t[0]=r-l+1;
41 for(i=1;i<=r-l+1;i++)
42 t[i]=s[r-i+1]-‘0‘;
43 memcpy(a,t,sizeof(t));
44 }
45 void plu(int a[])//高精度a+1
46 {
47 int i;
48 for(i=1,a[1]++;a[i]>=J;i++)
49 a[i]-=J,a[i+1]++;
50 a[0]=max(a[0],i);
51 }
52 void dec(int a[])//高精度a-1
53 {
54 int i;
55 for(i=1,a[1]--;a[i]<0;i++)
56 a[i]+=J,a[i+1]--;
57 while(!a[a[0]] && a[0]>1)a[0]--;
58 }
59 int gjdbigger(int a[],int b[])//高精度a和高精度b比大小 1:a>b 0:a<b -1:a=b
60 {
61 if(a[0]!=b[0])return a[0]>b[0];
62 int i;
63 for(i=a[0];i;i--)if(a[i]!=b[i])return a[i]>b[i];
64 return -1;
65 }
66 bool tou(int i,int len)//判断头部是否满足
67 {
68 int ii;
69 dec(a);
70 for(ii=1;ii<=a[0] && i-ii;ii++)
71 if(a[ii]!=s[i-ii]-‘0‘)return 0;
72 return 1;
73 }
74 void gjdjia(int a[],int b[],int c[])//高精度a+高精度b，答案存在高精度c里
75 {
76 int i;
77 mem(t,0);
78 t[0]=max(a[0],b[0]);
79 for(i=1;i<=t[0];i++)
80 t[i]=a[i]+b[i];
81 for(i=1;i<=t[0];i++)
82 t[i+1]+=t[i]/J,t[i]%=J;
83 while(t[t[0]+1])t[0]++;
84 while(!t[t[0]] && t[0]>1)t[0]--;
85 memcpy(c,t,sizeof(t));
86 }
87 void gjdjian(int a[],int b[],int c[])//高精度a-高精度b，答案存在高精度c里
88 {
89 int i;
90 mem(t,0);
91 t[0]=a[0];
92 for(i=1;i<=t[0];i++)
93 t[i]=a[i]-b[i];
94 for(i=1;i<=t[0];i++)
95 if(t[i]<0)t[i+1]--,t[i]+=J;
96 for(i=1;i<=t[0];i++)
97 t[i+1]+=t[i]/J,t[i]%=J;
98 while(t[t[0]+1])t[0]++;
99 while(!t[t[0]] && t[0]>1)t[0]--;
100 memcpy(c,t,sizeof(t));
101 }
102 void gjdchengdjd(int a[],int b,int c[])//高精度a*单精度b，答案存在高精度c里
103 {
104 int i;
105 mem(t,0);
106 t[0]=a[0]+3;
107 for(i=1;i<=a[0];i++)
108 t[i]=a[i]*b;
109 for(i=1;i<=t[0];i++)
110 t[i+1]+=t[i]/J,t[i]%=J;
111 while(t[t[0]+1])t[0]++;
112 while(!t[t[0]] && t[0]>1)t[0]--;
113 memcpy(c,t,sizeof(t));
114 }
115 void cal(int len,int head)//已知a为第一个完整出现的数字，len为长度，head为出现的位置
116 {
117 int i,j;
118 int tmp[N];
119 mem(tmp,0);
120 //gjdchange(a,head,head+len-1);
121 b[0]=1;b[1]=9;c[0]=1;c[1]=0;
122 for(i=1;i<a[0];i++)
123 {
124 gjdchengdjd(b,i,tmp);
125 gjdjia(c,tmp,c);
126 b[b[0]+1]=9;b[b[0]++]=0;
127 }
128 b[b[0]]=0;
129 b[0]=a[0];b[b[0]]=1;
130 for(i=1;i<a[0];i++)b[i]=0;
131 gjdjian(a,b,a);
132 gjdchengdjd(a,len,b);
133 gjdjia(c,b,c);
134 mem(b,0);
135 b[0]=1;b[1]=head-2;
136 while(b[b[0]]>=J)b[b[0]+1]=b[b[0]]/J,b[b[0]++]%=J;
137 gjdjian(c,b,c);
138 }
139 void work(int &len,int &head)//从中间断开分成两个数字的情况
140 {
141 int i,j,l;
142 if(len<=n)
143 gjdchange(a,head,head+len-1);
144 else a[0]=MAX;
145 b[0]=0;
146 for(l=0;l<=(n+1)/2+1;l++)//最终长度n-l
147 {
148 for(j=1;j<=l;j++)
149 if(s[j]!=s[n-l+j])break;
150 if(j<=l)continue;
151 for(i=l+2;i<=n-l;i++)
152 {
153 b[0]=0;
154 for(j=i-1;j;j--)b[++b[0]]=s[j]-‘0‘;
155 for(j=n-l;j>=i;j--)b[++b[0]]=s[j]-‘0‘;
156 plu(b);
157 if(b[b[0]]==0)continue;
158 for(j=0;j<b[0] && i+j<=n;j++)
159 if(b[b[0]-j]!=s[i+j]-‘0‘)break;
160 if(j<b[0] && i+j<=n)continue;
161 if(gjdbigger(a,b)){memcpy(a,b,sizeof(b));len=n-l;head=i;}
162 }
163 }
164 }
165 bool all9()//长度大于2且全是9的特殊情况 拆成8999999999...+9000000000...
166 {
167 int i;
168 gjdchange(b,1,n);
169 for(i=1;i<=b[0];i++)if(b[i]!=9)return 0;
170 if(b[0]==1)return 0;
171 b[b[0]]=8;plu(b);
172 memcpy(a,b,sizeof(b));
173 return 1;
174 }
175 bool all0()//全是0的特殊情况 1+00000000...
176 {
177 int i;
178 gjdchange(b,1,n);
179 for(i=1;i<=b[0];i++)if(b[i]!=0)return 0;
180 b[++b[0]]=1;
181 memcpy(a,b,sizeof(b));
182 return 1;
183 }
184 void gjdprint(int a[])//输出高精度a
185 {
186 int i;
187 printf("%d",a[a[0]]);
188 for(i=a[0]-1;i;i--)
189 printf("%d",a[i]);
190 puts("");
191 }
192 int main()
193 {
194 #ifndef ONLINE_JUDGE
195 freopen("1.txt","r",stdin);
196 // freopen("2.txt","w",stdout);
197 #endif
198 int i,j,k,l;
199 // for(scanf("%d",&cas);cas;cas--)
200 // for(scanf("%d",&cas),cass=1;cass<=cas;cass++)
201 while(~scanf("%s",s+1))
202 // while(~scanf("%d",&n))
203 {
204 s[0]=‘.‘;
205 n=strlen(s)-1;
206 for(l=1;l<=n;l++)
207 {
208 for(i=1;i<=l && i+l-1<=n;i++)
209 {
210 gjdchange(a,i,i+l-1);
211 if(a[a[0]]==0)continue;//前导0不符合要求
212 if(i!=1)
213 {
214 if(!tou(i,l))continue;
215 plu(a);
216 }
217 for(j=i+l;j<=n;j+=a[0])
218 {
219 plu(a);
220 for(k=0;j+k<=n && k<a[0];k++)
221 {
222 if(a[a[0]-k]!=s[j+k]-‘0‘)break;
223 }
224 if(j+k>n || k<a[0])break;
225 }
226 if(j+k>n)break;
227 if(j>n)break;
228 }
229 if(i<=l && i+l-1<=n)break;
230 }
231 //以上是将字符串拆成三段枚举，头，尾和中间的数字
232 work(l,i);
233 if(all9())cal(n,n);
234 else if(all0())cal(n+1,0);
235 else cal(l,i);
236 gjdprint(c);
237 }
238 return 0;
239 }
240 /*
241 //
242
243 //
244 */
时间： 2024-11-05 21:39:12