因為允許5%的誤差。。所以把>一定長度的一段看成一段近似計算就行了。。


1 #include<cstdio> 2 #include<cstdlib> 3 #include<cstring> 4 #include<ctime> 5 #include<cmath> 6 #include<iostream> 7 #include<algorithm> 8 #include<queue> 9 #include<stack> 10 #include<set> 11 #define clr(a,x) memset(a,x,sizeof(a)) 12 #define rep(i,l,r) for(int i=(l);i<(r);i++) 13 using namespace std; 14 typedef long long ll; 15 typedef pair<int,int> pii; 16 #define mkp(a,b) make_pair(a,b) 17 int read(){ 18 int ans=0,f=1; 19 char c=getchar(); 20 while(!isdigit(c)){ 21 if(c=='-') f=-1; 22 c=getchar(); 23 } 24 while(isdigit(c)){ 25 ans=ans*10+c-'0'; 26 c=getchar(); 27 } 28 return ans*f; 29 } 30 const int maxn=100009; 31 int n; 32 ll w[maxn]; 33 double a; 34 void solve(int k,int l,int r){ 35 double ans=0; 36 if(r-l<=500){ 37 rep(i,l,r+1) ans+=1.0*(w[i]-w[i-1])/(k-i); 38 }else{ 39 while(l<=r){ 40 int t=min(l+500,r); 41 ans+=1.0*(w[t]-w[l-1])/(k-(l+t)/2); 42 l=t+1; 43 } 44 } 45 printf("%.7lf\n",ans*(w[k]-w[k-1])); 46 } 47 int main(){ 48 n=read();scanf("%lf",&a); 49 rep(i,1,n+1) w[i]=w[i-1]+read(); 50 rep(i,1,n+1) solve(i,1,floor(a*i)); 51 return 0; 52 }
1011: [HNOI2008]遙遠的行星
Time Limit:?10 Sec??Memory Limit:?162 MBSec??Special JudgeSubmit:?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
3
5
6
2
4
Sample Output
0.000000
0.000000
0.000000
1.968750
2.976000
0.000000
0.000000
1.968750
2.976000
HINT
?
精確結果應該為0 0 0 2 3,但樣例輸出的結果誤差不超過5%,也算對
?