题目链接

思路分析

看到题目中 \(n,m \leq 1000\) ,故直接考虑 \(O(n^2)\) 级别做法。

我们先把所有的点按照 \(val\) 值从小到大排序,这样的话二维问题变成序列问题。

设 \(f_i\) 表示走到第 \(i\) 个点的价值的期望。

先列出裸的 \(dp\) 方程:(\(Num\) 表示符合条件的点的个数)

\[f_i =\frac{1}{Num} \sum_{a_i > a_j}(x_i-x_j)^2+(y_i-y_j)^2+f_j
\]

但是这个好像是 \(O(n^2m^2)\) 的优秀算法……

不要担心,我们把式子化简一下:

\[f_i=\frac{1}{Num}\sum_{a_i>a_j}x_i^2+2x_ix_j+x_j^2+y_i^2+2y_iy_j+y_j^2+f_j \\
f_i=\frac{1}{Num}\sum_{a_i>a_j}x_i^2+x_j^2+y_i^2+y_j^2+2(x_ix_j+y_iy_j)+f_j
\]

我们惊喜地发现,可以用前缀和优化。

设:

  • \(suma_i=\sum_{j=1}^ia_j\quad sumb_i=\sum_{j=1}^ib_j\)
  • \(sumpa_i=\sum_{j=1}^ia_j^2 \quad sumpb_i=\sum_{j=1}^ib_j^2\)
  • \(sumf_i=\sum_{j=1}^if_j\)

直接按照上述前缀和替换即可,在这里不写了。

Code

#include <bits/stdc++.h>

#define file(a) freopen(a".in", "r", stdin), freopen(a".out", "w", stdout)

#define Enter putchar('\n')
#define quad putchar(' ') #define int long long #define N 1005
#define mod 998244353 int n, m, a[N][N], tot, x, y, f[N * N], sumf[N * N];
int suma[N * N], sumb[N * N], sumpa[N * N], sumpb[N * N]; struct Node {
int x, y, num;
friend bool operator <(const Node &p, const Node &q) {
return p.num < q.num;
}
}node[N * N]; inline int ksm (int a, int n) {
int ret = 1;
while (n) {
if (n % 2 == 1) ret = (ret * a) % mod;
a = (a * a) % mod;
n /= 2;
}
return ret;
} inline void init() {
for (int i = 1; i <= tot; i++) {
sumpa[i] = (sumpa[i - 1] + node[i].x * node[i].x) % mod;
sumpb[i] = (sumpb[i - 1] + node[i].y * node[i].y) % mod;
suma[i] = (suma[i - 1] + node[i].x) % mod;
sumb[i] = (sumb[i - 1] + node[i].y) % mod;
}
return ;
} inline void solve(int pos, int id) {
// std::cout << pos << " " << id << std::endl;
f[pos] = (f[pos] + sumpa[id] + sumpb[id]) % mod;
f[pos] = f[pos] - 2 * node[pos].x * suma[id] - 2 * node[pos].y * sumb[id];
f[pos] = (f[pos] % mod + mod) % mod;
int px = node[pos].x * node[pos].x;
int py = node[pos].y * node[pos].y;
f[pos] = (f[pos] + id * px + id * py) % mod;
f[pos] = (f[pos] + sumf[id]) % mod;
f[pos] *= ksm(id, mod - 2); f[pos] %= mod;
sumf[pos] = (sumf[pos - 1] + f[pos]) % mod;
return ;
} signed main(void) {
std::cin >> n >> m;
for (int i = 1; i <= n; i++)
for (int j = 1; j <= m; j++)
scanf("%d", &a[i][j]);
for (int i = 1; i <= n; i++) {
for (int j = 1; j <= m; j++) {
node[++tot].x = i;
node[tot].y = j;
node[tot].num = a[i][j];
}
}
std::cin >> x >> y;
std::sort(node + 1, node + 1 + tot);
init();
int last;
node[0].num = -114514;
for (int i = 1; i <= tot; i++) {
if (node[i].num != node[i - 1].num) last = i - 1;
solve(i, last);
if (node[i].x == x && node[i].y == y) {
std::cout << f[i] << std::endl;
return 0;
}
}
}

CF1042E Vasya and Magic Matrix 题解的更多相关文章

  1. CF1042E Vasya and Magic Matrix

    感觉不会期望. 首先把所有格子按照权值从小到大排一下序,这样一共有$n * m$个元素,每个元素有三个属性$x, y, val$. 下文中的下标均为排序后的下标. 这样子我们就可以推出公式: $f_i ...

  2. CF 1042 E. Vasya and Magic Matrix

    E. Vasya and Magic Matrix http://codeforces.com/contest/1042/problem/E 题意: 一个n*m的矩阵,每个位置有一个元素,给定一个起点 ...

  3. Vasya and Magic Matrix CodeForces - 1042E (概率dp)

    大意:给定n*m矩阵, 初始位置(r,c), 每一步随机移动到权值小于当前点的位置, 得分为移动距离的平方, 求得分期望. 直接暴力dp的话复杂度是O(n^4), 把距离平方拆开化简一下, 可以O(n ...

  4. Educational Codeforces Round 9 F. Magic Matrix 最小生成树

    F. Magic Matrix 题目连接: http://www.codeforces.com/contest/632/problem/F Description You're given a mat ...

  5. Educational Codeforces Round 48 (Rated for Div. 2) D 1016D Vasya And The Matrix (构造)

    D. Vasya And The Matrix time limit per test 2 seconds memory limit per test 256 megabytes input stan ...

  6. codeforces1016 D. Vasya And The Matrix(思维+神奇构造)

    D. Vasya And The Matrix time limit per test 2 seconds memory limit per test 256 megabytes input stan ...

  7. Codeforces 632F Magic Matrix(bitset)

    题目链接  Magic Matrix 考虑第三个条件,如果不符合的话说明$a[i][k] < a[i][j]$ 或 $a[j][k] < a[i][j]$ 于是我们把所有的$(a[i][j ...

  8. D. Vasya And The Matrix(Educational Codeforces Round 48)

    D. Vasya And The Matrix time limit per test2 seconds memory limit per test256 megabytes inputstandar ...

  9. Vasya And The Matrix CodeForces - 1016D (思维+构造)

    Now Vasya is taking an exam in mathematics. In order to get a good mark, Vasya needs to guess the ma ...

随机推荐

  1. 详解javascript的eventloop(一)

    先看一段代码 console.log('Hi') setTimeOut(function cb1(){ console.log('cb1') },500) console.log('Bye') 第一步 ...

  2. BUUCTF-Web:[GXYCTF2019]Ping Ping Ping

    题目 解题过程 1.题目页面提示?ip=,猜测是让我们把这个当做变量上传参数,由此猜想是命令注入 2.用管道符加上linux常用命令ls(windwos可以尝试dir)试试 所谓管道符(linux)的 ...

  3. 关于Mysql索引的数据结构

    索引的数据结构 1.为什么使用索引 概念: 索引是存储索引用于快速找到数据记录的一种数据结构,就好比一本书的目录部分,通过目录中对应的文章的页码,便可以快速定位到需要的文章,Mysql 中也是一样的道 ...

  4. CoreWCF 1.0 正式发布,支持 .NET Core 和 .NET 5+ 的 WCF

    CoreWCF 1.0 正式发布,支持 .NET Core 和 .NET 5+ 的 WCF https://devblogs.microsoft.com/dotnet/corewcf-v1-relea ...

  5. [题解] 春荔(cut) | 贪心

    题目大意 有一个长度为 \(n\) 的非负整数序列 \(a_i\),每次可以选择一段区间减去 \(1\),要求选择的区间长度 \(\in[l,r]\),问最少多少次把每个位置减成 \(0\). 不保证 ...

  6. 公司官网建站笔记(一):腾讯云服务器装CentOS8.2系统、重置密码、远程ssh登陆、sftp传递文件以及新建开发者账户

    前言   本篇使用的是腾讯云服务器,讲解了部署安装服务器CentOS8.2系统,重置密码,添加用户,远程登陆,远程传递文件等基本流程.   前提条件   购买了腾讯云服务器,如下图:     云服务器 ...

  7. git提交代码到GitHub操作-简易版(后续完善)

    一.git上传代码到GitHub 1.远程仓库GitHub创建好一个新仓库注意仓库名 2.本地建一个目录写代码,目录名与仓库命名一致 3.在目录下右键 git Bash here 打开git终端命令行 ...

  8. 用NetworkX生成并绘制(带权)无向图

    NetworkX是一个非常强大的网络科学工具,它封装了图的数据结构和许多经典图算法,也内置了许多可视化函数可供调用. 1. 随机图生成 最经典的随机图当属我们在上一篇博客<Erdos-Renyi ...

  9. 深入C++06:深入掌握OOP最强大的机制

    深入掌握OOP最强大的机制 1. 继承的基本意义 类与类之间的关系:①组合:a part of ... 一部分的关系:②继承: a kind of ... 属于同一种的关系: 继承的本质:a. 代码的 ...

  10. Golang:手撸一个支持六种级别的日志库

    Golang标准日志库提供的日志输出方法有Print.Fatal.Panic等,没有常见的Debug.Info.Error等日志级别,用起来不太顺手.这篇文章就来手撸一个自己的日志库,可以记录不同级别 ...