CF1042E Vasya and Magic Matrix 题解
思路分析
看到题目中 \(n,m \leq 1000\) ,故直接考虑 \(O(n^2)\) 级别做法。
我们先把所有的点按照 \(val\) 值从小到大排序,这样的话二维问题变成序列问题。
设 \(f_i\) 表示走到第 \(i\) 个点的价值的期望。
先列出裸的 \(dp\) 方程:(\(Num\) 表示符合条件的点的个数)
\]
但是这个好像是 \(O(n^2m^2)\) 的优秀算法……
不要担心,我们把式子化简一下:
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 题解的更多相关文章
- CF1042E Vasya and Magic Matrix
感觉不会期望. 首先把所有格子按照权值从小到大排一下序,这样一共有$n * m$个元素,每个元素有三个属性$x, y, val$. 下文中的下标均为排序后的下标. 这样子我们就可以推出公式: $f_i ...
- CF 1042 E. Vasya and Magic Matrix
E. Vasya and Magic Matrix http://codeforces.com/contest/1042/problem/E 题意: 一个n*m的矩阵,每个位置有一个元素,给定一个起点 ...
- Vasya and Magic Matrix CodeForces - 1042E (概率dp)
大意:给定n*m矩阵, 初始位置(r,c), 每一步随机移动到权值小于当前点的位置, 得分为移动距离的平方, 求得分期望. 直接暴力dp的话复杂度是O(n^4), 把距离平方拆开化简一下, 可以O(n ...
- Educational Codeforces Round 9 F. Magic Matrix 最小生成树
F. Magic Matrix 题目连接: http://www.codeforces.com/contest/632/problem/F Description You're given a mat ...
- 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 ...
- 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 ...
- Codeforces 632F Magic Matrix(bitset)
题目链接 Magic Matrix 考虑第三个条件,如果不符合的话说明$a[i][k] < a[i][j]$ 或 $a[j][k] < a[i][j]$ 于是我们把所有的$(a[i][j ...
- 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 ...
- 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 ...
随机推荐
- 详解javascript的eventloop(一)
先看一段代码 console.log('Hi') setTimeOut(function cb1(){ console.log('cb1') },500) console.log('Bye') 第一步 ...
- BUUCTF-Web:[GXYCTF2019]Ping Ping Ping
题目 解题过程 1.题目页面提示?ip=,猜测是让我们把这个当做变量上传参数,由此猜想是命令注入 2.用管道符加上linux常用命令ls(windwos可以尝试dir)试试 所谓管道符(linux)的 ...
- 关于Mysql索引的数据结构
索引的数据结构 1.为什么使用索引 概念: 索引是存储索引用于快速找到数据记录的一种数据结构,就好比一本书的目录部分,通过目录中对应的文章的页码,便可以快速定位到需要的文章,Mysql 中也是一样的道 ...
- CoreWCF 1.0 正式发布,支持 .NET Core 和 .NET 5+ 的 WCF
CoreWCF 1.0 正式发布,支持 .NET Core 和 .NET 5+ 的 WCF https://devblogs.microsoft.com/dotnet/corewcf-v1-relea ...
- [题解] 春荔(cut) | 贪心
题目大意 有一个长度为 \(n\) 的非负整数序列 \(a_i\),每次可以选择一段区间减去 \(1\),要求选择的区间长度 \(\in[l,r]\),问最少多少次把每个位置减成 \(0\). 不保证 ...
- 公司官网建站笔记(一):腾讯云服务器装CentOS8.2系统、重置密码、远程ssh登陆、sftp传递文件以及新建开发者账户
前言 本篇使用的是腾讯云服务器,讲解了部署安装服务器CentOS8.2系统,重置密码,添加用户,远程登陆,远程传递文件等基本流程. 前提条件 购买了腾讯云服务器,如下图: 云服务器 ...
- git提交代码到GitHub操作-简易版(后续完善)
一.git上传代码到GitHub 1.远程仓库GitHub创建好一个新仓库注意仓库名 2.本地建一个目录写代码,目录名与仓库命名一致 3.在目录下右键 git Bash here 打开git终端命令行 ...
- 用NetworkX生成并绘制(带权)无向图
NetworkX是一个非常强大的网络科学工具,它封装了图的数据结构和许多经典图算法,也内置了许多可视化函数可供调用. 1. 随机图生成 最经典的随机图当属我们在上一篇博客<Erdos-Renyi ...
- 深入C++06:深入掌握OOP最强大的机制
深入掌握OOP最强大的机制 1. 继承的基本意义 类与类之间的关系:①组合:a part of ... 一部分的关系:②继承: a kind of ... 属于同一种的关系: 继承的本质:a. 代码的 ...
- Golang:手撸一个支持六种级别的日志库
Golang标准日志库提供的日志输出方法有Print.Fatal.Panic等,没有常见的Debug.Info.Error等日志级别,用起来不太顺手.这篇文章就来手撸一个自己的日志库,可以记录不同级别 ...