abc058d <公式化简>

合集 - ACM-abc(68)
1.abc306e <mutiset的使用>2.abc053d <简单贪心>3.abc042d <组合数>4.abc043d5.abc045d6.abc044d <根号分治>7.abc046d8.abc047d9.abc048d <博弈>10.abc049d <并查集>11.abc050d <???>12.abc051 <多源最短路>13.abc052d14.abc054d <dp, 背包>15.abc055d <枚举>16.abc056d <贪心 / 二分+DP+bitset>17.abc057d <贪心+组合计数>
18.abc058d <公式化简>
19.abc059d <博弈, 打表找规律>20.abc060d <dp, 背包>21.abc061d <单源最短路, spfa, 判断负环>22.abc062d <优先队列>23.abc063d <二分答案>24.abc064d <贪心/前缀和>25.abc309f <线段树 + 离散化 + 双指针>26.abc309e <dfs>27.abc065d <贪心+最小生成树> [lambda表达式]28.abc066d <组合>29.abc067d <博弈 + dfs>30.abc068d <思维 + 构造>31.abc069d <构造>32.abc070d <简单树上dfs>33.abc071d <递推>34.abc072d <贪心>35.abc073d <Floyed + 枚举排列>36.abc074d <Floyed 消除传递边>37.abc075d <暴力枚举 / 枚举+离散化+二维前缀和>38.abc076d <dp / 贪心>39.abc077d <思维 + 最短路 (将构造数字过程视为最短路)>40.abc078d <博弈>41.abc079d <Floyed>42.abc080d <区间重叠>43.abc081d <思维 构造>44.abc082d <bitset 状压dp>45.abc083d <思维 贪心>46.abc084d <素数筛 前缀和>47.abc085d <贪心>48.abc086d <二维前缀和 同余>49.abc087d <并查集 维护距离信息>50.abc088 <bfs 最短路>51.abc089 <前缀和>52.abc310d <dfs暴搜-分组方案数 / bitmask表示集合+dp>53.abc310e <公式递推(dp?)>54.abc310f <dp + bitmask>55.abc090d <枚举计数>56.abc312c <二分答案>57.abc312d <dp, 括号匹配方案数>58.abc312e <暴力>59.abc320f <dp >60.abc092d<构造，思维>61.abc094d<组合数>62.abc095d<思维>63.abc096d<素数筛，整除>64.abc097d<并查集，排列>65.abc098d<双指针，异或>66.abc099d<dfs，枚举排列方案>67.abc100d<枚举>68.abc101d<打表，数学>
收起

D - ###

原计算公式为:

\[\sum\limits_{1\le i<j\le n}\sum\limits_{1\le k<l\le m}(x_j-x_i)(y_l-y_k) \]

可将xy拆分:

\[\left(\sum\limits_{1\leq i<j\leq n}(x_j-x_i)\right)\left(\sum\limits_{1\leq k<l\leq m}(y_l-y_k)\right) \]

仅计算x侧可进一步化简:

\[\sum\limits_{1\le i<j\le n}(x_j-x_i)=\sum\limits_{1\le k\le n}((k-1)x_k-(n-k)x_k) \]

// https://atcoder.jp/contests/abc058/tasks/arc071_b
// 公式化简

#include <iostream>
#include <algorithm>
using namespace std;
typedef long long LL;
const int N = 1e5 + 10, mod = 1e9 + 7;
int n, m;
int x[N], y[N];

LL qpow(LL a, LL k)
{
    LL res = 1;
    a %= mod;
    while (k)
    {
        if (k & 1) res = res * a % mod;
        k >>= 1;
        a = a * a % mod;
    }
    return res;
}

LL C(LL n, LL m)
{
    if (m < 0 || m > n) return 0;
    if (m > n / 2) return C(n, n - m);
    LL res = 1;
    LL t = 1;
    for (int i = n, j = 1; j <= m; j ++, i --)
    {
        res = res * i % mod;
        t = t * j % mod;
    }
    return res * qpow(t, mod-2) % mod;
}

void solv()
{
    cin >> n >> m;
    for (int i = 1; i <= n; i ++) cin >> x[i];
    for (int i = 1; i <= m; i ++) cin >> y[i];
    sort(x+1, x+1+n); sort(y+1, y+1+m);
    LL xx = 0, yy = 0;
    for (LL i = 1; i <= n; i ++) xx = (xx + (i*2 - 1 - n) * x[i] % mod) % mod;
    for (LL i = 1; i <= m; i ++) yy = (yy + (i*2 - 1 - m) * y[i] % mod) % mod;
    LL ans = xx * yy % mod;
    cout << ans << endl;
}

int main()
{
    ios::sync_with_stdio(false); cin.tie(0); cout.tie(0);
    int T = 1;
    // cin >> T;
    while (T --)
    {
        solv();
    }
    return 0;
}
posted @   O2iginal  阅读(16)  评论(0)  编辑  收藏  举报
相关博文：
·  abc080d <区间重叠>
·  abc052d
·  AtCoder Grand Contest 025
·  D - I Wanna Win The Game
·  AtCoder Beginner Contest 281 (A-D)
阅读排行：
· 震惊！C++程序真的从main开始吗？99%的程序员都答错了
· 【硬核科普】Trae如何「偷看」你的代码？零基础破解AI编程运行原理
· 单元测试从入门到精通
· 上周热点回顾（3.3-3.9）
· winform 绘制太阳，地球，月球 运作规律
点击右上角即可分享
微信分享提示
AI IDE Trae