B P5 第十三届北航程序设计竞赛预赛
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 第十三届北航程序设计竞赛预赛的更多相关文章
- 第十三届北航程序设计竞赛决赛网络同步赛 B题 校赛签到(建树 + 打标记)
题目链接 校赛签到 对每个操作之间建立关系. 比较正常的是前$3$种操作,若第$i$个操作属于前$3$种,那么就从操作$i-1$向$i$连一条有向边. 比较特殊的是第$4$种操作,若第$i$个操作属 ...
- 《内蒙古自治区第十三届大学生程序设计竞赛试题_H 公孙玉龙》
这个题有点小坑,最坑的地方就是 输入的b 变量 有可能 是 负数 ! 负数 ! 负数! 对 ,你没有看错,就是负数,坑死我了, 一直都是 content.charAt(0) 判断 ...
- 2020年“感恩杯”台州学院第十三届大学生程序设计竞赛D、H、I题解(后续补充)
D题:小z与他的袜子 描述 小z每天会穿一双新袜子. 开始他的衣柜里有n双袜子,袜子会从1-n进行编号.每天早上他都会从衣柜里拿编号最小的袜子来穿.每天晚上他会把今天穿的袜子扔进篮子里,如果篮子里有n ...
- [题解]第十一届北航程序设计竞赛预赛——L.偶回文串
题目描述 长度为偶数的回文串被称为偶回文串.如果一个字符串重新排序之后能够成为一个偶回文串,则称为可回文的. 给一个字符串,求可回文的子串个数.字符串只含小写字母,单个字符串长度不超过10^5,所有数 ...
- [题解]第十一届北航程序设计竞赛预赛——I.神奇宝贝大师
题目描述 一张n*m的地图,每个格子里面有一定数量的神奇宝贝,求一个最优位置,使得所有神奇宝贝到该位置的曼哈顿距离最小. 一共有T组数据,每组数据包含两行,第一行是n和m(1<=n,m<= ...
- [题解]第十一届北航程序设计竞赛预赛——H.高中数学题
题目描述 解题思路 可以求得通项公式:an = 2n + 1,所以问题就变成等差数列求异或和,这个具体为什么对我还不能很好地解释清楚,先挖坑吧. 附:c++代码 1 #include <iost ...
- [题解]第十一届北航程序设计竞赛预赛——F.序列
题目描述 (1,--,n)的一个排列S,定义其对应的权值F[S]为:将S划分为若干段连续子序列,每个子序列都是上升序列,F[S]的值等于能划分出的最小段数. 求n的全排列的F[S]的和,答案mod(1 ...
- [题解]第十一届北航程序设计竞赛预赛——D.最大公约数
题目描述 给一个长度为n(1<=n<=100000)的正整数列,分成尽量多的非空段,使得每一段的最大公约数相等.一个数的最大公约数是它本身. 解题思路 要求每一段子列的gcd相等,不妨设为 ...
- [题解]第十一届北航程序设计竞赛预赛——A.模式
题目描述 输入一个学号,判断是计算机系or软件学院or其他院系. 解题思路 水题,直接判断or除以10000都可以.不废话,直接上代码. 1 #include <iostream> 2 # ...
随机推荐
- 抓包工具 Fiddler 使用:弱网络环境模拟限速测试流程
转自:http://www.51testing.com/html/80/n-3726980.html 抓包工具 Fiddler 使用:弱网络环境模拟限速测试流程 发表于:2018-6-06 11: ...
- Java代码执行机制
一.java源码编译机制 一个java文件编译为.class文件的流程图. Parse an Enter阶段: Parse:将代码字符串转换为token序列(com.sun.tools.javac.p ...
- WPARAM和LPARAM的含义
lParam 和 wParam 是宏定义,一般在消息函数中带这两个类型的参数,通常用来存储窗口消息的参数. LRESULT CALLBACK WindowProc(HWND hwnd, UINT uM ...
- Category 分类
1.Category 1)分类/类别(category): 允许以模块的方式向现有类定义添加新的方法(默认不能添加实例变量).用以扩展自己或他人以前实现的类,使它适合自己的需要. 分类的名称括在类名之 ...
- Total Commander的初次体验
从汉化新世纪下载到最新的TC张学思版后,运行文件只需依照其提示就可以完成该软件的安装.作为新手初次运行体验了以下功能: 一.目录跳转 1. 初次启动TC软件界面截图: 2. 按下Ctrl+d后,直接再 ...
- Spring core注解
1.@Autowired Autowired是用在JavaBean中的注解,通过byType形式,用来给指定的字段或方法注入所需的外部资源 Autowired注解来指定自动装配,可以修饰setter方 ...
- VIRT(虚拟内存)RES(常驻内存)和SHR(共享内存)
VIRT: 1.进程“需要的”虚拟内存大小,包括进程使用的库.代码.数据,以及malloc.new分配的堆空间和分配的栈空间等: 2.假如进程新申请10MB的内存,但实际只使用了1MB,那么它会增长1 ...
- jzoj4916. 【GDOI2017模拟12.9】完全背包问题 (背包+最短路)
题面 题解 考场上蠢了--这么简单的东西都想不到-- 首先排序加去重. 先来考虑一下,形如 \[a_1x_1+a_2x_2+...a_nx_n=w,a_1<a_2<...<a_n,x ...
- 各平台安装和使用Docker的差异
在Mac上运行Docker 使用Docker for Mac 系统要求 开启虚拟化的硬件支持,可以通过命令来检查:sysctl kern.hv_support OS X 10.10.3 或者更高版本 ...
- node.js 设置 淘宝 镜像
npm config set registry "https://registry.npm.taobao.org" npm info underscore (如果上面配置正确这个命 ...