#include <bits/stdc++.h>
using namespace std;
constexpr int BUF_SZ = 1 << 15;
inline namespace Input {
    char buf[BUF_SZ];
    int pos;
    int len;
    char next_char() {
        if (pos == len) {
            pos = 0;
            len = static_cast<int>(fread(buf, 1, BUF_SZ, stdin));
            if (!len) { return EOF; }
        }
        return buf[pos++];
    }
    int read_int() {
        char ch;
        int sgn = 1;
        while (!isdigit(ch = next_char())) {
            if (ch == '-') { sgn *= -1; }
        }
        int x = ch - '0';
        while (isdigit(ch = next_char())) { x = x * 10 + (ch - '0'); }
        return x * sgn;
    }
}
inline namespace Output {
    char buf[BUF_SZ];
    int pos;
    void flush_out() {
        fwrite(buf, 1, pos, stdout);
        pos = 0;
    }
    void write_char(const char c) {
        if (pos == BUF_SZ) { flush_out(); }
        buf[pos++] = c;
    }
    void write_int(int x) {
        static char num_buf[100];
        if (x < 0) {
            write_char('-');
            x *= -1;
        }
        int len = 0;
        for (; x >= 10; x /= 10) { num_buf[len++] = static_cast<char>('0' + (x % 10)); }
        write_char(static_cast<char>('0' + x));
        while (len) { write_char(num_buf[--len]); }
        write_char('\n');
    }
    void write_long_long(long long x) {
        static char num_buf[100];
        if (x < 0) {
            write_char('-');
            x *= -1;
        }
        int len = 0;
        for (; x >= 10; x /= 10) { num_buf[len++] = (char) ('0' + (x % 10)); }
        write_char((char) ('0' + x));
        while (len) { write_char(num_buf[--len]); }
        write_char('\n');
    }
    // auto-flush output when the program exits
    void init_output() { assert(atexit(flush_out) == 0); }
}
constexpr int MOD = 1e9 + 7;
long long mul(const long long a, const long long b) {
    return (a % MOD)*(b % MOD) % MOD;
}
int main() {//TODO verstehen
    init_output();
    int n = read_int();
    vector<long long> height(n + 1, 0), width(n + 1, 0);
    vector<long long> prefixSum(n + 1, 0); //Summe von der Breite
    vector<long long> dp(n + 1, 0);
    for (int i = 1; i <= n; ++i) {
        height[i] = read_int();
    }
    for (int i = 1; i <= n; ++i) {
        width[i] = read_int();
        prefixSum[i] = (prefixSum[i - 1] + width[i]) % MOD;
    }
    // absteigende Höhe
    stack<int> monoStack;
    monoStack.push(0);
    long long answer = 0;
    for (int i = 1; i <= n; ++i) {
        while (monoStack.size() > 1 && height[monoStack.top()] >= height[i]) {
            monoStack.pop();
        }
        const int j = monoStack.top();
        // Summe der Breiten von (j+1) .. i
        const long long len = (prefixSum[i] - prefixSum[j] + MOD) % MOD;
        //Summe der Breiten außer der aktuellen
        const long long L1 = (len - width[i] % MOD + MOD) % MOD;
        const long long triHeight = mul( mul(height[i], height[i] + 1), 500000004);
        const long long triWidth = mul( mul(width[i], width[i] + 1), 500000004);
        //T_A,B
        const long long term1 = mul( mul(L1, triHeight), width[i] );
        const long long term2 = mul(triHeight, triWidth);
        const long long A = (term1 + term2) % MOD;
        const long long B = mul(dp[j], width[i]);
        dp[i] = ( dp[j] + mul(len, triHeight) ) % MOD;
        monoStack.push(i);
        answer = (answer + A + B) % MOD;
    }
    write_long_long(answer);
    return 0;
}
| # | Verdict | Execution time | Memory | Grader output | 
|---|
| Fetching results... | 
| # | Verdict | Execution time | Memory | Grader output | 
|---|
| Fetching results... | 
| # | Verdict | Execution time | Memory | Grader output | 
|---|
| Fetching results... | 
| # | Verdict | Execution time | Memory | Grader output | 
|---|
| Fetching results... | 
| # | Verdict | Execution time | Memory | Grader output | 
|---|
| Fetching results... | 
| # | Verdict | Execution time | Memory | Grader output | 
|---|
| Fetching results... | 
| # | Verdict | Execution time | Memory | Grader output | 
|---|
| Fetching results... |