找出第二个文本抄袭第一个文本的所有位置和长度 后缀数组 UVA 10526 - Intellectual Property

该博客介绍了如何使用后缀数组解决UVA 10526 - Intellectual Property问题,即找出第二个文本抄袭第一个文本的所有位置和长度。通过将两个文本拼接,利用后缀数组和高度数组的特性,按长度降序和位置升序输出前k个抄袭片段。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

题目链接

题意:给定两个问题,要求找出第二个文本抄袭第一个文本的所有位置和长度,输出前k个,按长度从大到小先排,长度一样的按位置从小到大

思路:后缀数组,把两个文本拼接起来,记录下拼接位置为tdp,这样如果sa[i] < tdp就是前面的文本开头,如果sa[i] >= tdp就是后面的文本开头,拼接起来的求出height数组,利用该数组的性质,从前往后扫一遍,从后往前扫一遍,把所有位置的最大值保存下来,最后在扫描一遍位置,把答案记录下来

代码:

#include <cstdio>  
#include <cstring>  
#include <algorithm>  
  
using namespace std;  
  
const int MAXLEN = 200005;  
const int INF = 0x3f3f3f3f;  
  
char str[55555];  
int k, tdp, an, v[MAXLEN];  
  
struct Ans {  
    int len, pos;  
    Ans() {}  
    Ans(int len, int pos) {  
        this->len = len;  
        this->pos = pos;  
    }  
} ans[MAXLEN];  
  
bool cmp(Ans a, Ans b) {  
    if (a.len == b.len) return a.pos < b.pos;  
    return a.len > b.len;  
}  
  
struct Suffix {  
  
    int s[MAXLEN];  
    int sa[MAXLEN], t[MAXLEN], t2[MAXLEN], c[MAXLEN], n;  
    int rank[MAXLEN], height[MAXLEN];  
  
    void build_sa(int m) {  
        n++;  
        int i, *x = t, *y = t2;  
        for (i = 0; i < m; i++) c[i] = 0;  
        for (i = 0; i < n; i++) c[x[i] = s[i]]++;  
        for (i = 1; i < m; i++) c[i] += c[i - 1];  
        for (i = n - 1; i >= 0; i--) sa[--c[x[i]]] = i;  
        for (int k = 1; k <= n; k <<= 1) {  
            int p = 0;  
            for (i = n - k; i < n; i++) y[p++] = i;  
            for (i = 0; i < n; i++) if (sa[i] >= k) y[p++] = sa[i] - k;  
            for (i = 0; i < m; i++) c[i] = 0;  
            for (i = 0; i < n; i++) c[x[y[i]]]++;  
            for (i = 0; i < m; i++) c[i] += c[i - 1];  
            for (i = n - 1; i >= 0; i--) sa[--c[x[y[i]]]] = y[i];  
            swap(x, y);  
            p = 1; x[sa[0]] = 0;  
            for (i = 1; i < n; i++)  
                x[sa[i]] = (y[sa[i - 1]] == y[sa[i]] && y[sa[i - 1] + k] == y[sa[i] + k]) ? p - 1 : p++;  
            if (p >= n) break;  
            m = p;  
        }  
        n--;  
    }  
  
    void getHeight() {  
        int i, j, k = 0;  
        for (i = 1; i <= n; i++) rank[sa[i]] = i;  
        for (i = 0; i < n; i++) {  
            if (k) k--;  
            int j = sa[rank[i] - 1];  
            while (s[i + k] == s[j + k]) k++;  
            height[rank[i]] = k;  
        }  
    }  
  
    void init() {  
        tdp = 0; n = 0; an = 0;  
        gets(str);  
        while (gets(str)) {  
            if (strcmp(str, "END TDP CODEBASE") == 0) break;  
            int len = strlen(str);  
            str[len] = '\n';  
            for (int i = 0; i <= len; i++)  
                s[n++] = str[i];  
        }  
        tdp = n;  
        s[n++] = 260;  
        gets(str);  
        while (gets(str)) {  
            if (strcmp(str, "END JCN CODEBASE") == 0) break;  
            int len = strlen(str);  
            str[len] = '\n';  
            for (int i = 0; i <= len; i++)  
                s[n++] = str[i];  
        }  
        s[n] = 0;  
    }  
  
    void solve() {  
        init();  
        build_sa(261);  
        getHeight();  
        memset(v, 0, sizeof(v));  
        int Min = -1;  
        for (int i = 1; i <= n; i++) {  
            if (sa[i] < tdp) Min = INF;  
            else if (sa[i] > tdp) {  
                if (Min == -1) continue;  
                Min = min(height[i], Min);  
                v[sa[i] - tdp - 1] = max(Min, v[sa[i] - tdp - 1]);  
            }  
        }  
        Min = -1;  
        for (int i = n; i >= 1; i--) {  
            if (sa[i] < tdp) Min = INF;  
            else if (sa[i] > tdp) {  
                if (Min == -1) continue;  
                Min = min(height[i + 1], Min);  
                v[sa[i] - tdp - 1] = max(Min, v[sa[i] - tdp - 1]);  
            }  
        }  
        int r = -1;  
        for (int i = 0; i < n - tdp; i++) {  
            if (i + v[i] <= r) continue;  
            if (v[i] == 0) continue;  
            ans[an++] = Ans(v[i], i);  
            r = i + v[i];  
        }  
        sort(ans, ans + an, cmp);  
        for (int i = 0; i < min(an, k); i++) {  
            printf("INFRINGING SEGMENT %d LENGTH %d POSITION %d\n", i + 1, ans[i].len, ans[i].pos);  
            for (int j = ans[i].pos + tdp + 1; j < ans[i].pos + tdp + 1 + ans[i].len; j++)  
                printf("%c", s[j]);  
            printf("\n");  
        }  
    }  
  
} gao;  
  
int main() {  
    int bo = 0;  
    int cas = 0;  
    while (~scanf("%d%*c", &k) && k) {  
        if (bo) printf("\n");  
        else bo = 1;  
        printf("CASE %d\n", ++cas);  
        gao.solve();  
    }  
    return 0;  
}  


内容概要:本文详细介绍了扫描单分子定位显微镜(scanSMLM)技术及其在三维超分辨体积成像中的应用。scanSMLM通过电调透镜(ETL)实现快速轴向扫描,结合4f检测系统将不同焦平面的荧光信号聚焦到固定成像面,从而实现快速、大视场的三维超分辨成像。文章不仅涵盖了系统硬件的设计与实现,还提供了详细的软件代码实现,包括ETL控制、3D样本模拟、体积扫描、单分子定位、3D重建分子聚类分析等功能。此外,文章还比较了循环扫描与常规扫描模式,展示了前者在光漂白效应上的优势,并通过荧光珠校准、肌动蛋白丝、线粒体网络流感A病毒血凝素(HA)蛋白聚类的三维成像实验,验证了系统的性能应用潜力。最后,文章深入探讨了HA蛋白聚类与病毒感染的关系,模拟了24小时内HA聚类的动态变化,提供了从分子到细胞尺度的多尺度分析能力。 适合人群:具备生物学、物理学或工程学背景,对超分辨显微成像技术感兴趣的科研人员,尤其是从事细胞生物学、病毒学或光学成像研究的科学家技术人员。 使用场景及目标:①理解掌握scanSMLM技术的工作原理及其在三维超分辨成像中的应用;②学习如何通过Python代码实现完整的scanSMLM系统,包括硬件控制、图像采集、3D重建数据分析;③应用于单分子水平研究细胞内结构动态过程,如病毒入侵机制、蛋白质聚类等。 其他说明:本文提供的代码不仅实现了scanSMLM系统的完整工作流程,还涵盖了多种超分辨成像技术的模拟比较,如STED、GSDIM等。此外,文章还强调了系统在硬件改动小、成像速度快等方面的优势,为研究人员提供了从理论到实践的全面指导。
内容概要:本文详细介绍了基于Seggiani提出的渣层计算模型,针对Prenflo气流床气化炉中炉渣的积累流动进行了模拟。模型不仅集成了三维代码以提供气化炉内部的温度浓度分布,还探讨了操作条件变化对炉渣行为的影响。文章通过Python代码实现了模型的核心功能,包括炉渣粘度模型、流动速率计算、厚度更新、与三维模型的集成以及可视化展示。此外,还扩展了模型以考虑炉渣组成对特性的影响,并引入了Bingham流体模型,更精确地描述了含未溶解颗粒的熔渣流动。最后,通过实例展示了氧气-蒸汽流量增加2%时的动态响应,分析了温度、流动特性渣层分布的变化。 适合人群:从事煤气化技术研究的专业人士、化工过程模拟工程师、以及对工业气化炉操作优化感兴趣的科研人员。 使用场景及目标:①评估不同操作条件下气化炉内炉渣的行为变化;②预测并优化气化炉的操作参数(如温度、氧煤比等),以防止炉渣堵塞;③为工业气化炉的设计操作提供理论支持技术指导。 其他说明:该模型的实现基于理论公式经验数据,为确保模型准确性,实际应用中需要根据具体气化炉的数据进行参数校准。模型还考虑了多个物理场的耦合,包括质量、动量能量守恒方程,能够模拟不同操作条件下的渣层演变。此外,提供了稳态求解器动态模拟工具,可用于扰动测试工业应用案例分析。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值