bzoj1011

因为允许5%的误差。。所以把＞一定长度的一段看成一段近似计算就行了。。

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<ctime>
#include<cmath>
#include<iostream>
#include<algorithm>
#include<queue>
#include<stack>
#include<set>
#define clr(a,x) memset(a,x,sizeof(a))
#define rep(i,l,r) for(int i=(l);i<(r);i++)
using namespace std;
typedef long long ll;
typedef pair<int,int> pii;
#define mkp(a,b) make_pair(a,b)
int read(){
    int ans=0,f=1;
    char c=getchar();
    while(!isdigit(c)){
        if(c=='-') f=-1;
        c=getchar();
    }
    while(isdigit(c)){
        ans=ans*10+c-'0';
        c=getchar();
    }
    return ans*f;
}
const int maxn=100009;
int n;
ll w[maxn];
double a;
void solve(int k,int l,int r){
    double ans=0;
    if(r-l<=500){
        rep(i,l,r+1) ans+=1.0*(w[i]-w[i-1])/(k-i);
    }else{
        while(l<=r){
            int t=min(l+500,r);
            ans+=1.0*(w[t]-w[l-1])/(k-(l+t)/2);
            l=t+1;
        }
    }
    printf("%.7lf\n",ans*(w[k]-w[k-1]));
}
int main(){
    n=read();scanf("%lf",&a);
    rep(i,1,n+1) w[i]=w[i-1]+read();
    rep(i,1,n+1) solve(i,1,floor(a*i));
    return 0;
}

1011: [HNOI2008]遥远的行星

Time Limit: 10 Sec  Memory Limit: 162 MBSec  Special Judge
Submit: 2872  Solved: 1034
[Submit][Status][Discuss]

Description

直线上N颗行星，X=i处有行星i,行星J受到行星I的作用力，当且仅当i<=AJ.此时J受到作用力的大小为 Fi->j=Mi*Mj/(j-i) 其中A为很小的常量，故直观上说每颗行星都只受到距离遥远的行星的作用。请计算每颗行星的受力，只要结果的相对误差不超过5%即可.

Input

第一行两个整数N和A. 1<=N<=10^5.0.01< a < =0.35 
接下来N行输入N个行星的质量Mi,保证0<=Mi<=10^7

Output

N行，依次输出各行星的受力情况

Sample Input

5 0.3
3
5
6
2
4

Sample Output

0.000000
0.000000
0.000000
1.968750
2.976000

HINT

 

精确结果应该为0 0 0 2 3,但样例输出的结果误差不超过5%,也算对

 

Source

[Submit][Status][Discuss]