Solution -「ARC 126F」Affine Sort
\(\mathcal{Description}\)
Link.
给定 \(\{x_n\}\),令
\]
求出
\]
\(n\le10^3\),\(\sum x\le5\times10^5\)。
\(\mathcal{Solution}\)
太炫酷了叭。
令 \(g(k)=f(k)-f(k-1)\),断言这样一条性质:
\]
甚至可以认为给原式洛了一发。
证明
假定存在 $c=\lim_{k\rightarrow+\infty}\frac{g(k)}{k^2}$,则对于任意 $\varepsilon\in\mathbb R_+$,应有足够大的 $k$,使得 $(c-\eps)k^2\le g(k)\le (c+\eps)k^2$。继而,也应有足够大的 $k$ 使得 $(c-\eps)\sum_{i\le k}i^2\le f(k)\le(c+\eps)\sum_{i\le k}i^2$。考虑到 $\sum_{i\le k}i^2=\frac{k^3}{3}+\mathcal O(k^2)$,可以得到 $\frac{1}{3}c-\eps\le \frac{f(k)}{k^3}\le \frac{1}{3}c+\eps$。所以基于这个假设,得到 $\lim_{k\rightarrow+\infty}\frac{f(k)}{k^3}=\frac{1}{3}\lim_{k\rightarrow+\infty}\frac{g(k)}{k^2}$。
此后,考虑使得 \((ax_i+b)\bmod k\) 单增的 \((a,b)\),根据取模的性质,应该有 \(\left\{\frac{ax_i+b}{k}\right\}\) 单增,其中 \(\{x\}=x-\lfloor x\rfloor\) 表示取非负数 \(x\) 的小数部分。所以,令 \(\alpha=\frac{a}{k},\beta=\frac{b}{k}\),我们仅需研究使得 \(\{\alpha x_i+\beta\}\) 单增的 \((\alpha,\beta)\in[0,1)^2\)。具体地,令
\]
我们所求即为 \(D\) 的面积。
这步转化的具体数学原理我并没有学懂,这貌似也涉及了上文中“假定极限存在”的证明。当然存在一些感性理解方法也显得不够具有说服力。如果希望深刻研究请参考 官方题解 叭,这里先跳过,之后的步骤仍有思考价值。
我们把 \(\{y_i=\{\alpha x_i+\beta\}\}\) 放在 \([0,1)\) 这个关于 \(1\) 的“模域”上(Tiw: 你把“模”和“域”放一起是完全不合理的!我:?),进一步研究 \((\alpha,\beta)\) 的性质:
可见,\(\alpha\) 能决定 \(y_i\) 在圆周上的相对位置关系,而 \(\beta\) 仅能整体旋转 \(y\),也即是确定 \(0\) 的位置。显然,对于使得相对位置合法的 \(\alpha\),合法 \(\beta\) 取值区间的长度为 \(\{\alpha(x_1-x_n)\}\);另一方面,对于这样的 \(\alpha\),应有
\]
同时,\(\alpha\in[0,1)\),\(x_i\in\mathbb N\),则 \(\{\alpha x_{i+1}\}-\{\alpha x_i\}\equiv\{\alpha(x_{i+1}-x_i)\}\pmod 1\),可见在 \(\alpha\) 从 \(0\) 向 \(1\) 取值时,仅有经过某个 \(\frac{t}{|x_{i+1}-x_i|}\) 会改变左式的值,即这些点把 \([0,1)\) 分隔为若干区间,若一个区间合法,再积分求此时 \(\beta\) 的贡献,即 \(\int_{l}^r\{\alpha(x_1-x_n)\}\text d\alpha\)。注意 \(|x_n-x_{n+1}|\) 也在分隔点的分母中,所以这是单纯的一次函数积分。
最终,复杂度瓶颈为取出 \(\mathcal O(\sum x)\) 个区间排序,即 \(\mathcal O(\sum x\log\sum x)\)。
\(\mathcal{Code}\)
/*~Rainybunny~*/
#include <bits/stdc++.h>
#define rep( i, l, r ) for ( int i = l, rep##i = r; i <= rep##i; ++i )
#define per( i, r, l ) for ( int i = r, per##i = l; i >= per##i; --i )
const int MAXN = 1e3, MOD = 998244353, INV6 = 166374059;
int n, x[MAXN + 5];
inline int iabs( const int u ) { return u < 0 ? -u : u; }
inline int mul( const int u, const int v ) { return 1ll * u * v % MOD; }
inline int sub( int u, const int v ) { return ( u -= v ) < 0 ? u + MOD : u; }
inline int add( int u, const int v ) { return ( u += v ) < MOD ? u : u - MOD; }
inline int mpow( int u, int v = MOD - 2 ) {
int ret = 1;
for ( ; v; u = mul( u, u ), v >>= 1 ) ret = mul( ret, v & 1 ? u : 1 );
return ret;
}
struct Atom {
int a, b, c;
inline bool operator < ( const Atom& t ) const {
return 1ll * a * t.b < 1ll * t.a * b;
}
};
std::vector<Atom> sec;
int main() {
scanf( "%d", &n );
rep ( i, 1, n ) scanf( "%d", &x[i] );
x[n + 1] = x[1];
int circ = 0;
rep ( i, 1, n ) {
int k = iabs( x[i] - x[i + 1] ), d = x[i] < x[i + 1] ? -1 : 1;
rep ( j, 1, k - 1 ) sec.push_back( { j, k, d } );
circ += x[i] >= x[i + 1];
}
sec.push_back( { 0, 1, 0 } ), sec.push_back( { 1, 1, 0 } );
std::sort( sec.begin(), sec.end() );
int ans = 0;
for ( size_t i = 1; i < sec.size(); ++i ) {
if ( ( circ += sec[i - 1].c ) != 1 ) continue;
int dn = sub( mul( sec[i].a, mpow( sec[i].b ) ),
mul( sec[i - 1].a, mpow( sec[i - 1].b ) ) );
int le = ( 1ll * sec[i - 1].a * ( x[1] - x[n] ) % sec[i - 1].b
+ sec[i - 1].b ) % sec[i - 1].b;
if ( !le && x[1] < x[n] ) le = sec[i - 1].b;
le = mul( le, mpow( sec[i - 1].b ) );
int ri = ( 1ll * sec[i].a * ( x[1] - x[n] ) % sec[i].b
+ sec[i].b ) % sec[i].b;
if ( !ri && x[1] > x[n] ) ri = sec[i].b;
ri = mul( ri, mpow( sec[i].b ) );
ans = add( ans, mul( add( le, ri ), dn ) );
}
printf( "%d\n", mul( ans, INV6 ) );
return 0;
}
Solution -「ARC 126F」Affine Sort的更多相关文章
- Solution -「ARC 104E」Random LIS
\(\mathcal{Description}\) Link. 给定整数序列 \(\{a_n\}\),对于整数序列 \(\{b_n\}\),\(b_i\) 在 \([1,a_i]\) 中等概率 ...
- Solution -「ARC 101D」「AT4353」Robots and Exits
\(\mathcal{Description}\) Link. 有 \(n\) 个小球,坐标为 \(x_{1..n}\):还有 \(m\) 个洞,坐标为 \(y_{1..m}\),保证上述坐标 ...
- Solution -「ARC 125F」Tree Degree Subset Sum
\(\mathcal{Description}\) Link. 给定含有 \(n\) 个结点的树,求非负整数对 \((x,y)\) 的数量,满足存在 \(\exist S\subseteq V ...
- Solution -「ARC 125E」Snack
\(\mathcal{Description}\) Link. 把 \(n\) 种零食分给 \(m\) 个人,第 \(i\) 种零食有 \(a_i\) 个:第 \(i\) 个人得到同种零食数量 ...
- Solution -「ARC 110D」Binomial Coefficient is Fun
\(\mathcal{Description}\) Link. 给定非负整数序列 \(\{a_n\}\),设 \(\{b_n\}\) 是一个非负整数序列且 \(\sum_{i=1}^nb_i\ ...
- Solution -「ARC 063D」「AT 2149」Snuke's Coloring 2
\(\mathcal{Decription}\) Link. 平面上有一个左下角坐标 \((0,0)\) 右上角坐标 \((W,H)\) 的矩形,起初长方形内部被涂白. 现在给定 \(n\) ...
- Solution -「ARC 124E」Pass to Next
\(\mathcal{Description}\) Link. 有 \(n\) 个人站成一个环,初始时第 \(i\) 个人手里有 \(a_i\) 个球.第 \(i\) 个人可以将自己手中任意数 ...
- Solution -「ARC 126E」Infinite Operations
\(\mathcal{Description}\) Link. 给定序列 \(\{a_n\}\),定义一次操作为: 选择 \(a_i<a_j\),以及一个 \(x\in\mathbb R ...
- Solution -「ARC 058C」「AT 1975」Iroha and Haiku
\(\mathcal{Description}\) Link. 称一个正整数序列为"俳(pái)句",当且仅当序列中存在连续一段和为 \(x\),紧接着连续一段和为 \(y ...
随机推荐
- js 使用 substr 和substring 的区别记录
根据 提示,已经很清楚区别了 区别都是第二个参数 ,substr 是长度 ,而substring 仍是索引数
- 实验 1 :Mininet 源码安装和可视化拓扑
实验 1 : Mininet 源码安装和可视化拓扑工具 一 .实验目的 掌握 Mininet 的源码安装方法和 miniedit 可视化拓扑生成工具. 二 .实验任务 使用源码安装 Mininet 的 ...
- [转]webpack配置本地服务器
亲测,webpack打包vue项目之后生成的dist文件可以部署到 express 服务器上运行. 我的vue项目结构如下: 1. 进入该vue项目目录,打开git bash,执行:npm run b ...
- Typora中本地图片无法上传CSDN解决方案
解决方法 本地图片无法上传,我们可以选择使用在线免费图床 把想要使用的图先上传到图床后复制对应的MarkDown语句到typora即可 在这里,给大家推荐一个图床 [图床链接](Image Uploa ...
- 开发 IDEA Plugin 引入探针,基于字节码插桩获取执行SQL
作者:小傅哥 博客:https://bugstack.cn 沉淀.分享.成长,让自己和他人都能有所收获! 一.前言 片面了! 一月三舟,托尔斯泰说:"多么伟大的作家,也不过就是在书写自己的片 ...
- C++初始化列表各情况分析
今天回顾了下C++初始化列表的知识,接下来我对这一知识作一总结. 我们在定义了一个类的时候,需要对类的成员进行初始化.关于初始化,有两种方法,一种在初始化列表中进行,另一种就是在构造函数中进行,对于这 ...
- linux设置系统用户密码
目录 一:系统用户密码 1.设置用户密码 一:系统用户密码 1.设置用户密码 1.交互式方法 passwd [用户名] 2.免交互式 echo [设置密码] | passwd --stdin [用户名 ...
- python 小兵(8)闭包和装饰器
闭包"是什么,以及,更重要的是,写"闭包"有什么用处. (个人理解) 1."闭包"是什么 首先给出闭包函数的必要条件: 闭包函数必须返回一个函数对象 ...
- 用Json给表单赋值
$.extend({ setForm :function(frm,jsonValue) { var obj=$(frm); $.each(jsonValue, function (name, ival ...
- hihoCoder挑战赛1 毁灭者问题
题目链接:http://hihocoder.com/problemset/problem/1034 数据结构题,由于每个魔法单位有着不同的回复速度和上限,所以不能根据吸收时间点进行查询和更新.但是如果 ...