Description 给出n个数qi,给出Fj的定义如下: $$F_j = \sum_{i<j}\frac{q_i q_j}{(i-j)^2 }-\sum_{i>j}\frac{q_i q_j}{(i-j)^2 }$$ 令Ei=Fi/qi,求Ei. Input 第一行一个整数n. 接下来n行每行输入一个数,第i行表示qi. n≤100000,0<qi<1000000000 Output n行,第i行输出Ei.与标准答案误差不超过1e-2即可. Sample Input 5 400…

题面 事实说明只会FFT板子是没有用的,还要把式子推成能用FFT/转化一下卷积的方式 虽然这个题不算难的多项式卷积 稍微化简一下可以发现实际是$q_i$和$\frac{1}{(i-j)^2}$在卷,然后每两项是在向下标差值的那项做贡献,而直接卷是向两项下标和的那项做贡献.于是把前半部分的$\frac{1}{(i-j)^2}$做成负的,后半段的做成正的,这样卷完后半段就是题目要求的东西.当然把一个序列反过来再卷也是对的 #include<cmath> #include<cstdio>…

代换一下变成多项式卷积,这里是的答案是两个卷积相减,FFT求一下两个卷积就可以啦 详细的题解:http://www.cnblogs.com/iwtwiioi/p/4126284.html #include<cmath> #include<cstdio> #include<cstring> #include<algorithm> using namespace std; const int N = 500003; const double Pi = acos(…

Problem Description 给出 \(n\) 个数 \(q_i\),给出 \(F_j\) 的定义如下: \[F_j=\sum_{i<j} \frac{q_iq_j}{(i-j)^2} - \sum_{i>j} \frac{q_iq_j}{(i-j)^2}\] 令 \(E_i=F_i/q_i\),求 \(E_i\). Input Format 第一行一个整数\(n\). 接下来 \(n\) 行每行输入一个数,第 \(i\) 行表示 \(q_i\). Output Format \(n…

题意 https://loj.ac/problem/2201 思路 说白了就是一条路径上有 \(n\) 个二维坐标,求一条直线使得所有点到此直线的距离和最小. 设这条直线为 \(y=kx+b\) ,距离和为 \(\delta\). \[ \delta=\sum{(kx_i-y_i+b)^2\over k^2+1} \] \[ \delta={k^2\sum x_i^2+\sum y_i^2+nb^2-2k\sum x_iy_i+2kb\sum x_i-2b\sum y_i\over k^2+1}…

php大力力 [050节] 兄弟连高洛峰 PHP教程 2014年[数据库.PDO教程] 第14章 数据库252.[2014]兄弟连高洛峰 PHP教程14.1.1 复习数据库[已发布,点击下载]253.[2014]兄弟连高洛峰 PHP教程14.1.2 phpMyAdmin的使用[已发布,点击下载]254.[2014]兄弟连高洛峰 PHP教程14.1.3 php访问MySQL[已发布,点击下载]255.[2014]兄弟连高洛峰 PHP教程14.1.4 在PHP脚本中操作MySQL数据库1[已发布,点…

2015-08-25 php大力力016 兄弟连高洛峰php教程(2014年 14章数据库章节列表) [2014]兄弟连高洛峰 PHP教程14.1.1 复习数据库  15:58 [2014]兄弟连高洛峰 PHP教程14.1.2 phpMyAdmin的使用 15:59 [2014]兄弟连高洛峰 PHP教程14.1.3 php访问MySQL 17:27 [2014]兄弟连高洛峰 PHP教程14.1.4 在PHP脚本中操作MySQL数据库1  17:38 [2014]兄弟连高洛峰 PHP教程14.1.…

http://www.verycd.com/topics/2843130/ 第1部分 WEB开发入门篇第1章LAMP网站构建1.[2014]兄弟连高洛峰 PHP教程1.1.1 新版视频形式介绍[已发布,点击下载] 2.[2014]兄弟连高洛峰 PHP教程1.1.2 BS结构软件类型介绍[已发布,点击下载] 3.[2014]兄弟连高洛峰 PHP教程1.1.3 现在是Web2.0的时代[已发布,点击下载] 4.[2014][2014]兄弟连高洛峰 PHP教程1.1.4 Web开发标准[已发布,点击下…

无题面神题 原题意: 求所有的Ei=Fi/qi. 题解: qi被除掉了,则原式中的qj可以忽略. 用a[i]表示q[i],用b[j-i]来表示±1/((j-i)^2)(j>i时为正,j<i时为负) 则求E[j]就是多项式乘法了. 因为是FFT,所以b的下标要增加到0及以上. 这题时限有30s,比某题友好多了. 代码: type xs=record x,y:double; end; arr=..]of xs; var e,t:arr; a:..]of arr; n,m,i:longint; fu…

显然fft维护卷积就可以了 发现fft里面会改变很多东西 要还原一下 #include <bits/stdc++.h> #define dob complex<double> using namespace std; const int N=3e5; const double pi=acos(-1.0); dob a[N],a2[N],b[N]; int r[N],l; double ans1[N],sum[N]; int n,m; void fft(dob *a,int o) {…