[CF1110H]Modest Substrings
description
CodeForces
定义一个正整数\(x\)是合适的当且仅当\(l\le x\le r\),其中\(l,r\le 10^{800}\)。
找到一个长度为\(n\)的数字串,使其包含合适的数作为子串的次数最多,\(n\le 2000\)。
如果有多解,输出字典序最小的那个。
solution
如果模式串的个数不多,那么直接套用AC自动机上数位dp的方法即可。
关于这一方法可参见[SDOI2014]数数。
现在这个做法的缺陷是要放入的串太多。
考虑简化。
根据数位dp的思想,一个固定位数的数只要达到安全态,后面的数码可以随意选择,
因此我们将\(\ge l\)或\(\le r\)的达到安全态的前缀都放入\(AC\)自动机;
设计\(dp\)状态为\(f[i][u]\),表示考虑了前\(i\)位,目前匹配到\(AC\)自动机的节点\(u\)时,已经能够匹配的子串个数。
对于后面的随意数码,由于题目限制长度为\(n\),所以只需要判断其长度是否\(>n-i-1\)即可。
code
实现代码长度居然达到\(3.5k\)...
#include<bits/stdc++.h>
#define pb push_back
#define mp make_pair
#define FL "a"
using namespace std;
typedef long long ll;
typedef long double dd;
const int N=2e3+10;
const int M=2e4+10;
const int inf=2147483647;
const dd pi=acos(-1);
const ll INF=1ll<<60;
inline ll read(){
ll data=0,w=1;char ch=getchar();
while(ch!='-'&&(ch<'0'||ch>'9'))ch=getchar();
if(ch=='-')w=-1,ch=getchar();
while(ch<='9'&&ch>='0')data=data*10+ch-48,ch=getchar();
return data*w;
}
inline void file(){
freopen(FL".in","r",stdin);
freopen(FL".out","w",stdout);
}
char sl[N],sr[N];
int n,lenl,lenr,cnt,vis[M][10],fail[M],g[M][N],f[N][M];
inline void upd(int &a,int b){a=a>b?a:b;}
inline void getpre(){
scanf("%s",sl+1);lenl=strlen(sl+1);
scanf("%s",sr+1);lenr=strlen(sr+1);
n=read();
if(lenl==lenr){
int ul,ur;ul=ur=0;
for(int i=1;i<=lenl;i++)
if(ul==ur){
for(int c=sl[i]-48+1;c<sr[i]-48;c++){
if(!vis[ul][c])vis[ul][c]=++cnt;
g[vis[ul][c]][lenl-i]++;
}
if(!vis[ul][sl[i]-48])vis[ul][sl[i]-48]=++cnt;
if(!vis[ur][sr[i]-48])vis[ur][sr[i]-48]=++cnt;
ul=vis[ul][sl[i]-48];ur=vis[ur][sr[i]-48];
}
else{
for(int c=sl[i]-48+1;c<10;c++){
if(!vis[ul][c])vis[ul][c]=++cnt;
g[vis[ul][c]][lenl-i]++;
}
for(int c=0;c<sr[i]-48;c++){
if(!vis[ur][c])vis[ur][c]=++cnt;
g[vis[ur][c]][lenr-i]++;
}
if(!vis[ul][sl[i]-48])vis[ul][sl[i]-48]=++cnt;
if(!vis[ur][sr[i]-48])vis[ur][sr[i]-48]=++cnt;
ul=vis[ul][sl[i]-48];ur=vis[ur][sr[i]-48];
}
g[ul][0]++;if(ul!=ur)g[ur][0]++;
}
else{
int u;u=0;
for(int i=1;i<=lenl;i++){
for(int c=sl[i]-48+1;c<10;c++){
if(!vis[u][c])vis[u][c]=++cnt;
g[vis[u][c]][lenl-i]++;
}
if(!vis[u][sl[i]-48])vis[u][sl[i]-48]=++cnt;
u=vis[u][sl[i]-48];
}
g[u][0]++;u=0;
for(int i=1;i<=lenr;i++){
for(int c=0;c<sr[i]-48;c++){
if(!vis[u][c])vis[u][c]=++cnt;
g[vis[u][c]][lenr-i]++;
}
if(!vis[u][sr[i]-48])vis[u][sr[i]-48]=++cnt;
u=vis[u][sr[i]-48];
}
g[u][0]++;
for(int i=lenl+1;i<lenr;i++)
for(int c=1;c<10;c++){
if(!vis[0][c])vis[0][c]=++cnt;
g[vis[0][c]][i-1]++;
}
}
for(int i=0;i<=cnt;i++)vis[0][0]=0;
}
inline void getfail(){
static queue<int>Q;while(!Q.empty())Q.pop();
for(int c=0;c<10;c++)if(vis[0][c])Q.push(vis[0][c]);
while(!Q.empty()){
int u=Q.front();Q.pop();
for(int c=0;c<10;c++){
int &v=vis[u][c];
if(v){
Q.push(v);fail[v]=vis[fail[u]][c];
for(int i=0;i<=n;i++)g[v][i]+=g[fail[v]][i];
}
else v=vis[fail[u]][c];
}
}
for(int u=0;u<=cnt;u++)
for(int i=0;i<=n;i++)
g[u][i]+=g[u][i-1];
}
inline void init(){
getpre();
getfail();
}
bool could[N][M];
inline void solve(){
memset(f,128,sizeof(f));f[0][0]=0;
for(int i=0;i<=n;i++)
for(int u=0;u<=cnt;u++)
if(f[i][u]>=0){
f[i][u]+=g[u][n-i];
for(int c=0;c<10;c++)
upd(f[i+1][vis[u][c]],f[i][u]);
}
int ans=0;
for(int u=0;u<=cnt;u++)upd(ans,f[n][u]);
printf("%d\n",ans);
for(int u=0;u<=cnt;u++)if(f[n][u]==ans)could[n][u]=1;
for(int i=n-1;~i;i--)
for(int u=0;u<=cnt;u++)
if(f[i][u]>=0)
for(int c=0;c<10;c++)
if(could[i+1][vis[u][c]]&&
f[i+1][vis[u][c]]==f[i][u]+g[vis[u][c]][n-i-1]){
could[i][u]=1;break;
}
assert(could[0][0]==1);
int u=0;
for(int i=0;i<n;i++)
for(int c=0;c<10;c++)
if(could[i+1][vis[u][c]]&&f[i+1][vis[u][c]]==f[i][u]+g[vis[u][c]][n-i-1]){
putchar(48+c);u=vis[u][c];break;
}
}
int main()
{
init();
solve();
return 0;
}
[CF1110H]Modest Substrings的更多相关文章
- CF1110H Modest Substrings AC自动机、DP
传送门 如果\(r-l\)比较小,可以将所有满足条件的串扔进\(AC\)自动机然后在上面DP,从前往后确定字符串的每一位. 但是\(l,r \leq 10^{800}\)就十分不可行,所以需要优化这个 ...
- CodeForces 1110H. Modest Substrings
题目简述:给定$1 \leq l \leq r \leq 10^{800}$,求一个长度为$n \leq 2000$的数字串$s$,其含有最多的[好]子串.一个串$s$是[好]的,如果将其看做数字时无 ...
- CodeForces Global Round 1
CodeForces Global Round 1 CF新的比赛呢(虽然没啥区别)!这种报名的人多的比赛涨分是真的快.... 所以就写下题解吧. A. Parity 太简单了,随便模拟一下就完了. B ...
- [LeetCode] Unique Substrings in Wraparound String 封装字符串中的独特子字符串
Consider the string s to be the infinite wraparound string of "abcdefghijklmnopqrstuvwxyz" ...
- Leetcode: Unique Substrings in Wraparound String
Consider the string s to be the infinite wraparound string of "abcdefghijklmnopqrstuvwxyz" ...
- CSU-1632 Repeated Substrings (后缀数组)
Description String analysis often arises in applications from biology and chemistry, such as the stu ...
- CF451D Count Good Substrings (DP)
Codeforces Round #258 (Div. 2) Count Good Substrings D. Count Good Substrings time limit per test 2 ...
- LA4671 K-neighbor substrings(FFT + 字符串Hash)
题目 Source http://acm.hust.edu.cn/vjudge/problem/19225 Description The Hamming distance between two s ...
- 后缀数组---New Distinct Substrings
Description Given a string, we need to find the total number of its distinct substrings. Input T- nu ...
随机推荐
- linux(模糊批量删除文件)删除指定文件夹中某个文件除外的其他文件
# shopt -s extglob# rm -fr !(file1)如果是多个要排除的,可以这样:# rm -rf !(file1|file2)首先科普下shopt -s extglobBash S ...
- 浅谈ajax同步、异步的问题
最近实习的时候看到过firefox的同步.异步的警告,想着概念不是那么清楚,于是整理了一下ajax同步异步方面的知识.我是小白,做个笔记. 首先就是概念问题,ajax根据async进行区分同步和异步过 ...
- 苏州地区--校招IT公司
完整经历了苏州的秋招和春招,在本校和苏州大学跑了许多次的宣讲会,自认为对苏州IT企业的校招有一个充分的认知.原本打算在苏州找一份Java开发的工作,可是发现自己简历连那些公司的简历关都过不去(对双非学 ...
- maven:新建的maven工程需要添加一下插件
//下面的代码是改变maven的默认使用Java版本,本机是jdk1.7所以用1.7<plugins> <plugin> <groupId>org.apache.m ...
- Virtual DOM的简单实现
了解React的同学都知道,React提供了一个高效的视图更新机制:Virtual DOM,因为DOM天生就慢,所以操作DOM的时候要小心翼翼,稍微改动就会触发重绘重排,大量消耗性能. 1.Virtu ...
- Chapter 1 概述
软件分为通用软件和应用软件,具有复杂.不可见和不断变化的特点:从出现到现在共经历了四个重要的发展阶段,但现在的软件发展依然存在着许多问题. 软件工程是采用工程的概念.原理.技术和方法来开发与维护软件, ...
- “吃神么,买神么”的第一个Sprint计划
一.现状 我们这个团队刚接触网络网站的制作,前台后台的链接,数据库链接等还刚刚了解,在制作过程中药边学习边制作. 持着一个尽力做的心. 二.部分需求索引卡 主要的任务是把我们的主页面先大概做出来 三. ...
- jquery ajax 跨域问题
前言:关于跨域CORS 1.没有跨域时,ajax默认是带cookie的 2.跨域时,两种解决方案: 1)服务器端在filter中配置 详情:http://blog.csdn.net/wzl002/ar ...
- BETA-6
前言 我们居然又冲刺了·六 团队代码管理github 站立会议 队名:PMS 530雨勤(组长) 过去两天完成了哪些任务 新方案代码比之前的更简单,但是对场景的要求相应变高了,已经实现,误差感人 代码 ...
- iOS UIView性能最优的设计圆角并且绘制边框颜色
//以给cell切圆角为例- (void)collectionView:(UICollectionView *)collectionView willDisplayCell:(UICollection ...