联想杯 - G. Gentle Jena

最新推荐文章于 2024-11-08 14:08:42 发布

李峻枫

最新推荐文章于 2024-11-08 14:08:42 发布

阅读量100

点赞数

分类专栏： OI&ACM # 单调栈文章标签：数据结构单调栈

本文链接：https://blog.csdn.net/lijf2001/article/details/119217920

版权

OI&ACM 同时被 2 个专栏收录

397 篇文章 2 订阅

订阅专栏

单调栈

5 篇文章 0 订阅

订阅专栏

该博客介绍了如何利用单调栈解决一个关于数组前缀和与区间最值的问题。题目中要求计算特定函数的权值之和，并对一定模数取反，博主通过观察发现可以利用单调栈优化计算过程，将时间复杂度降低到线性。代码实现中展示了如何运用单调栈来快速计算每个位置的区间最值之和，从而高效地得出最终答案。

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

题目大意

给出一个长度为 $n$ 的数组 $b_1,b_2,\cdots,b_n$ 。
定义函数 $f(l,r)=\displaystyle\sum_{i=l}^r \min \{b_i\}$
定义权值 $A_i=\displaystyle\sum_{1\le l\le r\le i} f(l,r)$
$A_i$ 对 $998244353$ 取模。
求 $A_1\bigoplus A_2\bigoplus\cdots\bigoplus A_n$

时间限制

数据范围

$n\le 10^7$

题解

考虑如何求 $f$ ，如果弄个线段树或者别的，也十分不方便。
再来考虑如何求 $A_i$ ，不难发现 $A_i$ 只比 $A_{i-1}$ 多了以 $i$ 为右端点的区间的 $f$ 值。
发现了这一点，就非常好办了。

考虑如何快速求出 $f(1,i),f(2,i),\cdots,f(i,i)$
根据性质，就可以发现这些值是递减的，因此可以用一个单调栈维护。
整体时间复杂度就是 $O (n)$

Code

//#pragma GCC optimize (2)
//#pragma G++ optimize (2)
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <iostream>
#include <queue>
#include <map>
#define ll long long
#define G getchar
using namespace std;

ll read()
{
    char ch;
    for(ch = G();ch < '0' || ch > '9';ch = G());
    ll n = 0;
    for(;'0' <= ch && ch <= '9';ch = G())n = (n<<1)+(n<<3)+ch-48;
    return n;
}

void write(ll x)
{
    if (x > 9)
    {
        write(x / 10);
        putchar(x % 10 + 48);
    }
    else putchar(x + 48);
}

const int N = 10000007;
const int mo = 998244353;

ll n , p , x , y , z , b , a;
ll ans , sum , s;
ll zz[N] , id[N] , top;

int main()
{
    //freopen("i.in","r",stdin);
    //freopen("e.out","w",stdout);
    
    n = read();
    p = read();
    x = read();
    y = read();
    z = read();
    b = read();
    a = ans = b;
    zz[1] = b;
    id[1] = 1;
    top = 1;
    s = sum = b;

    for (int i = 2 ; i <= n ; i++)
    {
        b = (a * x + b * y + z) % p;
        for ( ; top && zz[top] > b ; top--)
        {
            sum = sum - zz[top] * (id[top] - id[top - 1]) % mo;
            if (sum < 0) sum = sum + mo;
        }

        top++;
        zz[top] = b;
        id[top] = i;
        sum = sum + b * (id[top] - id[top - 1]) % mo;
        if (sum >= mo) sum = sum - mo;
        s = s + sum;
        if (s >= mo) s = s - mo;
        a = s;

        ans = ans ^ a;
        //printf("%lld %lld\n", b , a);
    }

    printf("%lld\n", ans);

    return 0;
}