A dd爱科学1.0

主要思路

本质上求最长不下降子序列,dp解法,详细版本可以参考y总算法基础课

f[i]表示以a[i]结尾的最长不下降子序列长度,考虑如何更新f[i]

用j遍历1到i - 1,当w[i] >= w[j]时,更新f[i] = max(f[i], f[j] + 1)

  1. #include <iostream>
  2. #include <algorithm>
  3. using namespace std;
  4. const int N = 1010;
  5. int f[N], n;
  6. char a[N];
  7. int main(void)
  8. {
  9. cin >> n;
  10. for(int i = 1; i <= n; i++) cin >> a[i];
  11. for(int i = 1; i <= n; i++)
  12. {
  13. f[i] = 1;//f[i]最小值为1
  14. for(int j = 1; j < i; j++)
  15. {
  16. if(a[i] >= a[j]) f[i] = max(f[i], f[j] + 1);//更新f[i]
  17. }
  18. }
  19. int res = 0;
  20. for(int i = 1; i <= n; i++) res = max(res, f[i]);
  21. cout << n - res << endl;
  22. return 0;
  23. }

复杂度为O(n6会超时,考虑如何优化

首先考虑每一种长度的不下降子序列我们只需要存储该子序列末尾的最小值即可。

假设以两个元素a[i]和a[j]结尾的最长不下降子序列长度都为len,且a[i] < a[j],那么我们后续在更新其他元素时一定会用a[i]更新,不会用a[j]更新,所以每一种长度的不下降子序列我们只需要存储该子序列末尾的最小值即可。

我们用q[i]存储长度为i的最长不下降子序列结尾的最小值。

继续考虑q[i]数组的性质。

用反证法,假设q[i] > q[i + 1],也就是说长度为i的子序列末尾的最小值 > 长度为i + 1子序列末尾的最小值,那么

我们一定可以用q[i + 1]去更新q[i],所以一定有q[i] <= q[i + 1],所以q[i]数组具有单调性。

我们在更新f[i]时,只需要二分出q数组中第一个值小于等于a[i]的即可,也就是说找到值小于等于a[i]的最长长度。

  1. #include <iostream>
  2. #include <cstring>
  3. #include <string>
  4. #include <algorithm>
  5. #include <string>
  6. using namespace std;
  7. #define int long long
  8. #define debug(a) cout << #a << " = " << a << endl;
  9. #define x first
  10. #define y second
  11. typedef pair<int, int> P;
  12. const int N = 1000010;
  13. int n, f[N], q[N];
  14. signed main(void)
  15. {
  16. std::ios::sync_with_stdio(false);
  17. cin >> n;
  18. string t;
  19. cin >> t;
  20. int len = 0;
  21. for(int i = 0; i < n; i++)
  22. {
  23. int l = 0, r = len;//二分
  24. while(l < r)
  25. {
  26. int mid = l + r + 1 >> 1;
  27. if(q[mid] <= t[i]) l = mid;
  28. else r = mid - 1;
  29. }
  30. len = max(len, r + 1);//二分的结果r就是值小于等于t[i]的最长长度,更新最长长度
  31. q[r + 1] = t[i];//更新当前r+1长度的值。
  32. //因为如果q[r+1]是小于等于t[i]的那么二分的结果应该为r+1,不是r
  33. //所以q[r+1]一定大于t[i],直接赋值即可
  34. }
  35. cout << n - len << endl;
  36. return 0;
  37. }

B dd爱探险

主要思路

状态压缩dp

首先 n个点 n-1次跳完,明白这个跳跃过程其实是一条链,每个点只有被经过和没被经过两个状态,所以典型状压dp。

对于一个数i,我们把它转化为2进制,2进制一共n位,第i位为0表示第i个点没被经过,为1表示第i个点已经被经过

定义f [ i ] [ j ] [k]为当前经过点的状态为i,当前在j这个点上,k为0表示没经历过加速,k为1表示仅经历过重力加速,k为2表示仅经历过反重力加速,k为3表示经历过两次加速。

枚举最后停留的位置j,j范围为1~n,答案为min(f[(1 << n) - 1] [j] [3])。1 << n - 1表示所有点都被经过,有n个1的状态。

状态转移方程见代码

  1. #include <iostream>
  2. #include <cstring>
  3. #include <string>
  4. #include <algorithm>
  5. #include <stack>
  6. #include <string>
  7. #include <cmath>
  8. #include <unordered_map>
  9. #include <queue>
  10. using namespace std;
  11. #define debug(a) cout << #a << " = " << a << endl;
  12. #define x first
  13. #define y second
  14. typedef pair<int, int> P;
  15. const int N = 17;
  16. int n, p[N][N];
  17. int f[1 << N][N][4], st[N];
  18. int a[N], b[N];
  19. int main(void)
  20. {
  21. cin >> n;
  22. for(int i = 1; i <= n; i++)
  23. {
  24. for(int j = 1; j <= n; j++)
  25. {
  26. cin >> p[i][j];
  27. }
  28. }
  29. memset(f, 0x3f, sizeof(f));
  30. for(int i = 1; i <= n; i++)
  31. {
  32. f[1 << (i - 1)][i][0] = 0;
  33. }
  34. for(int t = 1; t < (1 << n); t++)
  35. {
  36. int x = t;
  37. for(int i = 1; i <= n; i++) st[i] = x & 1, x >>= 1;//更新当前状态上第i个点有没有被访问过
  38. int ta = 0, tb = 0;
  39. for(int i = 1; i <= n; i++)
  40. {
  41. if(st[i]) a[++ta] = i;//被访问过的存到a数组
  42. else b[++tb] = i;//没被访问过的存到b数组
  43. }
  44. for(int i = 1; i <= ta; i++)//枚举由哪个已经被经过的点跳到哪个没被经过的点
  45. {
  46. for(int j = 1; j <= tb; j++)
  47. {
  48. //t + (1 << (b[j] - 1))表示将b[j]位修改为1
  49. f[t + (1 << (b[j] - 1))][b[j]][0] = min(f[t + (1 << (b[j] - 1))][b[j]][0], f[t][a[i]][0] + p[a[i]][b[j]]);
  50. //没经历过加速的点只能由没经历过加速的点转移
  51. f[t + (1 << (b[j] - 1))][b[j]][1] = min(f[t + (1 << (b[j] - 1))][b[j]][1], min(f[t][a[i]][1] + p[a[i]][b[j]], f[t][a[i]][0]));
  52. //经过重力加速的点能由没经历过加速的点和经历过重力加速的点转移
  53. f[t + (1 << (b[j] - 1))][b[j]][2] = min(f[t + (1 << (b[j] - 1))][b[j]][2], min(f[t][a[i]][2] + p[a[i]][b[j]], f[t][a[i]][0] + 2 * p[a[i]][b[j]]));
  54. //经历过反重力加速的点能由没经历过加速的点和经历过反重力加速的点转移
  55. f[t + (1 << (b[j] - 1))][b[j]][3] = min(f[t + (1 << (b[j] - 1))][b[j]][3], min(f[t][a[i]][3] + p[a[i]][b[j]], min(f[t][a[i]][1] + 2 * p[a[i]][b[j]], f[t][a[i]][2])));
  56. //经历过两次加速的点能由经历过一次加速的点和经历过两次加速的点转移
  57. }
  58. }
  59. }
  60. int res = 0x3f3f3f3f;
  61. for(int i = 1; i <= n; i++)
  62. {
  63. res = min(res, f[(1 << n) - 1][i][3]);//更新答案
  64. }
  65. cout << res << endl;
  66. return 0;
  67. }

C dd爱科学2.0

f[i] [j]表示更改到第i个位置 更新成的字母为 j 的总偏移量和

f[i] [j]为min(f[i - 1] [j]), j从1开始,换句话说f[i] [j]由第i - 1位中小于等于j的更新而来

  1. #include <iostream>
  2. #include <cstring>
  3. #include <string>
  4. #include <algorithm>
  5. #include <stack>
  6. #include <string>
  7. #include <cmath>
  8. #include <unordered_map>
  9. #include <queue>
  10. using namespace std;
  11. #define int long long
  12. #define debug(a) cout << #a << " = " << a << endl;
  13. #define x first
  14. #define y second
  15. typedef pair<int, int> P;
  16. const int N = 1000010;
  17. int n, f[N][30];
  18. signed main(void)
  19. {
  20. std::ios::sync_with_stdio(false);
  21. cin >> n;
  22. string t;
  23. cin >> t;
  24. for(int i = 0; i < 26; i++) f[0][i] = abs(t[0] - 'A' - i);
  25. for(int i = 1; i < n; i++)
  26. {
  27. int minf = 0x3f3f3f3f;//维护f[i - 1][j]中更新成小于等于j的字母中的最小值
  28. for(int j = 0; j < 26; j++)
  29. {
  30. minf = min(minf, f[i - 1][j]);
  31. f[i][j] = f[i - 1][j] + abs(t[i] - 'A' - j);
  32. }
  33. }
  34. int res = 0x3f3f3f3f;
  35. for(int i = 0; i < 26; i++)
  36. {
  37. res = min(res, f[n - 1][i]);//遍历结尾更新成哪个字母结果最小
  38. }
  39. cout << res << endl;
  40. return 0;
  41. }

第二种dp方式

f[i] [j]表示改至第i位位置,最后一位小于等于j的最小值

则f[i] [j] = min(f[i] [j - 1], f[i - 1] [j] + abs(t[i] - ‘A’ - j));

状态转移方程的含义是f[i] [j]可以由改至第i位位置,最后一位小于等于j - 1的状态改至第i - 1位位置,最后一位小于等于j的状态更新

最后答案位f[n] [26]

  1. #include<bits/stdc++.h>
  2. using namespace std;
  3. int n,ans;
  4. char a[1000005];
  5. int f[1000005][30];
  6. int main(){
  7. scanf("%d",&n);
  8. scanf("%s",a+1);
  9. memset(f,0x3f,sizeof(f));
  10. for (int i=1;i<=26;++i)f[0][i]=0;
  11. for (int i=1;i<=n;++i)
  12. for (int j=1;j<=26;++j)f[i][j]=min(f[i][j-1],f[i-1][j]+abs(a[i]-(j-1+'A')));
  13. printf("%d",f[n][26]);
  14. }

D dd爱矩阵

贪心+优先队列

我们先把问题简化,如果两行,每行 n个数,怎么选和前n大的所有数

把两行分别降序 sort,如果令第一行为数组 a,第二行为数组 b

则可得到最大值为 a[1]+b[1]

那么如何快速求出最大的n个数呢

我们把a[1] + b[j]的和, 1 <= j <= n,和当前a的位置,a[1]记为1,放入一个结构体中存入优先队列中,每次取出和最大的一个数,然后找到a的位置,让a的位置+1的数来替代当前a,同时把更新过的和放入优先队列(因为更新过后下一个前n大的数一定在优先队列中),取n次最大值就能找到前n大的所有数

这是两行的操作,我们进行n-1次就能求出所有前n大的数

复杂度n^2logn

  1. #include <iostream>
  2. #include <cstring>
  3. #include <string>
  4. #include <algorithm>
  5. #include <stack>
  6. #include <string>
  7. #include <cmath>
  8. #include <unordered_map>
  9. #include <queue>
  10. using namespace std;
  11. #define int long long
  12. #define debug(a) cout << #a << " = " << a << endl;
  13. #define x first
  14. #define y second
  15. typedef pair<int, int> P;
  16. const int N = 1000010;
  17. int n;
  18. int a[N], b[N], c[N];
  19. struct node
  20. {
  21. int pos, w;
  22. friend bool operator < (node x,node y) {
  23. return x.w < y.w;
  24. }
  25. };
  26. bool cmp(int a, int b)
  27. {
  28. return a > b;
  29. }
  30. signed main(void)
  31. {
  32. priority_queue<node> q;
  33. cin >> n;
  34. for(int i = 0; i < n; i++) cin >> a[i];
  35. sort(a, a + n, cmp);
  36. for(int i = 2; i <= n; i++)
  37. {
  38. for(int j = 0; j < n; j++)
  39. {
  40. cin >> b[j];
  41. q.push({0, a[0] + b[j]});//把a数组的位置和a[0],b数组的和放入优先队列
  42. }
  43. for(int j = 0; j < n; j++)
  44. {
  45. node t = q.top();//取出队头和最大的元素
  46. q.pop();
  47. c[j] = t.w;//记录一下前n大的元素
  48. int pos = t.pos + 1;//当前位置以及a数组的下一个位置
  49. q.push({pos, t.w - a[t.pos] + a[pos]});
  50. }
  51. for(int j = 0; j < n; j++) a[j] = c[j];
  52. }
  53. for(int i = 0; i < n; i++) cout << c[i] << ' ';
  54. return 0;
  55. }

E dd爱旋转

操作1, 2进行两次都等同于没进行操作,模拟即可

  1. #include <iostream>
  2. #include <cstring>
  3. #include <string>
  4. #include <algorithm>
  5. #include <stack>
  6. #include <string>
  7. #include <cmath>
  8. #include <unordered_map>
  9. #include <queue>
  10. using namespace std;
  11. #define int long long
  12. #define debug(a) cout << #a << " = " << a << endl;
  13. #define x first
  14. #define y second
  15. typedef pair<int, int> P;
  16. const int N = 1010;
  17. int n, q;
  18. int a[N][N];
  19. int ans[N][N];
  20. void print()
  21. {
  22. for(int i = 0; i < n; i++)
  23. {
  24. for(int j = 0; j < n; j++)
  25. {
  26. cout << ans[i][j] << ' ';
  27. }
  28. cout << endl;
  29. }
  30. }
  31. signed main(void)
  32. {
  33. std::ios::sync_with_stdio(false);
  34. cin >> n;
  35. int t1 = 0, t2 = 0;
  36. for(int i = 0; i < n; i++)
  37. {
  38. for(int j = 0; j < n; j++)
  39. {
  40. cin >> a[i][j];
  41. ans[i][j] = a[i][j];
  42. }
  43. }
  44. cin >> q;
  45. for(int i = 0; i < q; i++)
  46. {
  47. int x;
  48. cin >> x;
  49. if(x == 1) t1++;
  50. else t2++;
  51. }
  52. t1 %= 2;
  53. t2 %= 2;
  54. if(t1)
  55. {
  56. int cnt1 = 0, cnt2 = 0;
  57. for(int i = n - 1; i >= 0; i--)
  58. {
  59. for(int j = n - 1; j >= 0; j--)
  60. {
  61. ans[cnt1][cnt2++] = a[i][j];
  62. }
  63. cnt1++;
  64. cnt2 = 0;
  65. }
  66. memcpy(a, ans, sizeof(ans));
  67. }
  68. if(t2)
  69. {
  70. int cnt1 = 0, cnt2 = 0;
  71. for(int i = n - 1; i >= 0; i--)
  72. {
  73. for(int j = 0; j < n; j++)
  74. {
  75. ans[cnt1][cnt2++] = a[i][j];
  76. }
  77. cnt1++;
  78. cnt2 = 0;
  79. }
  80. }
  81. print();
  82. return 0;
  83. }

F dd爱框框

思路1

前缀和+二分

首先二分出最小的区间长度,然后用前缀和遍历即可

  1. #include <iostream>
  2. #include <cstring>
  3. #include <string>
  4. #include <algorithm>
  5. #include <stack>
  6. #include <string>
  7. #include <cmath>
  8. #include <unordered_map>
  9. #include <queue>
  10. using namespace std;
  11. #define int long long
  12. #define debug(a) cout << #a << " = " << a << endl;
  13. #define x first
  14. #define y second
  15. typedef pair<int, int> P;
  16. const int N = 10000010;
  17. int a[N], s[N], n, m;
  18. bool check(int mid)
  19. {
  20. for(int i = 1; i + mid <= n; i++)
  21. {
  22. if(s[i + mid] - s[i - 1] >= m) return true;
  23. }
  24. return false;
  25. }
  26. signed main(void)
  27. {
  28. scanf("%lld%lld", &n, &m);
  29. for(int i = 1; i <= n; i++)
  30. {
  31. scanf("%lld", &a[i]);
  32. s[i] = s[i - 1] + a[i];
  33. }
  34. int l = 1, r = n;
  35. while(l < r)
  36. {
  37. int mid = l + r >> 1;
  38. if(check(mid)) r = mid;
  39. else l = mid + 1;
  40. }
  41. for(int i = 1; i + l <= n; i++)
  42. {
  43. if(s[i + l] - s[i - 1] >= m)
  44. {
  45. cout << i << ' ' << i + l << endl;
  46. break;
  47. }
  48. }
  49. return 0;
  50. }

思路二

前缀和+双指针

  1. #include <bits/stdc++.h>
  2. using namespace std;
  3. const int N = 10000009;
  4. const int mod = 1e9 + 7;
  5. int n, x, a[N];
  6. int main() {
  7. cin >> n >> x;
  8. for(int i = 1; i <= n; i ++) scanf("%d", a + i);
  9. int sum = 0;
  10. int l = 1, r = 1;
  11. int ans = N, ansl = N;
  12. while(r <= n) {
  13. if(sum >= x) {
  14. if(r - l < ans)
  15. ans = r - l, ansl = l;
  16. sum -= a[l ++];
  17. }
  18. else sum += a[r ++];
  19. }
  20. cout << ansl << ' ' << ansl + ans - 1;
  21. return 0;
  22. }

思路三

单调队列维护窗口移动

  1. #include<bits/stdc++.h>
  2. using namespace std;
  3. typedef long long ll;
  4. ll x, a[10000005], s;
  5. int n;
  6. int main(){
  7. scanf("%d%lld",&n, &x);
  8. for (int i = 1; i <= n; i++) scanf("%lld",&a[i]);
  9. int l = 0, r = n + 1;//表示答案
  10. int start = 1, end = 0;
  11. queue<int> q;//队列
  12. for(int i = 1; i <= n; i++)
  13. {
  14. s += a[i]; end = i;
  15. q.push(a[i]);
  16. while(q.size() && s - q.front() >= x) s -= q.front(), q.pop(), start++;
  17. //如果当前s - q.front() >= x那么将队头元素弹出队列,更新start,同时更新s
  18. if(s >= x)
  19. {
  20. if(r - l > end - start)//如果距离变小则更新
  21. {
  22. l = start;
  23. r = end;
  24. }
  25. }
  26. }
  27. printf("%d %d\n", l, r);
  28. }