2023-08-18：用go写算法。你会得到一个字符串 text，
你应该把它分成 k 个子字符串 (subtext1, subtext2，…， subtextk)。
要求满足:
subtexti 是 非空 字符串，
所有子字符串的连接等于 text ，
( 即subtext1 + subtext2 + … + subtextk == text )，
subtexti == subtextk - i + 1 表示所有 i 的有效值( 即 1 <= i <= k )。
返回k可能最大值。
输入：text = “ghiabcdefhelloadamhelloabcdefghi”。
输出：7。
解释：我们可以把字符串拆分成 “(ghi)(abcdef)(hello)(adam)(hello)(abcdef)(ghi)”。
来自小红书、谷歌、Bloomberg、微软、亚马逊、字节跳动、摩根大通、Uber。
答案2023-08-18：
这两种算法的步骤可以概括如下：
方法1：longestDecomposition1(text)
1.初始化变量 s 为文本的字节数组，变量 n 为文本长度，变量 l 和 r 分别指向字节数组的开头和末尾，变量 ans 初始化为 0。
2.使用双指针法，当 l 小于等于 r 时进行循环：
-
初始化变量
size为 1。 -
在一个循环内，通过比较子字符串是否相同来确定
size的值，直到size使得剩余的子字符串无法再被分割为相同的子字符串。 -
如果
size使得剩余的子字符串可以被分割为相同的子字符串，则将ans增加 2。 -
向右移动
l和向左移动r，使其分别指向剩余子字符串的新开头和新结尾。
3.如果 r 恰好等于 l-1，则返回 ans，否则返回 ans + 1。
方法2：longestDecomposition2(text)
1.初始化变量 s 为文本的字节数组，变量 n 为文本长度，通过调用 generateDC3 函数生成 dc3 结构体。
2.初始化变量 rank 为 dc3 结构体中的 Rank 切片，初始化变量 rmq 为通过调用 NewRMQ 函数生成的 RMQ 结构体，变量 l 和 r 分别指向字节数组的开头和末尾，变量 ans 初始化为 0。
3.使用双指针法，当 l 小于等于 r 时进行循环：
-
初始化变量
size为 1。 -
在一个循环内，通过比较子字符串是否相同来确定
size的值，直到size使得剩余的子字符串无法再被分割为相同的子字符串。 -
如果
size使得剩余的子字符串可以被分割为相同的子字符串，则将ans增加 2。 -
向右移动
l和向左移动r，使其分别指向剩余子字符串的新开头和新结尾。
4.如果 r 恰好等于 l-1，则返回 ans，否则返回 ans + 1。
复杂度分析：
最初的算法longestDecomposition1的时间复杂度和空间复杂度如下：
-
时间复杂度：O(n^2)，其中n为输入字符串text的长度。
-
空间复杂度：O(n)，需要额外的空间来存储中间计算结果。
优化后的算法longestDecomposition2使用了DC3和RMQ的数据结构，并进行了分块处理，因此时间复杂度和空间复杂度如下：
-
时间复杂度：O(n log n)，其中n为输入字符串text的长度。
-
空间复杂度：O(n)，需要额外的空间来存储中间计算结果。
总结：优化后的算法longestDecomposition2虽然时间复杂度更低，但空间复杂度和初次计算的算法longestDecomposition1相同，都为O(n)。
go语言完整代码如下：
package main
import (
"fmt"
"strings"
"time"
)
func longestDecomposition1(text string) int {
if len(text) == 1 {
return 1
}
s := []byte(text)
n := len(s)
l := 0
r := n - 1
ans := 0
for l <= r {
size := 1
for 2*size <= r-l+1 {
if same1(s, l, r, size) {
break
}
size++
}
if 2*size <= r-l+1 {
ans += 2
}
l += size
r -= size
}
if r == l-1 {
return ans
}
return ans + 1
}
func same1(s []byte, l, r, size int) bool {
for i, j := l, r-size+1; j <= r; i, j = i+1, j+1 {
if s[i] != s[j] {
return false
}
}
return true
}
func longestDecomposition2(str string) int {
if len(str) == 1 {
return 1
}
s := []byte(str)
n := len(s)
dc3 := generateDC3(s, n)
rank := dc3.Rank
rmq := NewRMQ(dc3.Height)
l := 0
r := n - 1
ans := 0
for l <= r {
size := 1
for ; 2*size <= r-l+1; size++ {
if same2(rank, rmq, l, r, size) {
break
}
}
if 2*size <= r-l+1 {
ans += 2
}
l += size
r -= size
}
if r == l-1 {
return ans
}
return ans + 1
}
func same2(rank []int, rmq *RMQ, l, r, size int) bool {
start1 := l
start2 := r - size + 1
minStart := getMin(rank[start1], rank[start2])
maxStart := getMax(rank[start1], rank[start2])
return rmq.Min(minStart+1, maxStart) >= size
}
func getMin(a, b int) int {
if a < b {
return a
} else {
return b
}
}
func getMax(a, b int) int {
if a > b {
return a
} else {
return b
}
}
func generateDC3(s []byte, n int) *DC3 {
min0 := 2147483647
max0 := -2147483648
for _, cha := range s {
min0 = getMin(min0, int(cha))
max0 = getMax(max0, int(cha))
}
arr := make([]int, n)
for i := 0; i < n; i++ {
arr[i] = int(s[i]) - min0 + 1
}
return NewDC3(arr, max0-min0+1)
}
type DC3 struct {
Sa []int
Rank []int
Height []int
}
func NewDC3(nums []int, max int) *DC3 {
res := &DC3{}
res.Sa = res.sa(nums, max)
res.Rank = res.rank()
res.Height = res.height(nums)
return res
}
func (dc3 *DC3) sa(nums []int, max int) []int {
n := len(nums)
arr := make([]int, n+3)
copy(arr, nums)
return dc3.skew(arr, n, max)
}
func (dc3 *DC3) skew(nums []int, n int, K int) []int {
n0 := (n + 2) / 3
n1 := (n + 1) / 3
n2 := n / 3
n02 := n0 + n2
s12 := make([]int, n02+3)
sa12 := make([]int, n02+3)
for i, j := 0, 0; i < n+(n0-n1); i++ {
if i%3 != 0 {
s12[j] = i
j++
}
}
dc3.radixPass(nums, s12, sa12, 2, n02, K)
dc3.radixPass(nums, sa12, s12, 1, n02, K)
dc3.radixPass(nums, s12, sa12, 0, n02, K)
name := 0
c0 := -1
c1 := -1
c2 := -1
for i := 0; i < n02; i++ {
if c0 != nums[sa12[i]] || c1 != nums[sa12[i]+1] || c2 != nums[sa12[i]+2] {
name++
c0 = nums[sa12[i]]
c1 = nums[sa12[i]+1]
c2 = nums[sa12[i]+2]
}
if sa12[i]%3 == 1 {
s12[sa12[i]/3] = name
} else {
s12[sa12[i]/3+n0] = name
}
}
if name < n02 {
sa12 = dc3.skew(s12, n02, name)
for i := 0; i < n02; i++ {
s12[sa12[i]] = i + 1
}
} else {
for i := 0; i < n02; i++ {
sa12[s12[i]-1] = i
}
}
s0 := make([]int, n0)
sa0 := make([]int, n0)
for i, j := 0, 0; i < n02; i++ {
if sa12[i] < n0 {
s0[j] = 3 * sa12[i]
j++
}
}
dc3.radixPass(nums, s0, sa0, 0, n0, K)
sasa := make([]int, n)
for p, t, k := 0, n0-n1, 0; k < n; k++ {
i := 0
if sa12[t] < n0 {
i = sa12[t]*3 + 1
} else {
i = (sa12[t]-n0)*3 + 2
}
j := sa0[p]
rr := false
if sa12[t] < n0 {
rr = dc3.leq1(nums[i], s12[sa12[t]+n0], nums[j], s12[j/3])
} else {
rr = dc3.leq(nums[i], nums[i+1], s12[sa12[t]-n0+1], nums[j], nums[j+1], s12[j/3+n0])
}
if rr {
sasa[k] = i
t++
if t == n02 {
for k++; p < n0; p, k = p+1, k+1 {
sasa[k] = sa0[p]
}
}
} else {
sasa[k] = j
p++
if p == n0 {
for k++; t < n02; t, k = t+1, k+1 {
if sa12[t] < n0 {
sasa[k] = sa12[t]*3 + 1
} else {
sasa[k] = (sa12[t]-n0)*3 + 2
}
}
}
}
}
return sasa
}
func (t *DC3) radixPass(nums []int, input []int, output []int, offset int, n int, k int) {
cnt := make([]int, k+1)
for i := 0; i < n; i++ {
cnt[nums[input[i]+offset]]++
}
for i, sum := 0, 0; i < len(cnt); i++ {
t := cnt[i]
cnt[i] = sum
sum += t
}
for i := 0; i < n; i++ {
output[cnt[nums[input[i]+offset]]] = input[i]
cnt[nums[input[i]+offset]]++
}
}
func (t *DC3) leq1(a1 int, a2 int, b1 int, b2 int) bool {
return a1 < b1 || (a1 == b1 && a2 <= b2)
}
func (t *DC3) leq(a1 int, a2 int, a3 int, b1 int, b2 int, b3 int) bool {
return a1 < b1 || (a1 == b1 && t.leq1(a2, a3, b2, b3))
}
func (t *DC3) rank() []int {
n := len(t.Sa)
ans := make([]int, n)
for i := 0; i < n; i++ {
ans[t.Sa[i]] = i
}
return ans
}
func (t *DC3) height(s []int) []int {
n := len(s)
ans := make([]int, n)
for i, k := 0, 0; i < n; i++ {
if t.Rank[i] != 0 {
if k > 0 {
k--
}
j := t.Sa[t.Rank[i]-1]
for i+k < n && j+k < n && s[i+k] == s[j+k] {
k++
}
ans[t.Rank[i]] = k
}
}
return ans
}
type RMQ struct {
min0 [][]int
}
func NewRMQ(arr []int) *RMQ {
res := &RMQ{}
n := len(arr)
k := power2(n)
res.min0 = make([][]int, n+1)
for i := 0; i <= n; i++ {
res.min0[i] = make([]int, k+1)
}
for i := 1; i <= n; i++ {
res.min0[i][0] = arr[i-1]
}
for j := 1; (1 << j) <= n; j++ {
for i := 1; i+(1<<j)-1 <= n; i++ {
res.min0[i][j] = getMin(res.min0[i][j-1], res.min0[i+(1<<(j-1))][j-1])
}
}
return res
}
func (t *RMQ) Min(l, r int) int {
l++
r++
k := power2(r - l + 1)
return getMin(t.min0[l][k], t.min0[r-(1<<k)+1][k])
}
func power2(m int) int {
ans := 0
for (1 << ans) <= (m >> 1) {
ans++
}
return ans
}
func generateS(a, b int) string {
ans := make([]byte, a+b)
for i := 0; i < a; i++ {
ans[i] = 'a'
}
for i, j := a, 0; j < b; i, j = i+1, j+1 {
ans[i] = 'b'
}
return string(ans)
}
func generateT(part string, n int) string {
builder := strings.Builder{}
for i := 0; i < n; i++ {
builder.WriteString(part)
}
return builder.String()
}
func main() {
fmt.Println("先展示一下DC3的用法")
test := "aaabaaa"
dc3 := generateDC3([]byte(test), len(test))
fmt.Println("sa[i]表示字典序排名第i的是什么位置开头的后缀串")
sa := dc3.Sa
for i := 0; i < len(test); i++ {
fmt.Printf("%d : %d\n", i, sa[i])
}
fmt.Println("rank[i]表示i位置开头的后缀串的字典序排多少名")
rank := dc3.Rank
for i := 0; i < len(test); i++ {
fmt.Printf("%d : %d\n", i, rank[i])
}
fmt.Println("height[i]表示字典序排名i的后缀串和前一个排名的后缀串，最长公共前缀是多长")
height := dc3.Height
for i := 0; i < len(test); i++ {
fmt.Printf("%d : %d\n", i, height[i])
}
fmt.Println("性能测试开始")
var start, end int64
s := generateS(300000, 1)
t := generateT(s, 2)
start = time.Now().UnixNano() / int64(time.Millisecond)
fmt.Println("方法1的结果 :", longestDecomposition1(t))
end = time.Now().UnixNano() / int64(time.Millisecond)
fmt.Println("方法1的运行时间 :", end-start, "毫秒")
start = time.Now().UnixNano() / int64(time.Millisecond)
fmt.Println("方法2的结果 :", longestDecomposition2(t))
end = time.Now().UnixNano() / int64(time.Millisecond)
fmt.Println("方法2的运行时间 :", end-start, "毫秒")
fmt.Println("性能测试结束")
}