首页 > 其他分享 >[ABC351F] Double Sum

[ABC351F] Double Sum

时间:2024-08-15 20:27:36浏览次数:12  
标签:200005 int Double Sum st next low now ABC351F

原题链接

题解

方法一:
双重循环, \(O(n^2)\)

方法二:

顺序遍历 \(i\),然后查找目前所有比 \(a_i\) 小的数,这是一个比较经典的树状数组的运用

时间复杂度 \(P(n\log A)\)

考虑优化,由于 \(A\) 可以达到 \(1e8\) ,而 \(n\) 只有 \(4e5\) ,所以我们可以对数据做离散化处理

code

#include<bits/stdc++.h>
using namespace std;
/*
mt19937_64 rnd(time(0));
#define double long double
const int inf=1e18;
const int mod=1e9+7;

const int N=4e5;
int qpow(int a,int n)
{
    int res=1;
    while(n)
    {
        if(n&1) res=res*a%mod;
        a=a*a%mod;
        n>>=1;
    }
    return res;
}
int inv(int x)
{
    return qpow(x,mod-2);
}
int fa[2000005];
int finds(int now) { return now == fa[now] ? now :fa[now]=finds(fa[now]); }

vector<int> G[200005];

int dfn[200005],low[200005];
int cnt=0,num=0;
int in_st[200005]={0};
stack<int> st;
int belong[200005]={0};

void scc(int now,int fa)
{
    dfn[now]=++cnt;
    low[now]=dfn[now];
    in_st[now]=1;
    st.push(now);

    for(auto next:G[now])
    {
        if(next==fa) continue;

        if(!dfn[next])
        {
            scc(next,now);
            low[now]=min(low[now],low[next]);
        }
        else if(in_st[next])
        {
            low[now]=min(low[now],dfn[next]);
        }
    }

    if(low[now]==dfn[now])
    {
        int x;
        num++;
        do
        {
            x=st.top();
            st.pop();
            in_st[x]=0;
            belong[x]=num;
        }while(x!=now);
    }
}
vector<int> prime;
bool mark[200005]={0};
void shai()
{
    for(int i=2;i<=200000;i++)
    {
        if(!mark[i]) prime.push_back(i);

        for(auto it:prime)
        {
            if(it*i>200000) break;

            mark[it*i]=1;
            if(it%i==0) break;
        }
    }
}
*/
#define lowbit(x) ((x)&(-x))
#define int long long

int tree1[400005],tree2[400005];
int n;
void update1(int x)
{
    while(x<=n)
    {
        tree1[x]++;;
        x+=lowbit(x);
    }
}

void update2(int x,int v)
{
    while(x<=n)
    {
        tree2[x]+=v;
        x+=lowbit(x);
    }
}

int query1(int x)
{
    int res=0;
    while(x)
    {
        res+=tree1[x];
        x-=lowbit(x);
    }
    return res;
}

int query2(int x)
{
    int res=0;
    while(x)
    {
        res+=tree2[x];
        x-=lowbit(x);
    }
    return res;
}
int a[400005],b[400005];
void solve()
{
    cin>>n;

    for(int i=1;i<=n;i++)
    {
        cin>>a[i];
        b[i]=a[i];
    }
    sort(b+1,b+1+n);
    int len=unique(b+1,b+1+n)-b-1;

    int ans=0;
    for(int i=1;i<=n;i++)
    {
        int pos=lower_bound(b+1,b+1+len,a[i])-b;
        ans+=query1(pos-1)*a[i]-query2(pos-1);
        update1(pos);
        update2(pos,a[i]);
    }

    cout<<ans;
}
signed main()
{
    ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
    int TT=1;
    //cin>>TT;
    while(TT--) solve();
    return 0;
}


标签:200005,int,Double,Sum,st,next,low,now,ABC351F
From: https://www.cnblogs.com/pure4knowledge/p/18361766

相关文章

  • md5sum+可执行文件 怎么理解?
    `md5sum`是一个在Unix、Linux以及其他类Unix系统中广泛使用的命令行工具,用于计算和校验文件的MD5哈希值。MD5哈希是一种广泛使用的加密哈希函数,可以产生一个128位(16字节)的哈希值(通常以32位的十六进制数表示),用于验证文件的完整性和一致性。 当你看到`md5sum`与“可......
  • 摘要生成—通过摘要风格控制摘要的生成/抽取,原文阅读与理解:GEMINI: Controlling The S
    GEMINI:ControllingTheSentence-LevelSummaryStyleinAbstractiveTextSummarizationGEMINI:在抽象文本摘要中控制句子级摘要风格paper:https://arxiv.org/abs/2304.03548github:https://github.com/baoguangsheng/gemini本文介绍了一种自适应摘要抽取/生成方......
  • AGC182 C - Sum of Number of Divisors of Product
    题目大意:求长度为1,2,...N,的好序列因数个数。好序列满足每个元素\(1\leqx\leqM\)\(N\leq1e18,M\leq16\)很容易想到维护所有好序列质因数的指数+1的乘积。\(\prodb_i,1\leqi\leq6\).考虑每个数对这个乘积的影响。假设我们只有2,3,5这三个质数,序列末端添加15,......
  • 题解:AT_abc351_f [ABC351F] Double Sum
    关于某c人士强制偷袭代码导致AT号被封,\({\color{red}\mathrm{警钟敲碎}}\)。题意一个长\(n\)的数组\(a\),求所有顺序对中两元素之差的和。分析听说有大佬2min切掉。很明显,暴力过不去。考虑计算每个元素的贡献。设\(id\)为该元素之前所有比它小的元素个数,\(sum\)表示这些......
  • 爱因斯坦求和约定einsum简单例题解读
    概论在爱因斯坦求和约定或einsum()格式字符串中,所有的索引都可以分为两类:自由索引集和求和索引集。它们的区别很简单:自由索引是用于输出规范中的索引。它们与外层for循环相关联。求和索引是所有其他索引:它们出现在参数规范中,但不出现在输出规范中。之所以称为求和索引,是因......
  • 举例说明二次型和用 einsum 计算
    什么是“二次型”在数学中,特别是线性代数中,二次型(quadraticform)是一个涉及向量和矩阵的表达式,其形式为:[Q(v)=v^TMv]其中:(v)是一个向量(长度为(n))。(M)是一个(n\timesn)的方阵(矩阵)。(v^T)表示向量(v)的转置(即一个列向量变为行向量)。这个表达......
  • [ABC366D] Cuboid Sum Query 题解
    [ABC366D]CuboidSumQuery题解原题传送门AT原题传送门题意翻译:给予一个\(N\timesN\timesN\)的三维矩阵,有\(Q\)次询问,对于每次询问,给与四个数,分别为\(L_1,R_1,L_2,R_2,L_3,R_3\)求在三维矩阵中\(a[L_1][L_2][L_3]\)到\(a[R_1][R_2][R_3]\)的区间和。三维前缀......
  • IT基础书籍汇集_sum
    希望STUDENT过软考,所以有些基础书籍还是需要看看--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------虽然书籍,标注“著”的书籍一般......
  • Python - sum()
     >>>bugs=["bug1","bug2"]>>>sum(bugs,[])Traceback(mostrecentcalllast):File"<stdin>",line1,in<module>TypeError:canonlyconcatenatelist(not"str")tolist>&g......
  • Range Minimum Sum
    非常经典的删数问题,见这篇题解我赛时的时候考虑的时候删除了\(a_i\)后,有哪些区间会被删除,哪些区间会被加入删除的区间:最小值是\(a_i\)的区间(\(O(1)\)计算)、\(a_i\)作为一个端点但是\(a_i\)不是最小值的区间(差分维护)加入的区间:左端点属于\((l_i,i)\)且右端点属于\((i,r_i)\)的区......