HEOI2016 题解

  • Q:为什么要在sdoi前做去年的heoi题
  • A:我省选药丸

http://cogs.pro/cogs/problem/index.php?key=heoi2016

D1T1 树

树剖直接做

或者操作逆序,并查集维护自己和祖先中第一个打标记的点,注意1号点本身有标记,mark[1]++

代码在最后

D1T2排序

计数排序拿80...

二分答案,\(\ge\)设为1,\(<\)设为0,排序就可以用线段树模拟了. 鏼鏼鏼!

代码还没写

D1T3序列

md一开始看错题了

写出\(O(n^2)\)的DP,发现裸三维偏序啊,\(j<i,\ max_j < a_i,\ a_j < min_i\)

注意最后调用cdq(mid+1, r)前按id排一下序

cdq真是太棒了!

代码在最后

D2T1游戏

二分图匹配...

硬石头把一行拆成了好几个点...

代码在最后

D2T2求和

以前做了

D2T3字符串

题意:多次询问\(S[a,b]\)的所有子串与\(S[c,d]\)的LCP

40分\(O(n^2)\)暴力 每次重构sam

一开始想的是,反转后建sam,建出Parent Tree,往上找第一个Right集合中有[a,b]的节点。问题在于限制了在[a,b]范围内,不一定是越深的祖先越好因为到b就截断了

考虑二分lcp长度,那么变成了判断Right有没有[a+len-1, b]范围内的点(注意我反转了字符串)

这样就可以倍增了,判断可以对dfs序建主席树权值是后缀

但这样是三个log

其实等价于我们先倍增到val>=len的最浅的点,看看Right集合里有没有,这样就是两个log了

代码

  1. #include <iostream>
  2. #include <cstdio>
  3. #include <cstring>
  4. #include <algorithm>
  5. #include <cmath>
  6. using namespace std;
  7. typedef long long ll;
  8. const int N=1e5+5;
  9. inline int read(){
  10. char c=getchar();int x=0,f=1;
  11. while(c<'0'||c>'9'){if(c=='-')f=-1;c=getchar();}
  12. while(c>='0'&&c<='9'){x=x*10+c-'0';c=getchar();}
  13. return x*f;
  14. }
  15. int n, Q, fa[N], u, v, f[N], ans[N], m, mark[N];
  16. int find(int x) {return x==fa[x] ? x : fa[x]=find(fa[x]);}
  17. struct meow{char s[5]; int x, id;} q[N];
  18. struct edge{int v, ne;} e[N];
  19. int cnt=1, h[N];
  20. inline void ins(int u, int v) {e[++cnt]=(edge){v, h[u]}; h[u]=cnt;}
  21. namespace graph {
  22. int q[N], head, tail;
  23. void bfs() {
  24. head=tail=1;
  25. q[tail++]=1;
  26. while(head != tail) {
  27. int u=q[head++];
  28. if(!fa[u]) fa[u] = find(f[u]);
  29. for(int i=h[u];i;i=e[i].ne) q[tail++] = e[i].v;
  30. }
  31. }
  32. }
  33. int main() {
  34. //freopen("in","r",stdin);
  35. freopen("heoi2016_tree.in","r",stdin);
  36. freopen("heoi2016_tree.out","w",stdout);
  37. n=read(); Q=read();
  38. fa[1]=1; mark[1]++;
  39. for(int i=1; i<n; i++) u=read(), v=read(), f[v]=u, ins(u, v);
  40. for(int i=1; i<=Q; i++) {
  41. scanf("%s", q[i].s), q[i].x=read();
  42. if(q[i].s[0] == 'C') fa[ q[i].x ] = q[i].x, mark[ q[i].x ]++;
  43. else q[i].id=++m;
  44. }
  45. graph::bfs();
  46. for(int i=Q; i>=1; i--) {
  47. if(q[i].s[0] == 'C') {
  48. mark[ q[i].x ]--;
  49. if(!mark[ q[i].x ]) fa[ q[i].x ] = find(f[q[i].x]);
  50. }else ans[q[i].id] = find(q[i].x);
  51. }
  52. for(int i=1; i<=m; i++) printf("%d\n", ans[i]);
  53. }
  1. #include <iostream>
  2. #include <cstdio>
  3. #include <cstring>
  4. #include <algorithm>
  5. #include <cmath>
  6. using namespace std;
  7. typedef long long ll;
  8. const int N=1e5+5;
  9. inline int read(){
  10. char c=getchar();int x=0,f=1;
  11. while(c<'0'||c>'9'){if(c=='-')f=-1;c=getchar();}
  12. while(c>='0'&&c<='9'){x=x*10+c-'0';c=getchar();}
  13. return x*f;
  14. }
  15. int n, Q, x, y, f[N];
  16. struct meow{int a, b, c, id;} a[N]; //max --> b, min --> c;
  17. inline bool cmpb(const meow &a, const meow &b) {return a.b < b.b;}
  18. inline bool cmpa(const meow &a, const meow &b) {return a.a < b.a;}
  19. inline bool cmpid(const meow &a, const meow &b) {return a.id < b.id;}
  20. namespace bit {
  21. int n, c[N], t[N], T;
  22. inline void ini() {T++; }
  23. inline void upd(int p, int v) { //printf("upd %d %d\n",p,v);
  24. for(; p<=n; p+=p&-p) {
  25. if(t[p] != T) t[p]=T, c[p]=v;
  26. else if(c[p] < v) c[p]=v;
  27. }
  28. }
  29. inline int max(int p) { //printf("max %d\n",p);
  30. int ans=0;
  31. for(; p; p-=p&-p) if(t[p] == T && c[p] > ans) ans = c[p];
  32. return ans;
  33. }
  34. }
  35. void cdq(int l, int r) { //printf("\n-------cdq %d %d\n",l,r);
  36. if(l==r) return;
  37. int mid=(l+r)>>1;
  38. cdq(l, mid);
  39. sort(a+l, a+mid+1, cmpb); sort(a+mid+1, a+r+1, cmpa);
  40. bit::ini();
  41. int p1=l, p2=mid+1;
  42. while(p1 <= mid || p2 <= r) {
  43. if(p2 > r || (p1 <= mid && a[p1].b <= a[p2].a) ) bit::upd(a[p1].a, f[a[p1].id]), p1++;
  44. else f[a[p2].id] = max(f[a[p2].id], bit::max(a[p2].c)+1), p2++; //printf("p2 %d %d %d\n", p2, a[p2].id, f[a[p2].id]),p2++;
  45. }
  46. sort(a+mid+1, a+r+1, cmpid);
  47. cdq(mid+1, r);
  48. }
  49. int main() {
  50. //freopen("in","r",stdin);
  51. freopen("heoi2016_seq.in","r",stdin);
  52. freopen("heoi2016_seq.out","w",stdout);
  53. n=read(); Q=read();
  54. for(int i=1; i<=n; i++) x=read(), a[i]=(meow){x, x, x, i}, f[i]=1, bit::n = max(bit::n, x);
  55. for(int i=1; i<=Q; i++) x=read(), y=read(), a[x].b = max(a[x].b, y), a[x].c = min(a[x].c, y), bit::n = max(bit::n, y);
  56. cdq(1, n);
  57. int ans=0;
  58. for(int i=1; i<=n; i++) ans = max(ans, f[i]);// printf("%d ",f[i]); puts(" f");
  59. printf("%d", ans);
  60. }
  1. #include <iostream>
  2. #include <cstdio>
  3. #include <cstring>
  4. #include <algorithm>
  5. #include <cmath>
  6. using namespace std;
  7. typedef long long ll;
  8. const int N=5100, M=1e5+5, INF=1e9+5;
  9. inline int read(){
  10. char c=getchar();int x=0,f=1;
  11. while(c<'0'||c>'9'){if(c=='-')f=-1;c=getchar();}
  12. while(c>='0'&&c<='9'){x=x*10+c-'0';c=getchar();}
  13. return x*f;
  14. }
  15. int n, m, s, t, X, Y, idx[55][55], idy[55][55];
  16. char str[55][55];
  17. namespace mf {
  18. struct edge{int v, ne, c, f;}e[M];
  19. int cnt=1, h[N];
  20. inline void ins(int u, int v, int c) {
  21. e[++cnt]=(edge){v, h[u], c, 0}; h[u]=cnt;
  22. e[++cnt]=(edge){u, h[v], 0, 0}; h[v]=cnt;
  23. }
  24. int q[N], head, tail, d[N], vis[N];
  25. bool bfs() {
  26. memset(vis, 0, sizeof(vis));
  27. head=tail=1;
  28. q[tail++]=s; d[s]=0; vis[s]=1;
  29. while(head != tail) {
  30. int u=q[head++];
  31. for(int i=h[u];i;i=e[i].ne)
  32. if(!vis[e[i].v] && e[i].c>e[i].f) {
  33. vis[e[i].v]=1;
  34. d[e[i].v] = d[u]+1;
  35. q[tail++]=e[i].v;
  36. if(e[i].v == t) return true;
  37. }
  38. }
  39. return false;
  40. }
  41. int cur[N];
  42. int dfs(int u, int a) {
  43. if(u==t || a==0) return a;
  44. int flow=0, f;
  45. for(int &i=cur[u];i;i=e[i].ne)
  46. if(d[e[i].v] == d[u]+1 && (f = dfs(e[i].v, min(a, e[i].c-e[i].f)) )>0) {
  47. flow += f;
  48. e[i].f += f;
  49. e[i^1].f -= f;
  50. a -= f;
  51. if(a==0) break;
  52. }
  53. if(a) d[u]=-1;
  54. return flow;
  55. }
  56. int dinic() {
  57. int flow=0;
  58. while(bfs()) {
  59. for(int i=s; i<=t; i++) cur[i]=h[i];
  60. flow += dfs(s, INF);
  61. }
  62. return flow;
  63. }
  64. int g[N][N];
  65. void build() {
  66. for(int i=1; i<=n; i++) {
  67. X++;
  68. for(int j=1; j<=m; j++) if(str[i][j] == '#') X++;
  69. else idx[i][j] = X;
  70. }
  71. Y=X;
  72. for(int j=1; j<=m; j++) {
  73. Y++;
  74. for(int i=1; i<=n; i++) if(str[i][j] == '#') Y++;
  75. else idy[i][j] = Y;
  76. }
  77. s=0; t=Y+1;
  78. //for(int i=1; i<=n; i++) for(int j=1; j<=m; j++) printf("idx %d %d %d\n",i,j,idx[i][j]);
  79. //for(int j=1; j<=m; j++) for(int i=1; i<=n; i++) printf("idy %d %d %d\n",i,j,idy[i][j]);
  80. int has[N];
  81. for(int i=1; i<=n; i++)
  82. for(int j=1; j<=m; j++) if(str[i][j] == '*')
  83. ins(idx[i][j], idy[i][j], 1), has[idx[i][j]]=1, has[idy[i][j]]=1;
  84. for(int i=1; i<=X; i++) if(has[i]) ins(s, i, 1);
  85. for(int i=X+1; i<=Y; i++) if(has[i]) ins(i, t, 1);
  86. }
  87. }
  88. int main() {
  89. //freopen("in","r",stdin);
  90. //freopen("heoi2016_game.in","r",stdin);
  91. //freopen("heoi2016_game.out","w",stdout);
  92. n=read(); m=read();
  93. for(int i=1; i<=n; i++) scanf("%s",str[i]+1);
  94. mf::build();
  95. int ans = mf::dinic();
  96. printf("%d", ans);
  97. }
  1. #include <iostream>
  2. #include <cstdio>
  3. #include <cstring>
  4. #include <algorithm>
  5. #include <cmath>
  6. using namespace std;
  7. typedef long long ll;
  8. #define lc(x) t[x].l
  9. #define rc(x) t[x].r
  10. const int N=2e5+5;
  11. inline int read(){
  12. char c=getchar();int x=0,f=1;
  13. while(c<'0'||c>'9'){if(c=='-')f=-1;c=getchar();}
  14. while(c>='0'&&c<='9'){x=x*10+c-'0';c=getchar();}
  15. return x*f;
  16. }
  17. int n, Q, pos[N], a, b, c, d;
  18. char s[N];
  19. namespace sam {
  20. struct meow{int ch[26], par, val;}t[N];
  21. int sz=1, last=1, root=1;
  22. int ri[N], suf[N];
  23. inline void extend(int c, int id) {
  24. int p=last, np=++sz; t[np].val = t[p].val+1;
  25. pos[id] = np; suf[np] = id; //printf("id %d %d\n",id,np);
  26. for(; p && !t[p].ch[c]; p=t[p].par) t[p].ch[c]=np;
  27. if(!p) t[np].par = root;
  28. else {
  29. int q = t[p].ch[c];
  30. if(t[q].val == t[p].val+1) t[np].par = q;
  31. else {
  32. int nq=++sz; t[nq]=t[q]; t[nq].val = t[p].val+1;
  33. t[q].par = t[np].par = nq;
  34. for(; p && t[p].ch[c]==q; p=t[p].par) t[p].ch[c]=nq;
  35. }
  36. }
  37. last=np;
  38. }
  39. } using sam::suf; using sam::extend;
  40. namespace tr {
  41. struct edge{int v, ne;} e[N];
  42. int cnt=1, h[N];
  43. inline void ins(int u, int v) { //printf("ins %d %d\n",u,v);
  44. e[++cnt]=(edge){v, h[u]}; h[u]=cnt;}
  45. struct meow{int l, r, size;} t[N*20];
  46. int sz=0, root[N];
  47. inline void insert(int &x, int l, int r, int p) {
  48. t[++sz]=t[x]; x=sz;
  49. t[x].size++;
  50. if(l==r) return;
  51. int mid=(l+r)>>1;
  52. if(p<=mid) insert(lc(x), l, mid, p);
  53. else insert(rc(x), mid+1, r, p);
  54. }
  55. inline bool exist(int x, int y, int l, int r, int ql, int qr) {
  56. if(ql<=l && r<=qr) return t[y].size - t[x].size > 0;
  57. if(t[y].size - t[x].size == 0) return 0;
  58. int mid = (l+r)>>1, ans=0;
  59. if(ql<=mid) ans |= exist(lc(x), lc(y), l, mid, ql, qr);
  60. if(mid<qr ) ans |= exist(rc(x), rc(y), mid+1, r, ql, qr);
  61. return ans;
  62. }
  63. int fa[N][20], deep[N], L[N], R[N], dfc, ver[N]; // dfc = sam::sz
  64. int tot;
  65. void dfs(int u) {
  66. L[u]=++dfc; ver[dfc]=u;
  67. for(int i=1; (1<<i)<=deep[u]; i++)
  68. fa[u][i] = fa[ fa[u][i-1] ][i-1];
  69. for(int i=h[u];i;i=e[i].ne)
  70. if(e[i].v != fa[u][0]) {
  71. fa[e[i].v][0] = u;
  72. deep[e[i].v] = deep[u]+1;
  73. dfs(e[i].v);
  74. }
  75. R[u]=dfc;
  76. }
  77. void build() {
  78. //for(int i=1; i<=dfc; i++) printf("%d ",ver[i]); puts("");
  79. //for(int i=1; i<=dfc; i++) printf("lr %d %d %d \n",i,L[i],R[i]);
  80. //for(int i=2; i<=dfc; i++) printf("val %d %d\n",i,sam::t[i].val);
  81. for(int i=1; i<=dfc; i++) {
  82. root[i] = root[i-1];
  83. int u=ver[i];
  84. if(suf[u]) insert(root[i], 1, n, suf[u]);// printf("suf %d --> %d --> %d\n", i, u, suf[u]);
  85. }
  86. }
  87. bool check(int len, int l, int r, int x) { //printf("check %d [%d, %d] %d\n", len, l, r, x);
  88. using sam::t;
  89. for(int i=17; i>=0; i--)
  90. if(fa[x][i] && t[ fa[x][i] ].val >= len) x=fa[x][i];
  91. //printf("x %d %d L %d R %d\n",x, t[x].val, L[x], R[x]);
  92. return exist(root[L[x]], root[R[x]], 1, n, l, r);
  93. }
  94. void solve(int a, int b, int c, int d) { //printf("\nsolve %d %d %d %d\n",a,b,c,d);
  95. int l=0, r=min(d-c+1, b-a+1), ans=-1;
  96. while(l<=r) {
  97. int mid=(l+r)>>1;
  98. if(check(mid, a+mid-1, b, pos[d])) ans=mid, l=mid+1;
  99. else r=mid-1;
  100. }
  101. printf("%d\n", ans);
  102. }
  103. }
  104. int main() {
  105. //freopen("in","r",stdin);
  106. //freopen("heoi2016_str.in","r",stdin);
  107. //freopen("heoi2016_str.out","w",stdout);
  108. n=read(); Q=read(); scanf("%s", s+1);
  109. reverse(s+1, s+1+n);
  110. for(int i=1; i<=n; i++) extend(s[i]-'a', i);
  111. for(int i=2; i<= sam::sz; i++) tr::ins(sam::t[i].par, i);
  112. tr::dfs(1); tr::build();
  113. for(int i=1; i<=Q; i++) {
  114. a=n-read()+1; b=n-read()+1; c=n-read()+1; d=n-read()+1;
  115. swap(a, b); swap(c, d);
  116. tr::solve(a, b, c, d);
  117. }
  118. }

sort的80分暴力

  1. #include <iostream>
  2. #include <cstdio>
  3. #include <cstring>
  4. #include <algorithm>
  5. #include <cmath>
  6. using namespace std;
  7. typedef long long ll;
  8. const int N=2e5+5;
  9. inline int read(){
  10. char c=getchar();int x=0,f=1;
  11. while(c<'0'||c>'9'){if(c=='-')f=-1;c=getchar();}
  12. while(c>='0'&&c<='9'){x=x*10+c-'0';c=getchar();}
  13. return x*f;
  14. }
  15. int n, Q, x, op, l, r, t1[N], t2[N];
  16. int *a = t1, *b = t2;
  17. int c[N], SIZE;
  18. void csort(int *a, int m) {
  19. memset(c, 0, SIZE);
  20. for(register int i=1; i<=m; i++) c[a[i]]++;// printf("%d ",a[i]); puts(" sort");
  21. for(register int i=1; i<=n; i++) c[i] += c[i-1];
  22. for(register int i=1; i<=m; i++) b[c[a[i]]--] = a[i];
  23. for(register int i=1; i<=m; i++) a[i]=b[i];// printf("b %d %d\n",i,b[i]);
  24. }
  25. void rever(int l, int r) { while(l<r) swap(a[l++], a[r--]); }
  26. int main() {
  27. //freopen("in","r",stdin);
  28. freopen("heoi2016_sort.in","r",stdin);
  29. freopen("heoi2016_sort.out","w",stdout);
  30. n=read(); Q=read(); SIZE = sizeof(int)*(n+1);
  31. for(int i=1; i<=n; i++) a[i]=read();
  32. for(int i=1; i<=Q; i++) {
  33. op=read(); l=read(); r=read();
  34. csort(a+l-1, r-l+1);
  35. if(op==1) rever(l, r);
  36. //for(int i=1; i<=n; i++) printf("%d ", a[i]); puts("");
  37. }
  38. printf("%d", a[read()]);
  39. }

HEOI2016 题解的更多相关文章

  1. 【题解】P4091 [HEOI2016/TJOI2016]求和

    [题解]P4091 [HEOI2016/TJOI2016]求和 [P4091 HEOI2016/TJOI2016]求和 可以知道\(i,j\)从\(0\)开始是可以的,因为这个时候等于\(0\).这种 ...

  2. 题解 P4093 【[HEOI2016/TJOI2016]序列】

    这道题原来很水的? noteskey 一开始以为是顺序的 m 个修改,然后选出一段最长子序列使得每次修改后都满足不降 这 TM 根本不可做啊! 于是就去看题解了,然后看到转移要满足的条件的我发出了黑人 ...

  3. 【题解】 [HEOI2016]排序题解 (二分答案,线段树)

    题目描述 在2016年,佳媛姐姐喜欢上了数字序列.因而他经常研究关于序列的一些奇奇怪怪的问题,现在他在研究一个难题,需要你来帮助他.这个难题是这样子的:给出一个1到n的全排列,现在对这个全排列序列进行 ...

  4. 题解 P4092 【[HEOI2016/TJOI2016]树】

    参考了皎月半洒花的博客 看到树想到树剖,由于要取距自己到根离自己最近的标记点,刚开始想到线段树里存节点深度,查询时返回最大值.但是这样的话只能得到节点深度,无法得知节点编号,就想倍增乱搞一下,求出标记 ...

  5. BZOJ4553:[HEOI2016/TJOI2016]序列——题解

    https://www.lydsy.com/JudgeOnline/problem.php?id=4553 佳媛姐姐过生日的时候,她的小伙伴从某宝上买了一个有趣的玩具送给他.玩具上有一个数列,数列中某 ...

  6. BZOJ4552:[HEOI2016/TJOI2016]排序——题解

    https://www.lydsy.com/JudgeOnline/problem.php?id=4552 https://www.luogu.org/problemnew/show/P2824 在2 ...

  7. [题解向] Luogu4092 [HEOI2016/TJOI2016]树

    #\(\mathcal{\color{red}{Description}}\) \(Link\) 给定一棵以\(1\)为根的树,有两种操作: \(C: \ \ x\)给点\(x\)打上花标记. \(Q ...

  8. 【题解】Luogu P4091 [HEOI2016/TJOI2016]求和

    原题传送门 \[\begin{aligned} a n s &=\sum_{i=0}^{n} \sum_{j=0}^{i}\left\{\begin{array}{c}{i} \\ {j}\e ...

  9. luogu题解 P4092 【[HEOI2016/TJOI2016]树】树链剖分

    题目链接: https://www.luogu.org/problemnew/show/P4092 瞎扯--\(O(Q \log^3 N)\)解法 这道先yy出了一个\(O(Q \log^3 N)\) ...

随机推荐

  1. Shortest path of the king

    必须要抄袭一下这个代码 The king is left alone on the chessboard. In spite of this loneliness, he doesn't lose h ...

  2. HDU--1213并查集

    题目传送门:HDU--1213 //题意:ignatius过生日,客人来到,他想知道他需要准备多少张桌子.然而一张桌子上面只能坐上相互熟悉的人, //其中熟悉可定义成为A与B认识,B与C认识,我们就说 ...

  3. [国嵌笔记][004][Linux快速体验]

    Linux文件系统 bin目录:可执行的程序 boot目录:与Linux启动相关的文件 dev目录:设备以文件的方式存放 etc目录:配置文件 home目录:用户文件 lib目录:与库相关的文件 ro ...

  4. UE4 分层材质 Layerd Materials

    在UE4中最正规的材质制作流程就像: 建立新材质,并将其调整为达至完美. 在内容浏览器中,建立新材质函数,并将所有材质函数节点复制/粘贴到其中. 将网络连接到新的 Make Material Attr ...

  5. 阿里云邮件服务器怎么设置才能在QQ邮箱访问,互发邮件?

    必须要在阿里云邮上打开IMAP和SMTP IMAP能够代发代收.在线更改.垃圾拦截,比POP3好: 记住打开的协议号: IMAP:143 带SSL:993 SMTP: 25 带SSL:465 前提是能 ...

  6. sublime如何实现函数折叠

    Ctrl+Shift+[ 折叠代码 Ctrl+Shift+] 展开代码 Ctrl+KT 折叠属性 Ctrl+K0 展开所有

  7. Varnish的vcl子程序

    Varnish的VCL子程序 以下内容参考: http://book.varnish-software.com/4.0/ VCL子进程,在其中定制Varnish的行为.VCL子例程可用于:添加自定义标 ...

  8. 【绘图技巧】ps快捷键的用法

    Ctrl+N:新建画布          Ctrl+O:打开对话框 F: 在三种画布中切换 Z:缩放工具(临时)      Ctrl+0:满画面显示 空格:切换到手(临时) Ctrl+":网 ...

  9. Android调试错误-No resource identifier found for attribute 'showAsAction'

    转载自:http://www.bubuko.com/infodetail-498830.html 1.问题描述: 24\YoumiAndroidSdk\demo\offers\res\menu\mai ...

  10. js的join和split

    1,split   将字符串转为数组 第一个参数必须,可以是字符串和正则表达式,表示从该地方开始分割字符串:第二个参数可选,表示返回数组的最大长度. split():把整个字符串塞入一个数组中 spl ...