题面

题解

恩, 我们首先有这两个关系

\[\displaystyle\begin{aligned} F_j &= \sum_{i < j}\frac{q_iq_j}{(i - j)^2} - \sum_{i > j}\frac{q_iq_j}{(i - j)^2}\\ &= q_j\cdot(\sum_{i < j}\frac{q_i}{(i - j)^2}-\sum_{i>j}\frac{q_i}{(i - j)^2})\end{aligned}
\]

然后由

\[\displaystyle E_i = \frac{F_i}{q_i}
\]

我们可以推出这样的式子

\[\displaystyle\begin{aligned}
E_i&=\frac{F_i}{q_i}\\&=\sum_{j<i}\frac{q_j}{(i - j)^2}-\sum_{j>i}\frac{q_j}{(j - i)^2}\\
\end{aligned}
\]

我们来看一下 n = 3 的情况

我们将红色的线的两端乘起来, 再把三条红线的和相加就得到了\(E_0\), 同理我们通过蓝线, 绿线可以得到\(E_1\), \(E_2\)

我们定义\(a_i = p_i\), \(b_i = \frac{1}{i ^ 2}\), 特别的, 我们定义\(b_{-i} = -\frac{1}{i ^ 2}\), 先别管为什么下标为负

\[\displaystyle\begin{aligned}
E_0&=a_0b_0+a_1b_{-1}+a_2b_{-2}\\
E_1&=a_0b_1+a_1b_0+a_2b_{-1}\\
E_2&=a_0b_2+a_1b_1+a_2b_0
\end{aligned}
\]

恩, 我们来讨论数组下标小于零的问题, 整体平移就行了嘛, 注意到\(b\)最小的下标为\((-3 + 1)\), 所以我们平移\((3 - 1)\)位就可以了, 所以有:

\[\begin{aligned}
E_0&=a_0b_2+a_1b_1+a_2b_0\Leftrightarrow E_0 = \sum_{i+j=2}a_ib_j\\
E_1&=a_0b_3+a_1b_2+a_2b_1\Leftrightarrow E_1 = \sum_{i+j=3}a_ib_j\\
E_2&=a_0b_4+a_1b_3+a_2b_2\Leftrightarrow E_2 = \sum_{i+j=4}a_ib_j
\end{aligned}
\]

恩, 我们来看一下卷积的形式

\[c_k = \sum_{i + j = k}a_ib_j
\]

我们会发现这里的\(E\)数组的下标离可以卷积差了一个\((3 - 1)\), 用\(E_2\)代替\(E_0\)就可以用\(a\)与\(b\)的卷积来算了

我们将特殊情况扩展到一般情况, 设\(E\)共有\(n\)项, 自己瞎猜一下可以发现\(b\)数组下标平移了\(n - 1\)位, \(E\)数组的下标平移了\(n - 1\)位, 所以, 最后我们输出\(E\)的第\([n - 1, (n - 1) + n - 1]\)项即可

这篇题解就是写给我们这种蒟蒻看的, 至于卷积, 我猜各位是为了练习FFT才来做这道题的吧

Code

#include <algorithm>
#include <iostream>
#include <cstring>
#include <cstdlib>
#include <cstdio>
#include <vector>
#include <cmath>
#define itn int
#define reaD read
#define N 1000005
const double pi = acos(-1);
using namespace std; int n, m, cnt, rev[N];
struct complex
{
double x, y;
complex(double xx = 0, double yy = 0) { x = xx; y = yy; }
complex operator + (complex p) { return complex(x + p.x, y + p.y); }
complex operator - (complex p) { return complex(x - p.x, y - p.y); }
complex operator * (complex p) { return complex(x * p.x - y * p.y, x * p.y + y * p.x); }
} a[N], b[N]; inline int read()
{
int x = 0, w = 1; char c = getchar();
while(c < '0' || c > '9') { if (c == '-') w = -1; c = getchar(); }
while(c >= '0' && c <= '9') { x = x * 10 + c - '0'; c = getchar(); }
return x * w;
} void FFT(complex* p, int opt)
{
for(int i = 0; i < m; i++) if(i < rev[i]) swap(p[i], p[rev[i]]);
for(int i = 1; i < m; i <<= 1)
{
complex cur = complex(cos(pi / i), opt * sin(pi / i));
for(int j = 0; j < m; j += (i << 1))
{
complex w = complex(1, 0);
for(int k = 0; k < i; k++, w = w * cur)
{
complex l = p[j + k], r = w * p[i + j + k];
p[j + k] = l + r; p[i + j + k] = l - r;
}
}
}
} int main()
{
n = read();
for(int i = 0; i < n; i++) scanf("%lf", &a[i].x);
for(int i = 1; i < n; i++)
{
b[i + n - 1].x = 1 / (1.0 * (double) i * i);
b[n - 1 - i].x = -b[i + n - 1].x;
}
for(m = 1; m <= 3 * n - 2; m <<= 1, cnt++);
for(int i = 0; i < m; i++) rev[i] = (rev[i >> 1] >> 1) | ((i & 1) << (cnt - 1));
FFT(a, 1); FFT(b, 1);
for(int i = 0; i <= m; i++) a[i] = a[i] * b[i];
FFT(a, -1);
for(int i = n - 1; i <= 2 * n - 2; i++) printf("%.3lf\n", a[i].x / m);
return 0;
}

这可能是我写得最认真的一篇题解吧, 希望看这篇题解的泥萌可以康懂

[题解] [ZJOI2014] 力的更多相关文章

  1. [ZJOI3527][Zjoi2014]力

    [ZJOI3527][Zjoi2014]力 试题描述 给出n个数qi,给出Fj的定义如下: 令Ei=Fi/qi.试求Ei. 输入 包含一个整数n,接下来n行每行输入一个数,第i行表示qi. 输出 有n ...

  2. [洛谷P3338] [ZJOI2014]力

    洛谷题目链接:P3338 [ZJOI2014]力 题目描述 给出n个数qi,给出Fj的定义如下: \[F_j = \sum_{i<j}\frac{q_i q_j}{(i-j)^2 }-\sum_ ...

  3. bzoj3527: [Zjoi2014]力 fft

    bzoj3527: [Zjoi2014]力 fft 链接 bzoj 思路 但是我们求得是 \(\sum\limits _{i<j} \frac{q_i}{(i-j)^2}-\sum_{i> ...

  4. 洛谷 P3338 [ZJOI2014]力 解题报告

    P3338 [ZJOI2014]力 题目描述 给出n个数qi,给出Fj的定义如下: \(F_j = \sum_{i<j}\frac{q_i q_j}{(i-j)^2 }-\sum_{i>j ...

  5. 【BZOJ 3527】 3527: [Zjoi2014]力 (FFT)

    3527: [Zjoi2014]力 Time Limit: 30 Sec  Memory Limit: 256 MBSec  Special JudgeSubmit: 2003  Solved: 11 ...

  6. P3338 [ZJOI2014]力(FFT)

    题目 P3338 [ZJOI2014]力 做法 普通卷积形式为:\(c_k=\sum\limits_{i=1}^ka_ib_{k-i}\) 其实一般我们都是用\(i=0\)开始的,但这题比较特殊,忽略 ...

  7. [Luogu P3338] [ZJOI2014]力 (数论 FFT 卷积)

    题面 传送门: 洛咕 BZOJ Solution 写到脑壳疼,我好菜啊 我们来颓柿子吧 \(F_j=\sum_{i<j}\frac{q_i*q_j}{(i-j)^2}-\sum_{i>j} ...

  8. 笔记-[ZJOI2014]力

    [ZJOI2014]力 \[\begin{split} E_j=&\sum_{i=1}^{j-1}\frac{q_i}{(i-j)^2}-\sum_{i=j+1}^{n}\frac{q_i}{ ...

  9. [ZJOI2014]力 题解

    题目地址 洛谷P3338 Solution 第一道FFT的应用AC祭! 我们要求: \[E_j=\frac{F_j}{q_j}=\sum_{i<j}\frac{q_i}{(i-j)^2}-\su ...

随机推荐

  1. Where is __dso_handle defined?

    Where is __dso_handle defined? 来源  https://stackoverflow.com/questions/34308720/where-is-dso-handle- ...

  2. PHP之50个开源项目

    GitHub上50个最受欢迎的PHP开源项目[2019] 1.Laravel Laravel是一个为Web开发者打造的PHP开发框架. GitHub Stars: 43.5k+ 网址: https:/ ...

  3. 公众平台第三方平台 .NET开发

    前言:本博客借鉴了很多三方内容整理的,参考博客:竹叶苿. 一.开发的目的(以下是引用官方的话) 公众平台第三方平台 是为了让公众号或小程序运营者,在面向垂直行业需求时,可以一键授权给第三方平台(并且可 ...

  4. MySQL 5.7.18 zip版本的安装使用方法

    转自:https://www.cnblogs.com/nepulgh/p/7152618.html MySQL 5.7.18 zip版本的安装使用方法 这个版本的MySQL不像那种点击就可以立即安装, ...

  5. MySQL表内更新时,自动记录时间

    1.创建表: create table test_time(id int primary key not null,status  varchar(24),create_time datetime d ...

  6. SpringBoot static修饰的字段/方法如何获取application.yml配置

    SpringBoot的application.yml一种特殊的应用场景,一般我们获取application.yml的配置文件只要@Value就可以获取到值了,但是如果是static修饰的字段肯定就不能 ...

  7. C语言求π的方法

    #include <stdio.h> #include <math.h> int main() { int r; double PI,s; scanf("%d&quo ...

  8. qunee 缩略图

    <!DOCTYPE html> <html> <head> <title>Hello Qunee for HTML5</title> < ...

  9. Pyspark常用API总结

    DF 类似于二维表的数据结果 mame age 狗山石 23 获取df的列名: df.columns 显示当前值 打印 df.show() show(2) show括号里面传入参数可以显示查看几行 s ...

  10. 9.RNN应用

    import numpy as np from keras.datasets import mnist from keras.utils import np_utils from keras.mode ...