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. linux(模糊批量删除文件)删除指定文件夹中某个文件除外的其他文件

    # shopt -s extglob# rm -fr !(file1)如果是多个要排除的,可以这样:# rm -rf !(file1|file2)首先科普下shopt -s extglobBash S ...

  2. 浅谈ajax同步、异步的问题

    最近实习的时候看到过firefox的同步.异步的警告,想着概念不是那么清楚,于是整理了一下ajax同步异步方面的知识.我是小白,做个笔记. 首先就是概念问题,ajax根据async进行区分同步和异步过 ...

  3. 苏州地区--校招IT公司

    完整经历了苏州的秋招和春招,在本校和苏州大学跑了许多次的宣讲会,自认为对苏州IT企业的校招有一个充分的认知.原本打算在苏州找一份Java开发的工作,可是发现自己简历连那些公司的简历关都过不去(对双非学 ...

  4. maven:新建的maven工程需要添加一下插件

    //下面的代码是改变maven的默认使用Java版本,本机是jdk1.7所以用1.7<plugins> <plugin> <groupId>org.apache.m ...

  5. Virtual DOM的简单实现

    了解React的同学都知道,React提供了一个高效的视图更新机制:Virtual DOM,因为DOM天生就慢,所以操作DOM的时候要小心翼翼,稍微改动就会触发重绘重排,大量消耗性能. 1.Virtu ...

  6. Chapter 1 概述

    软件分为通用软件和应用软件,具有复杂.不可见和不断变化的特点:从出现到现在共经历了四个重要的发展阶段,但现在的软件发展依然存在着许多问题. 软件工程是采用工程的概念.原理.技术和方法来开发与维护软件, ...

  7. “吃神么,买神么”的第一个Sprint计划

    一.现状 我们这个团队刚接触网络网站的制作,前台后台的链接,数据库链接等还刚刚了解,在制作过程中药边学习边制作. 持着一个尽力做的心. 二.部分需求索引卡 主要的任务是把我们的主页面先大概做出来 三. ...

  8. jquery ajax 跨域问题

    前言:关于跨域CORS 1.没有跨域时,ajax默认是带cookie的 2.跨域时,两种解决方案: 1)服务器端在filter中配置 详情:http://blog.csdn.net/wzl002/ar ...

  9. BETA-6

    前言 我们居然又冲刺了·六 团队代码管理github 站立会议 队名:PMS 530雨勤(组长) 过去两天完成了哪些任务 新方案代码比之前的更简单,但是对场景的要求相应变高了,已经实现,误差感人 代码 ...

  10. iOS UIView性能最优的设计圆角并且绘制边框颜色

    //以给cell切圆角为例- (void)collectionView:(UICollectionView *)collectionView willDisplayCell:(UICollection ...