2020.7.15 --差分

最新推荐文章于 2023-10-13 16:31:22 发布

naygnay

最新推荐文章于 2023-10-13 16:31:22 发布

阅读量131

点赞数

分类专栏：差分

本文链接：https://blog.csdn.net/lingdie_yy/article/details/107351041

版权

差分专栏收录该内容

1 篇文章 0 订阅

订阅专栏

本文探讨了差分操作，它是前缀和的逆运算，允许在进行m次区间加减修改后，以O(n)的时间复杂度获取所有数值。通过分析经典题目如P4552 IncDec Sequence和Namomo contest 2 C.序列，阐述如何将区间加减转化为选择正数减1和负数加1的操作，并讨论如何找到使序列相同数的最小操作数以及在最小操作数下可能序列的数量。

摘要由CSDN通过智能技术生成

两道经典题目

差分,前缀和的逆操作
可以在m次区间+ - 修改后,O(n)的得到所有的值

P4552 IncDec Sequence

排列不等式:
在这里插入图片描述
区间加减操作转换成选一个正数 -1 选一个负数 +1

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;

const int maxl = 5e5 + 10;
const int mod = 1e9 + 7;

int n, m, cas, k, cnt, tot, ans, top;
int a[maxl], b[maxl], c[maxl], s[maxl];
bool in[maxl];
ll ansmi, ansmx;

inline void prework()
{
    scanf("%d", &n);
    for (int i = 1; i <= n; i++)
    {
        scanf("%d", &a[i]);
        c[i] = b[i] = a[i] - a[i - 1];
    }
    c[n + 1] = b[n + 1] = -a[n];
}

inline ll sqr(ll x)
{
    return x * x % mod;
}

inline void mainwork()
{
    ansmi = 0;
    for (int i = n + 1; i >= 1; i--)
        if (b[i] < 0)
            s[++top] = i;
    for (int i = 1; i <= n; i++)
        if (b[i] > 0)
        {
            while (b[i] > 0)
            {
                if (-b[s[top]] > b[i])
                {
                    ansmi = (ansmi + sqr(s[top] - i) * b[i] % mod) % mod;
                    b[s[top]] += b[i];
                    b[i] = 0;
                    break;
                }
                else
                {
                    ansmi = (ansmi + sqr(s[top] - i) * (-b[s[top]]) % mod) % mod;
                    b[i] += b[s[top]];
                    b[s[top]] = 0;
                    top--;
                }
            }
        }
    top = 0;
    ansmx = 0;
    for (int i = 1; i <= n + 1; i++)
        if (c[i] > 0)
            s[++top] = i;
        else
        {
            while (c[i] < 0)
            {
                if (c[s[top]] > -c[i])
                {
                    ansmx = (ansmx + sqr(i - s[top]) * (-c[i])) % mod;
                    c[s[top]] += c[i], c[i] = 0;
                    break;
                }
                else
                {
                    ansmx = (ansmx + sqr(i - s[top]) * c[s[top]]) % mod;
                    c[i] += c[s[top]];
                    c[s[top]] = 0;
                    top--;
                }
            }
        }
}

inline void print()
{
    printf("%lld %lld", ansmi, ansmx);
}

int main()
{
    int t = 1;
    //scanf("%d",&t);
    for (cas = 1; cas <= t; cas++)
    {
        prework();
        mainwork();
        print();
    }
    return 0;
}