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.
1< a < =3.5 接下来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%,也算对
题解
这题考的应该是对数学式子的转化,坑爹地卡精度。不看题解真没想出来。这里讲得很详细:
http://hi.baidu.com/zeonsgtr/item/789da6f2838a3dc742c36ab7
写了好几版代码。
第一版:4884ms。
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<cmath>
#include<algorithm>
#define err 1e-8
#define t 1000
using namespace std;
int n;
double A,m[100005],f[100005];
void init()
{
scanf("%d%lf",&n,&A);
int i;
for(i=1;i<=n;i++) scanf("%lf",&m[i]);
}
void work()
{
int i,j,l,r;
double w;
for(j=1;j<=t;j++)
{r=(int) (err+floor(A*j));
for(i=1;i<=r;i++)
{f[j]+=(double)m[i]*m[j]/double(j-i);}
}
/*for(i=1;i<=t;i++)
printf("%.6lf\n",f[i]);*/
//--------------------------------------------------以上是暴力
for(j=t+1;j<=n;j++)
{l=(int) (err+floor(A*(j-t)));
r=(int) (err+floor(A*j));
for(i=l+1;i<=r;i++)
f[j]+=m[j]*m[i]/(j-i);
f[j]+=(double)m[j]*(f[j-t]/m[j-t])*(j-t-(double)l/2)/double(j-(double)l/2);
}
for(int i=1;i<=n;i++) printf("%.6lf\n",f[i]);
}
int main()
{
init(); work();
return 0;
}
第二版:把暴力的部分t改到100:1936ms。
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<cmath>
#include<algorithm>
#define err 1e-8
#define t 100
using namespace std;
int n;
double A,m[100005],f[100005];
void init()
{
scanf("%d%lf",&n,&A);
int i;
for(i=1;i<=n;i++) scanf("%lf",&m[i]);
}
void work()
{
int i,j,l,r;
double w;
for(j=1;j<=t;j++)
{r=(int) (err+floor(A*j));
for(i=1;i<=r;i++)
{f[j]+=(double)m[i]*m[j]/double(j-i);}
}
for(j=t+1;j<=n;j++)
{l=(int) (err+floor(A*(j-t)));
r=(int) (err+floor(A*j));
for(i=l+1;i<=r;i++)
f[j]+=m[j]*m[i]/(j-i);
f[j]+=(double)m[j]*(f[j-t]/m[j-t])*(j-t-(double)l/2)/double(j-(double)l/2);
}
for(int i=1;i<=n;i++) printf("%.6lf\n",f[i]);
}
int main()
{
init(); work();
return 0;
}
第三版:t=15 1664ms
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<cmath>
#include<algorithm>
#define err 1e-8
#define t 15
using namespace std;
int n;
double A,m[100005],f[100005];
void init()
{
scanf("%d%lf",&n,&A);
int i;
for(i=1;i<=n;i++) scanf("%lf",&m[i]);
}
void work()
{
int i,j,l,r;
double w;
for(j=1;j<=t;j++)
{r=(int) (err+floor(A*j));
for(i=1;i<=r;i++)
{f[j]+=(double)m[i]*m[j]/double(j-i);}
}
for(j=t+1;j<=n;j++)
{l=(int) (err+floor(A*(j-t)));
r=(int) (err+floor(A*j));
for(i=l+1;i<=r;i++)
f[j]+=m[j]*m[i]/(j-i);
f[j]+=(double)m[j]*(f[j-t]/m[j-t])*(j-t-(double)l/2)/double(j-(double)l/2);
}
for(int i=1;i<=n;i++) printf("%.6lf\n",f[i]);
}
int main()
{
init(); work();
return 0;
}