惯性聚合 高效追踪和阅读你感兴趣的博客、新闻、科技资讯
阅读原文 在惯性聚合中打开

推荐订阅源

阮一峰的网络日志
阮一峰的网络日志
D
Darknet – Hacking Tools, Hacker News & Cyber Security
S
Schneier on Security
The Last Watchdog
The Last Watchdog
Cyberwarzone
Cyberwarzone
S
Securelist
Threat Intelligence Blog | Flashpoint
Threat Intelligence Blog | Flashpoint
C
Cyber Attacks, Cyber Crime and Cyber Security
L
Lohrmann on Cybersecurity
钛媒体:引领未来商业与生活新知
钛媒体:引领未来商业与生活新知
博客园 - 司徒正美
The Cloudflare Blog
V
V2EX
博客园_首页
博客园 - 聂微东
Vercel News
Vercel News
人人都是产品经理
人人都是产品经理
G
GRAHAM CLULEY
T
Tenable Blog
Last Week in AI
Last Week in AI
Y
Y Combinator Blog
L
LINUX DO - 最新话题
cs.CL updates on arXiv.org
cs.CL updates on arXiv.org
SecWiki News
SecWiki News
博客园 - 三生石上(FineUI控件)
S
Secure Thoughts
N
News | PayPal Newsroom
T
The Blog of Author Tim Ferriss
The GitHub Blog
The GitHub Blog
T
Troy Hunt's Blog
博客园 - 【当耐特】
Forbes - Security
Forbes - Security
H
Hacker News: Front Page
A
About on SuperTechFans
B
Blog RSS Feed
Engineering at Meta
Engineering at Meta
MongoDB | Blog
MongoDB | Blog
CTFtime.org: upcoming CTF events
CTFtime.org: upcoming CTF events
freeCodeCamp Programming Tutorials: Python, JavaScript, Git & More
罗磊的独立博客
D
DataBreaches.Net
P
Privacy & Cybersecurity Law Blog
Schneier on Security
Schneier on Security
Application and Cybersecurity Blog
Application and Cybersecurity Blog
Google DeepMind News
Google DeepMind News
奇客Solidot–传递最新科技情报
奇客Solidot–传递最新科技情报
OSCHINA 社区最新新闻
OSCHINA 社区最新新闻
Jina AI
Jina AI
D
Docker
P
Proofpoint News Feed

博客园 - xmx

Amazon 2面杯具 N皇后回溯 "中航文化杯" 2007 ACM/ICPC 国际大学生程序设计竞赛亚洲区域赛(南京) 一个用来练dfs的简单迷宫问题 pku 1662 还是找规律的 pku 1806 Manhattan 2025(找规律) 今天西华的比赛,啥都不说啦,相当的nice!~~~ 今天北京赛区的比赛 pku 1505 copying books(DP) 最近看的一些东西 The 2007 ACM Asia Programming Contest Changchun Site Internet Preliminary Contest nice 位运算果真是好东西,今天算是学到点啦^_^ FOJ月赛-2007年9月 pku 1850 前面一直没注意到某个不规范的情况,导致结果一直比标准的大...调了好久... The 2007 ACM Asia Programming Contest - Nanjing Preliminary pku 3219 人家居然用几十B就过了,肯定有超强的规律,可是我自己找了个,挂了...只能老实算... 一道双向dp,差点超时^_^||| dp pku 1050 N和素数P,求杨辉三角第N行中能被P整除的数的个数
终于有算最长重复子串(数)的后缀数组啦,nice
xmx · 2007-09-30 · via 博客园 - xmx

#include <cstdio>
#include <string>
#include <algorithm>
using namespace std;

const int MAX = 25000;

  1#include <cstdio>
  2#include <string>
  3#include <algorithm>
  4using namespace std;
  5
  6const int MAX = 25000;
  7int num[MAX];
  8int mem[3][MAX], c[MAX];
  9int * SA, * Rank, * nSA, * nRank;
 10int n,kk,nmax;
 11
 12void make_sa()
 13{
 14    int i,j;
 15    SA = mem[0]; Rank = mem[1]; nSA = mem[2];
 16    memset(c,0,sizeof(c));
 17    for(i=0;i<n;i++{
 18        c[ num[i] +1++;
 19    }

 20    for(i=1;i<=nmax+1;i++{
 21        c[i] += c[i-1];
 22    }

 23    for(i=0;i<n;i++{
 24        SA[ --c[ num[i] +1] ] = i;
 25    }

 26    Rank[ SA[0] ] = 0;
 27    for(i=1;i<n;i++{
 28        Rank[SA[i]] = Rank[SA[i-1]];
 29        if(num[SA[i]] != num[SA[i-1]]) {
 30            Rank[SA[i]] ++;
 31        }

 32    }

 33
 34    int k;
 35    for(k=1;k<&& Rank[ SA[n-1] ] < n-1;k*=2{
 36        memset(c,0,sizeof(c));
 37        for(i=0;i<n;i++{
 38            c[ Rank[SA[i]] ] ++;
 39        }

 40        for(i=1;i<n;i++{
 41            c[i] += c[i-1];
 42        }

 43        for(i=n-1;i>=0;i--{
 44            if(SA[i] >= k) {
 45                nSA[ -- c[ Rank[SA[i]-k] ] ] = SA[i] - k;
 46            }

 47        }

 48        for(i=n-k;i<n;i++{
 49            nSA[ -- c[Rank[i]] ] = i;
 50        }

 51        nRank = SA;
 52        nRank[ nSA[0] ] = 0;
 53        for(i=1;i<n;i++{
 54            nRank[ nSA[i] ] = nRank[ nSA[i-1] ];
 55            if(Rank[nSA[i]] != Rank[nSA[i-1]] || Rank[nSA[i]+k] != Rank[nSA[i-1]+k]) {
 56                nRank[ nSA[i] ] ++;
 57            }

 58        }

 59        SA = nSA;
 60        nSA = Rank;
 61        Rank = nRank;
 62    }

 63}

 64
 65int h[MAX];
 66int get_lcp()
 67{
 68    int i,j,k;
 69    
 70    for(i=0,k=0;i<n;i++{
 71        if(Rank[i] == n-1{
 72            h[ Rank[i] ] = k = 0;
 73        }

 74        else {
 75            if(k > 0{
 76                k --;
 77            }

 78            j = SA[Rank[i] +1];
 79            while(num[i+k] == num[j+k]) {
 80                k ++;
 81            }

 82            h[ Rank[i] ] = k;
 83        }

 84    }

 85    int ret = 0;
 86    kk --;
 87    for(i=0;i<n-kk;i++{
 88        int t = INT_MAX;
 89        for(j=0;j<kk;j++{
 90            t = min(t, h[i+j]);
 91        }

 92        ret = max(t, ret);
 93    }

 94    return ret;
 95}

 96
 97int main()
 98{
 99    int i,j;
100    while(scanf("%d %d"&n,&kk)==2{
101        nmax = -1;
102        for(i=0;i<n;i++{
103            scanf("%d", num+i);
104            nmax = max(nmax, num[i]);
105        }

106        num[n ++= -1;
107        make_sa();
108        printf("%d\n", get_lcp());
109    }

110}

111