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的更多相关文章

  1. CF1110H Modest Substrings AC自动机、DP

    传送门 如果\(r-l\)比较小,可以将所有满足条件的串扔进\(AC\)自动机然后在上面DP,从前往后确定字符串的每一位. 但是\(l,r \leq 10^{800}\)就十分不可行,所以需要优化这个 ...

  2. CodeForces 1110H. Modest Substrings

    题目简述:给定$1 \leq l \leq r \leq 10^{800}$,求一个长度为$n \leq 2000$的数字串$s$,其含有最多的[好]子串.一个串$s$是[好]的,如果将其看做数字时无 ...

  3. CodeForces Global Round 1

    CodeForces Global Round 1 CF新的比赛呢(虽然没啥区别)!这种报名的人多的比赛涨分是真的快.... 所以就写下题解吧. A. Parity 太简单了,随便模拟一下就完了. B ...

  4. [LeetCode] Unique Substrings in Wraparound String 封装字符串中的独特子字符串

    Consider the string s to be the infinite wraparound string of "abcdefghijklmnopqrstuvwxyz" ...

  5. Leetcode: Unique Substrings in Wraparound String

    Consider the string s to be the infinite wraparound string of "abcdefghijklmnopqrstuvwxyz" ...

  6. CSU-1632 Repeated Substrings (后缀数组)

    Description String analysis often arises in applications from biology and chemistry, such as the stu ...

  7. CF451D Count Good Substrings (DP)

    Codeforces Round #258 (Div. 2) Count Good Substrings D. Count Good Substrings time limit per test 2 ...

  8. LA4671 K-neighbor substrings(FFT + 字符串Hash)

    题目 Source http://acm.hust.edu.cn/vjudge/problem/19225 Description The Hamming distance between two s ...

  9. 后缀数组---New Distinct Substrings

    Description Given a string, we need to find the total number of its distinct substrings. Input T- nu ...

随机推荐

  1. nGrinder性能测试平台的安装部署

    1.从GitHub下载war包: https://github.com/naver/ngrinder/releases 2.把ngrinder-controller-3.4.2.war重命名为ngri ...

  2. mysql面试常见题目2

    Sutdent表的定义 字段名 字段描述 数据类型 主键 外键 非空 唯一 自增 Id 学号 INT(10) 是 否 是 是 是 sName 姓名 VARCHAR(20) 否 否 是 否 否 Sex ...

  3. 利用VS2015开发python版本的caffe应用

    打开VS2015,选择“新建项目”->“其它语言”->“python”,VS会提示你安装PTVS(Python Tools for Visual Studio)插件,安装完毕后即可开始py ...

  4. Shuffle Bags让你的随机不那么随机

    前言 当我最初写游戏时,我经常使用标准Random()函数,然后写一堆if和else条件来我获得预期结果.如果结果不太好,我会写更多的条件进行过滤或者筛选,直到我觉得游戏变得有趣.最近我发现有更好的方 ...

  5. vue-scroller实现vue单页面的上拉加载和下拉刷新问题

    在vue中如何简单的实现页面的上拉加载和下拉刷新,在这里我推荐使用vue-scrolle插件. vue-scrolle的基本使用方法: 1.下载 npm i vue-scroller -D 2.导包 ...

  6. JDBC Mysql 驱动连接异常

    在做JDBC连接Mysql的时候遇到了三个异常: 第一个是:mysql8.0 caching_sha2_password 这个异常是由于是因为在mysql8.0之前的密码规则是mysql_native ...

  7. socket 编程 : shutdown vs close

    TCP/IP 四次挥手 首先作者先描述一下TCP/IP 协议中四次挥手的过程,如果对此已经熟悉的读者可以跳过本节. 四次挥手 这是一个很经典的示例图,众所周知tcp socket 在一个生命周期中有很 ...

  8. Dede织梦验证码不显示,织梦后台登陆验证码不显示解决方法

    关于"织梦验证码不显示"的解决方法 "织梦验证码无法显示出来"的问题分析? 1.之前显示正常,但是换了服务器后就不能够正常显示:(这种通常是网站程序经过迁移后所 ...

  9. Hands on Machine Learning with Sklearn and TensorFlow学习笔记——机器学习概览

    一.什么是机器学习? 计算机程序利用经验E(训练数据)学习任务T(要做什么,即目标),性能是P(性能指标),如果针对任务T的性能P随着经验E不断增长,成为机器学习.[这是汤姆米切尔在1997年定义] ...

  10. openssl在多平台和多语言之间进行RSA加解密注意事项

    首先说一下平台和语言: 系统平台为CentOS6.3,RSA加解密时使用NOPADDING进行填充 1)使用C/C++调用系统自带的openssl 2)Android4.2模拟器,第三方openssl ...