Math

题目

题目描述
构造一个长度为 $n$ 的序列 $a_1,a_2,a_3,\dots ,a_n$ ，使得其满足如下条件。

• $10\le n\le 10^3$
• $\forall x\in [1,n],|a_x|\le 10^9$
• $\mu =\frac{{\sum }_{x=1}^n{a}_x}{n}$ ，即平均数为 $\mu$
• $\nu =\frac{{\sum }_{x=1}^n(a_x-\mu {)}^2}{n}$ ，即方差为 $\nu$

输出任意一组解即可。

数据范围与约定
数据组数 $T\le 10,\mu \le 10^6,0\le \nu \le 10^9$

思路

首先将所有的 $a$ 减去一个 $\mu$ ，可以简化一下题意：

• 一次方之和为零。即 ${\sum }_{x=1}^n{a}_x=0$ 。
• 二次方之和为 $n\nu$ 。即 ${\sum }_{x=1}^n{a}_x^2=n\nu$ 。

为了满足第一条，我们规定 $2|n$ ，且 $\forall k\in [1,\frac{n}{2}],a_{2k-1}+a_{2k}=0$ 。

此时，我们只需要求一个长度为 $\frac{n}{2}$ 的数列 $⟨b_i⟩$，并且满足 ${\sum }_{x=1}^{n/2}{b}_x^2=\frac{n\nu }{2}$

这并不难，只需要每次将 $b_x$ 取的尽量大即可。因为最后一定有解（只要取 $b_x=1$ 就行了）。

$n$ 可以任取一个足够大的偶数，取 $n=10^3$ 是最好的。

代码

#include <cstdio>
#include <iostream>
#include <vector>
#include <algorithm>
#include <cmath>
using namespace std;
inline int readint(){
	int a = 0; char c = getchar(), f = 1;
	for(; c<'0' or c>'9'; c=getchar())
		if(c == '-') f = -f;
	for(; '0'<=c and c<='9'; c=getchar())
		a = (a<<3)+(a<<1)+(c^48);
	return a*f;
}
inline void writeint(long long x){
	if(x < 0) putchar('-'), x = -x;
	if(x > 9) writeint(x/10);
	putchar((x%10)^48);
}
# define MB template < class T >
MB void getMax(T &a,const T &b){ if(a < b) a = b; }
MB void getMin(T &a,const T &b){ if(b < a) a = b; }

const int MaxN = 1005;
int a[MaxN];

int main(){
	freopen("math.in","r",stdin);
	freopen("math.out","w",stdout);
	for(int T=readint(); T; --T){
		int u = readint();
		long long v = readint();
		int n = 1000; v *= n;
		for(int i=1; i<=(n>>1); ++i){
			a[i] = int(sqrt(v>>1));
			v -= 2ll*a[i]*a[i];
		}
		printf("%d\n",n);
		for(int i=1; i<=(n>>1); ++i)
			printf("%d %d ",a[i]+u,-a[i]+u);
		putchar('\n');
	}
	return 0;
}