• D. Make Them Equal(dp + 范围优化 )


    Problem - 1633D - Codeforces

    你有一个大小为n的整数数组a,初始时,数组中的所有元素都等于1,你可以进行以下操作:选择两个整数i(1≤i≤n)和x(x>0),然后将ai的值增加⌊aix⌋(即使ai=ai+⌊aix⌋)。

    进行所有操作后,你将为所有ai=bi的这种i收到ci币。

    你的任务是确定执行不超过k个操作所能得到的最大数量的硬币。

    输入
    第一行包含一个整数t(1≤t≤100)--测试案例的数量。

    每个测试案例的第一行包含两个整数n和k(1≤n≤103;0≤k≤106)--分别是数组的大小和最大操作数。

    第二行包含n个整数b1,b2,...,bn(1≤bi≤103)。

    第三行包含n个整数c1,c2,...,cn(1≤ci≤106)。

    所有测试案例的n之和不超过103。

    输出
    对于每个测试案例,打印一个整数--你通过执行不超过k个操作所能得到的最大硬币数。

    例子
    输入复制
    4
    4 4
    1 7 5 2
    2 6 5 2
    3 0
    3 5 2
    5 4 7
    5 9
    5 2 5 6 3
    5 9 1 9 7
    6 14
    11 4 6 2 8 16
    43 45 9 41 15 38
    输出拷贝
    9
    0
    30
    167
    题解:
    首先我们要考虑每个1最少要经过几步才能够到b[i],这其实也是个dp

    由于n只有1000,直接枚举1000内所有i

    转移方程为

    注意初始化

        if(i + i/j <= 1000)
         {
                 v[i + i/j] = min(v[i + i/j] ,v[i] + 1);
         }

     剩下其实就是一个01背包,k内范围如何得到最多代价,但是k有点大有1e6(CF平台有点nb,这1e3*1e6都能跑完)

    但我们还是想想如何优化吧,仔细看看我们发现v[i]最大值不超过13,所以我们直接

    k = min(k,n*13)比n*13大的是不需要的

    这样我们时间复杂度就成了1e3*1e4

    1. #include<iostream>
    2. #include<algorithm>
    3. #include<cstring>
    4. #include<string>
    5. #include<map>
    6. #include<vector>
    7. #include<queue>
    8. using namespace std;
    9. #define int long long
    10. int b[1004];
    11. int c[1004];
    12. int v[1004];
    13. int dp[1000040];
    14. void solve()
    15. {
    16. int n,k;
    17. cin >> n >> k;
    18. for(int i = 1;i <= n;i++)
    19. {
    20. cin >> b[i];
    21. }
    22. for(int i = 1;i <= n;i++)
    23. cin >> c[i];
    24. for(int i = 0;i <= k;i++)
    25. dp[i] = 0;
    26. k = min(k,n*13);
    27. for(int i = 1;i <= n;i++)
    28. {
    29. for(int j = k;j >= v[b[i]];j--)
    30. {
    31. dp[j] = max(dp[j],dp[j - v[b[i]]] + c[i]);
    32. }
    33. }
    34. cout<< dp[k] << "\n";
    35. }
    36. signed main()
    37. {
    38. // ios::sync_with_stdio(false);
    39. // cin.tie(0);
    40. // cout.tie(0);
    41. int t = 1;
    42. cin >> t;
    43. memset(v,0x3f,sizeof v);
    44. v[1] = 0;
    45. for(int i = 1;i <= 1000;i++)
    46. {
    47. for(int j = 1;j <= i;j++)
    48. {
    49. if(i + i/j <= 1000)
    50. {
    51. v[i + i/j] = min(v[i + i/j] ,v[i] + 1);
    52. }
    53. }
    54. }
    55. while(t--)
    56. {
    57. solve();
    58. }
    59. }
    60. //1 10 11
    61. //001
    62. //010
    63. //011
    64. //100

  • 相关阅读:
    iOS——引用计数(一)
    混合网状防火墙的兴起如何彻底改变网络安全
    怎样给蔬菜图片抠图?看看这个马路中央的紫茄子
    带你区分几种并行
    SpringCloud微服务实战——搭建企业级开发框架(三十八):搭建ELK日志采集与分析系统
    nginx反向代理,负载均衡配置
    跨域的五种解决方案
    Python爬虫之xpath语法及案例使用
    Apache Spark 的基本概念
    7-15 求矩阵的局部极大值
  • 原文地址:https://blog.csdn.net/m0_64158084/article/details/128091978