51nod1712 区间求和

http://www.51nod.com/Challenge/Problem.html#problemId=1712
先考虑题面中的简化问题。
对于\(i\in [1,n]\)\(a_i\)的贡献为\(a_i*(i-1)-a_i*(n-i)\)
那么对于\(i\in [l,r](a_l=a_r)\),贡献为\(a_i*(i-l)-a_i*(r-i)=2i*a_i-a_i(l+r)\)
这个式子只需要统计4个东西就可以\(O(n)\)计算了。

  • \(i\)左侧\(a_l\)的个数\(A\)
  • \(i\)右侧\(a_r\)的个数\(B\)
  • \(i\)左侧\(a_l\)\(l\)之和\(C\)
  • \(i\)右侧\(a_r\)\(r\)之和\(D\)

那么对答案的贡献就是\(2*i*a_i*A*B-a_i*(BC+AD)\)
统计一下\(A,B,C,D\)即可。

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

inline int read() {
    int x = 0, f = 1; char c = getchar();
    while(c < '0' || c > '9') {if(c == '-') f = -1; c = getchar();}
    while(c >= '0' && c <= '9') {x = x * 10 + c - '0'; c = getchar();}
    return x * f;
}

typedef unsigned int uint;
const int N = 1000010;

int n, lim, a[N];
uint num[N], sum[N], numl[N], numr[N], suml[N], sumr[N];
uint Sum, Num;

int main() {
    n = read();
    for(int i = 1; i <= n; ++i) {
        a[i] = read();
        numr[a[i]]++; sumr[a[i]] += i; 
    }
    uint ans = 0;
    for(int i = 1; i <= n; ++i) {
        numl[a[i]]++; suml[a[i]] += i;
        Num -= num[a[i]];
        num[a[i]] = numl[a[i]] * numr[a[i]];
        Num += num[a[i]];
        Sum -= sum[a[i]];
        sum[a[i]] = suml[a[i]] * numr[a[i]] + sumr[a[i]] * numl[a[i]];
        Sum += sum[a[i]];
        
        ans += 2 * i * a[i] * Num - a[i] * Sum;
        
        numr[a[i]]--; sumr[a[i]] -= i;
        Num -= num[a[i]];
        num[a[i]] = numl[a[i]] * numr[a[i]];
        Num += num[a[i]];
        Sum -= sum[a[i]];
        sum[a[i]] = suml[a[i]] * numr[a[i]] + sumr[a[i]] * numl[a[i]];
        Sum += sum[a[i]];
    }
    printf("%u\n", ans);
}

转载于:https://www.cnblogs.com/henry-1202/p/11494374.html

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值