• 差分,前缀和,离散化——模板


    1.一维差分 

    我们给定一个数组a[2,4 ,6 ,8 ,10],
    则它的差分数组b[2,2,2,2,2,2]

    假设我们让数组区间[2,4]元素都增加1
    则数组a[2,5,7,9,10]
    差分数组b[2,3,2,2,1]

    我们可以发现规律,对于差分数组b,对于任意区间[l,r],增加x
    它的数组变化为b[l]+=x, b[r+1]-=x
    ;

    1. #include
    2. using namespace std;
    3. //#define get_sum(L,R) (L?sum[R]-sum[L-1]:sum[R])
    4. const int n=6;
    5. int d[6]={0};
    6. void add(int L, int R,int V)
    7. {
    8. d[L]+=V;
    9. d[R+1]-=V;
    10. }
    11. int main()
    12. {
    13. int arr[n]={1,3,7,5,2};
    14. add(2,4,5);
    15. add(1,3,2);
    16. add(0,2,-3);
    17. for(int i=1;i
    18. {
    19. d[i]=d[i-1]+d[i];
    20. }
    21. for(int i=0;i<5;i++)
    22. {
    23. arr[i]=arr[i]+d[i];
    24. cout<" ";
    25. }
    26. return 0;
    27. }

    2.二维差分

    来使被选中的子矩阵中的每个元素的值加上c
    b[x1][y1] + = c;
    b[x1,][y2+1] - = c;
    b[x2+1][y1] - = c;
    b[x2+1][y2+1] + = c;

    每次对b数组执行以上操作,等价于:

    1. for(int i=x1;i<=x2;i++)
    2. for(int j=y1;j<=y2;j++)
    3. a[i][j]+=c;

     

    1. #include
    2. using namespace std;
    3. const int n=4,m=5;
    4. int g[n][m]={{0,0,0,0,0},
    5. {0,1,5,6,8},
    6. {0,9,6,7,3},
    7. {0,5,3,2,4}};
    8. //#define out_sum(x1,y1,x2,y2) (sum[x1,y1]-sum[x2][y1-1]-sum[x1-1][y2]+sum[x1-1][y1-1]);//
    9. void add_d(int x1,int y1,int x2,int y2,int v,int (*d)[5])
    10. {
    11. d[x1][y1]=d[x1][y1]+v;
    12. if(y2+1<5)
    13. d[x1][y2+1]-=v;//注意越界的影响,数组分配的内存是一段连续的空间
    14. d[x2+1][y1]-=v;
    15. d[x2+1][y2+1]+=0;
    16. }
    17. int main()
    18. {
    19. int d[n][m]={0};
    20. int sum[n][m]={0};
    21. add_d(1,1,3,4,2,d);
    22. for(int i=0;i<4;i++)
    23. {
    24. for(int j=0;j<5;j++)
    25. {
    26. cout<" ";
    27. }
    28. cout<
    29. }
    30. cout<
    31. for(int i=1;i<4;i++)
    32. {
    33. sum[i][0]=sum[i-1][0]+d[i][0];
    34. }
    35. for(int i=1;i<5;i++)
    36. {
    37. sum[0][i]=sum[0][i-1]+d[0][i];
    38. }
    39. for(int i=1;i<4;i++)
    40. {
    41. for(int j=1;j<5;j++)
    42. {
    43. sum[i][j]=sum[i][j-1]+sum[i-1][j]-sum[i-1][j-1]+d[i][j];
    44. cout<" ";
    45. }
    46. cout<
    47. }
    48. cout<
    49. for(int i=1;i<4;i++)
    50. {
    51. for(int j=1;j<5;j++)
    52. {
    53. g[i][j]=g[i][j]+sum[i][j];
    54. cout<" ";
    55. }
    56. cout<
    57. }
    58. return 0;
    59. }

    3.一维前缀和

    1. #include
    2. using namespace std;
    3. //#define get_sum(L,R) (L?sum[R]-sum[L-1]:sum[R])
    4. const int n=5;
    5. int sum[n];
    6. int get_sum(int L, int R)
    7. {
    8. if(L>0)
    9. {
    10. return sum[R]-sum[L-1];
    11. }
    12. else
    13. {
    14. return sum[R];
    15. }
    16. }
    17. int main()
    18. {
    19. int arr[n]={1,3,7,5,2};
    20. sum[0]=arr[0];
    21. for(int i=1;i
    22. {
    23. sum[i]=sum[i-1]+arr[i];//公式
    24. }
    25. cout<<get_sum(2,4);
    26. return 0;
    27. }

    4.二维前缀和

    1. #include
    2. using namespace std;
    3. const int n=4,m=5;
    4. int g[n][m]={{0,0,0,0,0},
    5. {0,1,5,6,8},
    6. {0,9,6,7,3},
    7. {0,5,3,2,4}};
    8. //#define out_sum(x1,y1,x2,y2) (sum[x1,y1]-sum[x2][y1-1]-sum[x1-1][y2]+sum[x1-1][y1-1]);//
    9. int out_sum(int x1,int y1,int x2,int y2,int (*sum)[5])
    10. {
    11. return sum[x2][y2]-sum[x2][y1-1]-sum[x1-1][y2]+sum[x1-1][y1-1];
    12. }
    13. int main()
    14. {
    15. int sum[n][m]={0};
    16. for(int i=1;i<4;i++)
    17. {
    18. for(int j=1;j<5;j++)
    19. {
    20. sum[i][j]=sum[i][j-1]+sum[i-1][j]-sum[i-1][j-1]+g[i][j];;
    21. cout<" ";
    22. }
    23. cout<
    24. }
    25. cout<<out_sum(2,2,3,3,sum);
    26. return 0;
    27. }

    5.离散化

    1. #include
    2. #include
    3. using namespace std;
    4. int main()
    5. {
    6. int a[200000+10];
    7. int n;
    8. cin>>n;
    9. for(int i=1;i<=n;i++)
    10. {
    11. cin>>a[i];
    12. }
    13. //排序,左闭右开
    14. sort(a+1,a+n+1);
    15. //unique函数去重,本质是把重复的值放到数组末尾
    16. //函数返回一个迭代器,指向去重后不重复序列的最后那一个元素
    17. int len=unique(a+1,a+n+1)-(a+1);//得到去重后数组长度
    18. /*手动去重
    19. int cnt=0;
    20. for(int i=1;i<=n;i++)
    21. {
    22. if(i==0||a[i]!=a[i-1])
    23. {
    24. a[++cnt]=a[i];
    25. }
    26. }*/
    27. for(int i=1;i<=len;i++)
    28. {
    29. cout<" ";
    30. //在有序数组中查找大于等于啊a[i]的第一个数的下标
    31. a[i]=lower_bound(a+1,a+n+1,a[i])-a;
    32. cout<'\n';
    33. }
    34. }

  • 相关阅读:
    xss的绕过
    SLAM视觉里程计的理解
    ESPHome 通过http获取设备的状态
    平安城市解决方案-最新全套文件
    PyTorch 2.0 重磅发布:一行代码提速 30%
    4 Paimon数据湖之Hive Catalog的使用
    k8s的“前端”——ingress
    SpringMVC(四万五字超详细笔记)
    2.K3S+Rancher之在线安装
    1.77亿美元,安世被迫出售晶圆大厂NWF | 百能云芯
  • 原文地址:https://blog.csdn.net/weixin_62599885/article/details/126166412