[HNOI2008]遥远的行星

嘟嘟嘟

 

此题我想了半个点，得出一个重要结论：我除了暴力以外啥也想不出来。

然后看了一下题解，彻底懵了……

对于每一个数 j，令x = j * a，因为题中说了一句“只要结果的相对误差不超过5%即可”，所以用分块解决这[1, x]。正好一块的部分[L, R]就都除以 j - (L +R) / 2……没错中位数……多出来的部分就暴力统计……【真德秀】。（块大小100左右就行）

#include<cstdio>
#include<iostream>
#include<cmath>
#include<algorithm>
#include<cstring>
#include<cstdlib>
#include<cctype>
#include<vector>
#include<stack>
#include<queue>
using namespace std;
#define enter printf("
")
#define space printf(" ")
#define Mem(a) memset(a, 0, sizeof(a))
typedef long long ll;
typedef double db;
const int INF = 0x3f3f3f3f;
const int eps = 1e-8;
const int maxn = 1e5 + 5;
inline ll read()
{
    ll ans = 0;
    char ch = getchar(), last = ' ';
    while(!isdigit(ch)) {last = ch; ch = getchar();}
    while(isdigit(ch))
    {
        ans = ans * 10 + ch - '0'; ch = getchar();
    }
    if(last == '-') ans = -ans;
    return ans;
}
inline void write(ll x)
{
    if(x < 0) x = -x, putchar('-');
    if(x >= 10) write(x / 10);
    putchar(x % 10 + '0');
}

const int T = 100;
int n;
db a;
ll m[maxn], sum[maxn];

int main()
{
    n = read(); scanf("%lf", &a);
    for(int i = 1; i <= n; ++i) {m[i] = read(); sum[i] = sum[i - 1] + m[i];}
    for(int i = 1; i <= n; ++i)
    {
        int x = i * a; 
        db ans = 0;
        int L = 1, R = T;
        while(R <= x)    //在一块 
        {
            ans += (sum[R] - sum[L - 1]) * m[i] / (i - ((L + R) >> 1));
            L += T; R += T;
        }
        for(int j = L; j <= x; ++j) ans += m[j] * m[i] / (i - j);    //块外暴力统计 
        printf("%.6lf
", ans);
    }
    return 0;
}