https://buaacoding.cn/contest-ng/index.html#/188/problems

其实这题挺简单的。

注意到答案的大小最多是22

二分，check长度是mid的不同子串有多少个，然后，如果不是1 << mid个，那么肯定是不行的。

想到了这个，我居然用了字符串hash去判，傻逼了，注意到是01串，相当于二进制。

#include <bits/stdc++.h>

#define IOS ios::sync_with_stdio(false)

using namespace std;

#define inf (0x3f3f3f3f)

typedef long long int LL;

const int maxn = 5e6 + ;

typedef unsigned long long int ULL;

char str[maxn];

int lenstr;

int cnt[ << ], DFN;

int po[maxn];

bool check(int mid) {

    DFN++;

    int val = ;

    for (int i = ; i <= mid; ++i) {

        val = val *  + str[i] - '';

    }

    cnt[val] = DFN;

    for (int i = mid + ; i <= lenstr; ++i) {

        val -= po[mid - ] * (str[i - mid] - '');

        val *= ;

        val += str[i] - '';

        cnt[val] = DFN;

    }

    int en = ( << mid) - ;

    for (int i = ; i <= en; ++i) {

        if (cnt[i] != DFN) {

            return true;

        }

    }

    return false;

}

void work() {

    scanf("%s", str + );

    lenstr = strlen(str + );

    int be = , en = min(, lenstr);

    while (be <= en) {

        int mid = (be + en) >> ;

        if (check(mid)) en = mid - ;

        else be = mid + ;

    }

    printf("%d\n", be);

}

int main() {

#ifdef local

    freopen("data.txt", "r", stdin);

//    freopen("data.txt", "w", stdout);

#endif

    po[] = ;

    for (int i = ; i <= ; ++i) po[i] = po[i - ] * ;

    int t;

    scanf("%d", &t);

    while (t--) work();

    return ;

}

很奇怪为什么能卡我dc3，我能算出长度是i的不同子串的个数，

思路就是，每一个sa[i]，有lenstr - sa[i] + 1个前缀，然后和前面的sa[i - 1] 会有重叠的地方，也就是height[i]个，然后减去即可。

区间减法打标记。但是，，TLE额。

why

可能多了点常数

/*

 Author: stupid_one

 Result: TLE    Submission_id: 506947

 Created at: Sun Dec 17 2017 10:45:12 GMT+0800 (CST)

 Problem_id: 959    Time: 1000    Memory: 0

*/

#include <bits/stdc++.h>

#define IOS ios::sync_with_stdio(false)

using namespace std;

#define inf (0x3f3f3f3f)

typedef long long int LL;

const int maxn = 9e6 + ;

typedef unsigned long long int ULL;

char str[maxn];

int lenstr;

int cnt[maxn];

const int N = maxn;

#define F(x) ((x)/3+((x)%3==1?0:tb))

#define G(x) ((x)<tb?(x)*3+1:((x)-tb)*3+2)

int r[maxn];

int wa[maxn],wb[maxn],wv[maxn],WS[maxn];

int sa[maxn];

int c0(int *r,int a,int b) {

    return r[a]==r[b]&&r[a+]==r[b+]&&r[a+]==r[b+];

}

int c12(int k,int *r,int a,int b) {

    if(k==) return r[a]<r[b]||r[a]==r[b]&&c12(,r,a+,b+);

    else return r[a]<r[b]||r[a]==r[b]&&wv[a+]<wv[b+];

}

void sort(int *r,int *a,int *b,int n,int m) {

    int i;

    for(i=; i<n; i++) wv[i]=r[a[i]];

    for(i=; i<m; i++) WS[i]=;

    for(i=; i<n; i++) WS[wv[i]]++;

    for(i=; i<m; i++) WS[i]+=WS[i-];

    for(i=n-; i>=; i--) b[--WS[wv[i]]]=a[i];

    return;

}

void dc3(int *r,int *sa,int n,int m) { //涵义与DA 相同

    int i,j,*rn=r+n,*san=sa+n,ta=,tb=(n+)/,tbc=,p;

    r[n]=r[n+]=;

    for(i=; i<n; i++) if(i%!=) wa[tbc++]=i;

    sort(r+,wa,wb,tbc,m);

    sort(r+,wb,wa,tbc,m);

    sort(r,wa,wb,tbc,m);

    for(p=,rn[F(wb[])]=,i=; i<tbc; i++)

        rn[F(wb[i])]=c0(r,wb[i-],wb[i])?p-:p++;

    if(p<tbc) dc3(rn,san,tbc,p);

    else for(i=; i<tbc; i++) san[rn[i]]=i;

    for(i=; i<tbc; i++) if(san[i]<tb) wb[ta++]=san[i]*;

    if(n%==) wb[ta++]=n-;

    sort(r,wb,wa,ta,m);

    for(i=; i<tbc; i++) wv[wb[i]=G(san[i])]=i;

    for(i=,j=,p=; i<ta && j<tbc; p++)

        sa[p]=c12(wb[j]%,r,wa[i],wb[j])?wa[i++]:wb[j++];

    for(; i<ta; p++) sa[p]=wa[i++];

    for(; j<tbc; p++) sa[p]=wb[j++];

    return;

}

int RANK[maxn], height[maxn];

void calheight(int *r,int *sa,int n) { // 此处N为实际长度

    int i,j,k=;        // height[]的合法范围为 1-N, 其中0是结尾加入的字符

    for(i=; i<=n; i++) RANK[sa[i]]=i; // 根据SA求RANK

    for(i=; i<n; height[RANK[i++]] = k ) // 定义：h[i] = height[ rank[i] ]

        for(k?k--:,j=sa[RANK[i]-]; r[i+k]==r[j+k]; k++); //根据 h[i] >= h[i-1]-1 来优化计算height过程

}

void work() {

    scanf("%s", str);

    lenstr = strlen(str);

    for (int i = ; i < lenstr; ++i) r[i] = str[i];

    r[lenstr] = ;

    dc3(r, sa, lenstr + , );

    calheight(r, sa, lenstr);

//    printf("%d\n", sa[1]);

    for (int i = ; i <= lenstr; ++i) {

        cnt[]++;

        cnt[lenstr - sa[i] + ]--;

        cnt[]--;

        cnt[height[i] + ]++;

    }

    for (int i = ; i < ; ++i) cnt[i] += cnt[i - ];

    for (int i = ; i <= ; ++i) {

        if (cnt[i] != ( << (i))) {

            printf("%d\n", i);

            for (int k = ; k < ; ++k) cnt[k] = ;

            return;

        }

    }

}

int main() {

#ifdef local

    freopen("data.txt", "r", stdin);

//    freopen("data.txt", "w", stdout);

#endif

    int t;

    scanf("%d", &t);

    while (t--) work();

    return ;

}
/*

 Author: stupid_one

 Result: TLE    Submission_id: 506775

 Created at: Sun Dec 17 2017 01:21:28 GMT+0800 (CST)

 Problem_id: 959    Time: 1000    Memory: 0

*/

#include <bits/stdc++.h>

#define IOS ios::sync_with_stdio(false)

using namespace std;

#define inf (0x3f3f3f3f)

typedef long long int LL;

const int maxn = 5e6 + ;

typedef unsigned long long int ULL;

ULL po[maxn];

char str[maxn];

int lenstr;

ULL hs[maxn];

const int seed = ;

const int mod = ;

struct MP {

    int first[mod + ], num;

    ULL eval[maxn];

    int nxt[maxn];

    void init() {

        memset(first, -, sizeof first);

        num = ;

    }

    bool addEdge(ULL val) {

        int u = val % mod;

        for (int i = first[u]; ~i; i = nxt[i]) {

            if (val == eval[i]) return false;

        }

        eval[num] = val;

        nxt[num] = first[u];

        first[u] = num++;

        return true;

    }

} h;

bool check(int val) {

    int t =  << val;

    if (lenstr - val +  < t) return false;

    h.init();

    int ans = ;

    for (int i = val; i <= lenstr; ++i) {

        ans += h.addEdge(hs[i] - po[val] * hs[i - val]);

        if (ans == t) return true;

        if (ans + (lenstr - i) < t) return false;

    }

    return false;

}

int sa[maxn], x[maxn], y[maxn], book[maxn]; //book[]大小起码是lenstr，book[rank[]]

bool cmp(int r[], int a, int b, int len) { //这个必须是int r[]，

    return r[a] == r[b] && r[a + len] == r[b + len];

}

void da(char str[], int sa[], int lenstr, int mx) {

    int *fir = x, *sec = y, *ToChange;

    for (int i = ; i <= mx; ++i) book[i] = ; //清0

    for (int i = ; i <= lenstr; ++i) {

        fir[i] = str[i];  //开始的rank数组，只保留相对大小即可，开始就是str[]

        book[str[i]]++; //统计不同字母的个数

    }

    for (int i = ; i <= mx; ++i) book[i] += book[i - ]; //统计 <= 这个字母的有多少个元素

    for (int i = lenstr; i >= ; --i) sa[book[fir[i]]--] = i;

// <=str[i]这个字母的有x个，那么，排第x的就应该是这个i的位置了。

    //倒过来排序，是为了确保相同字符的时候，前面的就先在前面出现。

    //p是第二个关键字0的个数

    for (int j = , p = ; p <= lenstr; j <<= , mx = p) {  //字符串长度为j的比较

        //现在求第二个关键字，然后合并（合并的时候按第一关键字优先合并）

        p = ;

        for (int i = lenstr - j + ; i <= lenstr; ++i) sec[++p] = i;

//上面的位置，再跳j格就是越界了的，所以第二关键字是0，排在前面

        for (int i = ; i <= lenstr; ++i)

            if (sa[i] > j)  //如果排名第i的起始位置在长度j之后

                sec[++p] = sa[i] - j;

//减去这个长度j，表明第sa[i] - j这个位置的第二个是从sa[i]处拿的，排名靠前也//正常，因为sa[i]排名是递增的

        //sec[]保存的是下标，现在对第一个关键字排序

        for (int i = ; i <= mx; ++i) book[i] = ; //清0

        for (int i = ; i <= lenstr; ++i) book[fir[sec[i]]]++;

        for (int i = ; i <= mx; ++i) book[i] += book[i - ];

        for (int i = lenstr; i >= ; --i) sa[book[fir[sec[i]]]--] = sec[i];

//因为sec[i]才是对应str[]的下标

//现在要把第二关键字的结果，合并到第一关键字那里。同时我需要用到第一关键//字保存的记录，所以用指针交换的方式达到快速交换数组中的值

        ToChange = fir, fir = sec, sec = ToChange;

        fir[sa[]] = ; //固定的是0 因为sa[1]固定是lenstr那个0

        p = ;

        for (int i = ; i <= lenstr; ++i)  //fir是当前的rank值，sec是前一次的rank值

            fir[sa[i]] = cmp(sec, sa[i - ], sa[i], j) ? p -  : p++;

    }

    return ;

}

int height[maxn], RANK[maxn];

void calcHight(char str[], int sa[], int lenstr) {

    for (int i = ; i <= lenstr; ++i) RANK[sa[i]] = i; //O(n)处理出rank[]

    int k = ;

    for (int i = ; i <= lenstr - ; ++i) {

//最后一位不用算，最后一位排名一定是1，然后sa[0]就尴尬了

        k -= k > ;

        int j = sa[RANK[i] - ]; //排名在i前一位的那个串，相似度最高

        while (str[j + k] == str[i + k]) ++k;

        height[RANK[i]] = k;

    }

    return ;

}

int cnt[maxn];

void work() {

    scanf("%s", str + );

    lenstr = strlen(str + );

    for (int i = ; i <= lenstr; ++i) cnt[i] = ;

    str[lenstr + ] = '$';

    str[lenstr + ] = ;

    da(str, sa, lenstr + , );

    calcHight(str, sa, lenstr + );

    for (int i = ; i <= lenstr + ; ++i) {

        cnt[]++;

        cnt[lenstr - sa[i] + ]--;

        cnt[]--;

        cnt[height[i] + ]++;

    }

    for (int i = ; i <= lenstr; ++i) cnt[i] += cnt[i - ];

    for (int i = ; i <= lenstr; ++i) {

        if (cnt[i] != ( << i)) {

            printf("%d\n", i);

            return;

        }

    }

//    printf("\n");

//    for (int i = 1; str[i]; ++i) hs[i] = hs[i - 1] * seed + str[i];

//    for (int i = 1; i <= 22; ++i) {

//        if (!check(i)) {

//            printf("%d\n", i);

//            return ;

//        }

//    }

}

int main() {

#ifdef local

    freopen("data.txt", "r", stdin);

//    freopen("data.txt", "w", stdout);

#endif

    po[] = ;

    for (int i = ; i <= maxn - ; ++i) po[i] = po[i - ] * seed;

    int t;

    scanf("%d", &t);

    while (t--) work();

    return ;

}

da

B P5 第十三届北航程序设计竞赛预赛的更多相关文章

  1. 第十三届北航程序设计竞赛决赛网络同步赛 B题 校赛签到（建树 + 打标记）

    题目链接  校赛签到 对每个操作之间建立关系. 比较正常的是前$3$种操作,若第$i$个操作属于前$3$种,那么就从操作$i-1$向$i$连一条有向边. 比较特殊的是第$4$种操作,若第$i$个操作属 ...

  2. 《内蒙古自治区第十三届大学生程序设计竞赛试题_H 公孙玉龙》

    这个题有点小坑,最坑的地方就是  输入的b 变量  有可能 是  负数  !   负数 !  负数!    对 ,你没有看错,就是负数,坑死我了, 一直都是  content.charAt(0) 判断 ...

  3. 2020年“感恩杯”台州学院第十三届大学生程序设计竞赛D、H、I题解（后续补充）

    D题:小z与他的袜子 描述 小z每天会穿一双新袜子. 开始他的衣柜里有n双袜子,袜子会从1-n进行编号.每天早上他都会从衣柜里拿编号最小的袜子来穿.每天晚上他会把今天穿的袜子扔进篮子里,如果篮子里有n ...

  4. [题解]第十一届北航程序设计竞赛预赛——L.偶回文串

    题目描述 长度为偶数的回文串被称为偶回文串.如果一个字符串重新排序之后能够成为一个偶回文串,则称为可回文的. 给一个字符串,求可回文的子串个数.字符串只含小写字母,单个字符串长度不超过10^5,所有数 ...

  5. [题解]第十一届北航程序设计竞赛预赛——I.神奇宝贝大师

    题目描述 一张n*m的地图,每个格子里面有一定数量的神奇宝贝,求一个最优位置,使得所有神奇宝贝到该位置的曼哈顿距离最小. 一共有T组数据,每组数据包含两行,第一行是n和m(1<=n,m<= ...

  6. [题解]第十一届北航程序设计竞赛预赛——H.高中数学题

    题目描述 解题思路 可以求得通项公式:an = 2n + 1,所以问题就变成等差数列求异或和,这个具体为什么对我还不能很好地解释清楚,先挖坑吧. 附:c++代码 1 #include <iost ...

  7. [题解]第十一届北航程序设计竞赛预赛——F.序列

    题目描述 (1,--,n)的一个排列S,定义其对应的权值F[S]为:将S划分为若干段连续子序列,每个子序列都是上升序列,F[S]的值等于能划分出的最小段数. 求n的全排列的F[S]的和,答案mod(1 ...

  8. [题解]第十一届北航程序设计竞赛预赛——D.最大公约数

    题目描述 给一个长度为n(1<=n<=100000)的正整数列,分成尽量多的非空段,使得每一段的最大公约数相等.一个数的最大公约数是它本身. 解题思路 要求每一段子列的gcd相等,不妨设为 ...

  9. [题解]第十一届北航程序设计竞赛预赛——A.模式

    题目描述 输入一个学号,判断是计算机系or软件学院or其他院系. 解题思路 水题,直接判断or除以10000都可以.不废话,直接上代码. 1 #include <iostream> 2 # ...

随机推荐

  1. WordCount编码测试

    Github项目地址:https://github.com/LantyrLYL/WordCount PSP表格: PSP2.1 PSP阶段 预估耗时 (分钟) 实际耗时 (分钟) Planning 计 ...

  2. LeetCode第136题：只出现一次的数字

    问题描述 给定一个非空整数数组,除了某个元素只出现一次以外,其余每个元素均出现两次.找出那个只出现了一次的元素. 说明: 你的算法应该具有线性时间复杂度. 你可以不使用额外空间来实现吗? 示例 1: ...

  3. Android之ContextMenu的使用方法以及与OptionMenu的区别(转)

    >> ContextMenu是android的context menu上下文菜单,选择某项VIEW后长按menu键,就会显示出来.比如EditeText就可以通过长按来弹出拥有“cut”, ...

  4. register的功能

    register:这个关键字请求编译器尽可能 的将变量存在CPU内部寄存器中,而不是通过内存寻址访问,以提高效率.这里注意是尽可能,不是绝对.你想想,一个CPU拥有 的寄存器也就那么几个或几十个,你要 ...

  5. PHP 中for循环的一个小小改进

    注意 : 1, $size这个值可以放在for循环中的第一个 ';' 前声明, 因为这个声明只会执行一次; 2, 第二个 ';' 中的内容, 会重复运行, 所以$i < $size 这个判断会每 ...

  6. c#委托与事件、消息、WndProc用法(转)

    c#委托与事件 心得 c#用委托来实现事件通知机制.委托相当与c++函数指针.整个过程涉及一个呼叫者,一个被呼叫者,还有就是这个委托. - 实现步骤有以下几步: 1. 申明委托, 2.定义呼叫者和调用 ...

  7. Copy拷贝

    前言 copy:需要先实现 NSCopying 协议,创建的是不可变副本. mutableCopy:需要实现 NSMutableCopying 协议,创建的是可变副本. 浅拷贝:指针拷贝,源对象和副本 ...

  8. Python学习过程（二）

    条件判断和循环 条件判断 age = 20 if age >= 18: print 'your age is',age print 'adult' elif age >=6 : print ...

  9. git 本地分支与远程分支相关操作记录

    1.远程分支中有新增分支,但自己的本地分支没有对应同步 git checkout -b [remote-branch-name] origin/[remote-branch-name] 2. 查看本地 ...

  10. Geometry - DbGeometry的使用说明一

    说明:工作中使用过但是没有详细的研究过,使用c#语言进行逻辑处理.分享出来希望各位发表意见 geometry是arcgis的空间对象 dbgeometry是微软的空间对象 geometry对象转换为d ...