原文链接https://www.cnblogs.com/zhouzhendong/p/UOJ172.html

题解

首先,这个问题显然是个背包问题。

然后,可以证明:一个字符串的 border 长度可以划分成 $O(log |S|)$ 个等差数列。

(以下图片摘自  金策 - 《字符串算法选讲》)

由于长度 n 可以随便取,所以我们可以在对n取模的意义下做背包,设 dis[i] 为 占用背包容量%n = i 时至少要占用多少背包容量,那么直接建 $n^2$ 条边跑一下 dijkstra 就可以在 $O(n^2\log n )$ 的时间复杂度内解决此题。

对于所有的长度,我们将其分组,保证每一组是一个等差数列。

我们枚举等差数列,依次将其加入背包。

假设最终背包中最多包含当前等差数列的一个元素,那么比较好做,直接单调队列就好了。

但是可以有多个当前等差数列的元素计入最终背包容量。也就是说,假设当前等差数列用 $v + kd$ 表示,那么可能会有 $tv + \sum k_i d$ 这种情况。这在对 n 取模的意义下非常难做。

于是我们考虑把它转化成对 $v$ 取模的意义下的结果,这样就方便多了。

考虑如何转化:先把原先在对 n 取模意义下的结果放到对 v 取模意义下的数组里,然后考虑每一个 $dis'[i]$ 都可以更新 $dis'[(i+n)\mod v ]$,即 $dis'[(i+n)\mod v] =min(dis'[(i+n)\mod v],dis'[i] + n)$ 。

把转移形成的环搞出来,dp转移一下就可以得到新的dis数组。

剩下的单调队列+dp也类似。

于是时间复杂度是 $O(n\log n)$ 。

代码

#pragma GCC optimize("Ofast","inline")
#include <bits/stdc++.h>
#define clr(x) memset(x,0,sizeof (x))
#define fi first
#define se second
using namespace std;
typedef long long LL;
LL read(){
LL x=0,f=0;
char ch=getchar();
while (!isdigit(ch))
f|=ch=='-',ch=getchar();
while (isdigit(ch))
x=(x<<1)+(x<<3)+(ch^48),ch=getchar();
return f?-x:x;
}
const int N=500005;
const LL INF=2e18;
int T,n,m;
LL w;
int Fail[N],nf[N];
struct brd{
int v,d,c;
}b[N];
char s[N];
LL dis[N],rng;
LL gcd(LL a,LL b){
return b?gcd(b,a%b):a;
}
void Trans(LL nr){
static LL d[N];
for (int i=0;i<nr;i++)
d[i]=INF;
for (int i=0;i<rng;i++)
d[dis[i]%nr]=min(d[dis[i]%nr],dis[i]);
int g=gcd(rng,nr);
for (int t=0;t<g;t++){
int e=nr/g*2,_i=t,rg=rng%nr;
for (LL i=t,j=0;j<e;i+=rng,j++){
d[_i]=min(d[_i],d[_i<rg?_i-rg+nr:_i-rg]+rng);
if ((_i+=rg)>=nr)
_i-=nr;
}
}
rng=nr;
for (int i=0;i<rng;i++)
dis[i]=d[i];
}
pair <LL,LL> Q[N*2];
int head,tail;
void DP(int d,int c){
int g=gcd(rng,d);
for (int t=0;t<g;t++){
head=1,tail=0;
int e=rng/g*2,_i=t;
for (LL i=t,j=0;j<e;i+=d,j++){
if (head<=tail&&Q[head].fi+(LL)c*d<=i)
head++;
if (head<=tail)
dis[_i]=min(dis[_i],Q[head].se+rng+i);
while (head<=tail&&Q[tail].se>=dis[_i]-i)
tail--;
Q[++tail]=make_pair(i,dis[_i]-i);
if ((_i+=d)>=rng)
_i-=rng;
}
}
}
int main(){
T=read();
while (T--){
n=read(),w=read();
scanf("%s",s+1);
Fail[1]=0;
for (int i=2;i<=n;i++){
int k=Fail[i-1];
while (k>0&&s[k+1]!=s[i])
k=Fail[k];
if (s[k+1]==s[i])
k++;
Fail[i]=k;
}
m=0;
for (int i=n;i>0;i=Fail[i])
nf[Fail[i]]=i;
for (int j=0;j<n;j=nf[j]){
int i=n-j;
if (j==0)
m++,b[m].v=i,b[m].d=0,b[m].c=1;
else if (!b[m].d)
b[m].d=b[m].v-i,b[m].c++,b[m].v=i;
else if (b[m].v-b[m].d==i)
b[m].v=i,b[m].c++;
else
m++,b[m].v=i,b[m].d=0,b[m].c=1;
}
if (!b[m].d)
b[m].d=1;
clr(dis);
rng=n;
for (int i=0;i<n;i++)
dis[i]=INF;
dis[0]=n;
for (int i=1;i<=m;i++){
Trans(b[i].v);
DP(b[i].d,b[i].c);
}
LL ans=0;
for (int i=0;i<rng;i++)
if (dis[i]<=w)
ans+=w/rng-dis[i]/rng+(w%rng>=dis[i]%rng?1:0);
cout<<ans<<endl;
}
return 0;
}

  

UOJ#172. 【WC2016】论战捆竹竿 字符串 KMP 动态规划 单调队列 背包的更多相关文章

  1. bzoj4406: [Wc2016]论战捆竹竿&&uoj#172. 【WC2016】论战捆竹竿

    第二次在bzoj跑进前十竟然是因为在UOJ卡常致死 首先这个题其实就是一个无限背包 一般做法是同余最短路,就是bzoj2118: 墨墨的等式可以拿到30分的好成绩 背包是个卷积就分治FFT优化那么下面 ...

  2. Luogu4156 WC2016 论战捆竹竿 KMP、同余类最短路、背包、单调队列

    传送门 豪华升级版同余类最短路-- 官方题解 主要写几个小trick: \(1.O(nm)\)实现同余类最短路: 设某一条边长度为\(x\),那么我们选择一个点,在同余类上不断跳\(x\),可以形成一 ...

  3. luogu P4156 [WC2016]论战捆竹竿

    传送门 官方题解(证明都在这) 神仙题鸭qwq 转化模型,发现这题本质就是一个集合,每次可以加上集合里的数,问可以拼出多少不同的数 首先暴力需要膜意义下的最短路,例题戳这 然后这个暴力可以优化成N^2 ...

  4. BZOJ4406 WC2016 论战捆竹竿

    Problem BZOJ Solution 显然是一个同余系最短路问题,转移方案就是所有|S|-border的长度,有 \(O(n)\) 种,暴力跑dijkstra的复杂度为 \(O(n^2\log ...

  5. 「WC2016」论战捆竹竿

    「WC2016」论战捆竹竿 前置知识 参考资料:<论战捆竹竿解题报告-王鉴浩>,<字符串算法选讲-金策>. Border&Period 若前缀 \(pre(s,x)​\ ...

  6. 【WC2016】论战捆竹竿

    已经快三周了啊--终于把挖的坑填了-- 首先显然是把除了自身的所有border拿出来,即做 \(\left\{ n - b_1, n - b_2, \dots, n - b_k, n \right\} ...

  7. 【LuoguP4156】论战捆竹竿

    题目链接 题意简述 你有一个长度为 n 的字符串 , 将它复制任意次 , 复制出的串的前缀可以与之前的串的后缀重叠在一起 , 问最后总共可能的长度数目 , 长度不能超过 \(w\) 多组数据. \(n ...

  8. UOJ#172. 【WC2016】论战捆竹竿

    传送门 首先这个题目显然就是先求出所有的 \(border\),问题转化成一个可行性背包的问题 一个方法就是同余类最短路,裸跑 \(30\) 分,加优化 \(50\) 分 首先有个性质 \(borde ...

  9. hdu 5510 Bazinga(字符串kmp)

    Bazinga Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Others)Total Sub ...

随机推荐

  1. Flask 蓝图,数据库链接

    蓝图 使用场景 如果代码非常多,要进行归类.不同的功能放在不同的文件,把相关的视图函数也放进去. 蓝图也就是对flask的目录结构进行分配(应用于小,中型的程序) 当然对于大型项目也可以通过   ur ...

  2. Android popupwindow 失去焦点或者点击空白区域时消失的解决方法

    先来看下Android API 的这个Methods: public void setOutsideTouchable (boolean touchable) Controls whether the ...

  3. I2C 总线原理与架构

    一.I2C总线原理 I2C是一种常用的串行总线,由串行数据线SDA 和串行时钟线SCL组成.I2C是一种多主机控制总线,它和USB总线不同,USB是基于master-slave机制,任何设备的通信必须 ...

  4. Linux基本命令总结(四)

    接上篇: 16,locate 让使用者可以很快速的搜寻档案系统内是否有指定的档案.其方法是先建立一个包括系统内所有档案名称及路径的数据库,之后当寻找时就只需查询这个数据库,而不必实际深入档案系统之中了 ...

  5. Maven pom文件标签解析大全

    <span style="padding:0px; margin:0px"><project xmlns="http://maven.apache.or ...

  6. python之shelve模块详解

    一.定义 Shelve是对象持久化保存方法,将对象保存到文件里面,缺省(即默认)的数据存储文件是二进制的. 二.用途 可以作为一个简单的数据存储方案. 三.用法 使用时,只需要使用open函数获取一个 ...

  7. 驱动调试(一)-printk

    目录 驱动调试(一)-printk 引入 框架 入口console_setup add_preferred_console register_console s3c24xx_serial_initco ...

  8. SpringBoot项目@RestController使用 redirect 重定向无效

    Spring MVC项目中页面重定向一般使用return "redirect:/other/controller/";即可. 而Spring Boot当我们使用了@RestCont ...

  9. 我对SAP Business One 项目实施的理解

    一.什么是SAP: 大家都知道ERP是什么,ERP是企业资源计划管理系统.是指建立在信息技术基础上,集信息技术与先进管理思想于一身,以系统化的管理思想,为企业员工及决策层提供决策手段的管理平台.那么问 ...

  10. [Everyday Mathematics]20150305

    设 $f\in C^2[0,1]$, $$\bex f(0)=-1,\quad f'(1)=3,\quad \int_0^1 xf''(x)\rd x=1. \eex$$ 试求 $f(1)$. 解答: ...