F-Pairwise Modulo

d数组是来算（x整除y）*y中y比x小的数

s数组是算（x整除y）*y中y比x大的数

（x整除y）*y 看x对于前面大于他的数是枚举，对于前面小于他的数是d树状数组储存。

d中 x整除y表示x中有多少个y 所以每个y对于他的倍数都加上y 比如d【y】+y ，d【y*2】+y...d【y*n】+y；当一个比y大的x要求（x整除y）*y的时候就可以根据d【x】来求

{要还是不明白的话就设y=3，x=10 （10/3)*3=9=3+3+3 =d【3】+d【6】+d【9】这真的很巧妙。这样每个数字对于他后面的数字的贡献就可以确定了。就是这个意思。}

那样就看代码吧。

#include<iostream>
#include<cstdio>
using namespace std;typedef long long ll;
const int N=1e6+10;
ll d[N],s[N],ans,sum,num;     //d  s 是上述的意思；  void add(ll *x,int k,int p)
{while(k<N)x[k]+=p,k+=k&-k;
}
ll ask(ll *x,int k)
{ll ans=0;while(k)ans+=x[k],k-=k&-k;return ans;
}int main()
{int n;scanf("%d",&n);for(int i=0;i<n;i++){scanf("%lld",&num);ans+=1ll*num*i+sum;     sum+=num;ans-=ask(d,num);  //这里是求（x整除y）*y中 小与num的数； for(int j=num;j+num<N;j+=num) //枚举num的倍数 {ans-=1ll*(ask(s,j+num-1)-ask(s,j-1))*j; //枚举num的倍数  s存的是数量，求j-num-1到j-1的数量，这段中那个式子等于j； add(d,j,num);//到时候不懂再说； }add(s,num,1); //num的数量+1； printf("%lld ",ans);}return 0;}

F-Pairwise Modulo相关推荐

codeforces 7.22 F Pairwise Modulo
codeforces 7.22 F Pairwise Modulo 给出n个数的数列a,每个数互不相同且都小于3e5,求出qk=∑1<=i,j<=kaimodajq_k=\sum_{1&l ...
Harbour.Space Scholarship Contest 2021-2022 F. Pairwise Modulo 逆向思维 + 树状数组
传送门文章目录题意思路: 题意给定一个不同数组成的序列aaa,定义pkp_kpk为pk=∑i=1k∑j=1kaimodajp_k=\sum_{i=1}^k\sum_{j=1}^ka_i\bm ...
codeforces1553 F. Pairwise Modulo（数学）
F. Pairwise Modulo 想到了,但又没完全想到..wtcl 首先 pk=pk−1+∑1≤i<kakmodai+∑1≤i<kaimodakp_k=p_{k-1}+\sum_{1 ...
CodeForces - 1553F Pairwise Modulo(数论+树状数组)
题目链接:点击查看题目大意:给出一个长度为 nnn 的序列,求 pk=∑1≤i,j≤kaimodajp_k = \sum_{1 \le i, j \le k} a_i \bmod a_jpk=∑1 ...
Harbour.Space Scholarship Contest 2021-2022 (open for everyone, rated, Div. 1 + Div. 2)
Harbour.Space Scholarship Contest 2021-2022 (open for everyone, rated, Div. 1 + Div. 2) 题号题目知识点 A ...
2017 ACM-ICPC 亚洲区（西安赛区）网络赛 F
f(cos(x))=cos(n∗x) holds for all xxx. Given two integers nnn and mmm, you need to calculate the coef ...
245热图展示微生物组的物种和功能丰度或有无、距离矩阵
245热图展示微生物组的物种和功能丰度或有无本节作者:吴一磊中科院微生物所版本1.0.7,更新日期:2020年8月13日本项目永久地址:https://github.com/YongxinLi ...
Tutorial of Codeforces Round 729 (Div.2)
文章目录 A 题目 C++代码 B 题目分析 C++代码 C 题目 C++代码 A 题目 Example input 5 2 2 3 4 5 3 2 3 4 5 5 5 1 2 4 1 2 3 4 ...
2018 11.1 PION 模拟赛
期望:250 100+100+50 实际:210 80+100+30 期望:100 实际:80 最后:两个点T了.可能是求逆元的方法太慢了,也可能是闲的又加了一个快速乘的原因. #inclu ...
矩阵快速幂2Jzzhu and Sequences Recurrences Contemplation! Algebra Reading comprehension
来!接着写题解,立个flag:矩阵快速幂题解不写完今晚不睡觉! 此题解是接着矩阵快速幂1写的,此篇有哪里写不清楚的地方导致看不懂的可以先看1:https://blog.csdn.net/qq_4536 ...

F-Pairwise Modulo

F-Pairwise Modulo相关推荐

最新文章

热门文章