网站联盟三要素,wordpress开发用什么软件,沧州网站运营,网站如何建设手机版Description 直线上N颗行星#xff0c;Xi处有行星i,行星J受到行星I的作用力#xff0c;当且仅当iAJ.此时J受到作用力的大小为 Fi-jMi*Mj/(j-i) 其中A为很小的常量#xff0c;故直观上说每颗行星都只受到距离遥远的行星的作用。请计算每颗行星的受力 #xff0c;只要…Description 直线上N颗行星Xi处有行星i,行星J受到行星I的作用力当且仅当iAJ.此时J受到作用力的大小为 Fi-jMi*Mj/(j-i) 其中A为很小的常量故直观上说每颗行星都只受到距离遥远的行星的作用。请计算每颗行星的受力 只要结果的相对误差不超过5%即可. Input 第一行两个整数N和A. 1N10^5.0.01 a 0.35,接下来N行输入N个行星的质量Mi,保证0Mi10^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 Solution 神仙操作。。 注意到只要结果的相对误差不超过5%即可可以考虑近似计算。 对于第\(i\)个位置答案显然就是\[ \sum_{k1}^{\lfloor a\cdot i\rfloor}\frac{m_im_k}{i-k} \] 如果\(a\cdot i\)较小直接计算就好了。 否则可以考虑分成若干段近似计算对于一个区间\([a,b]\)可以把答案近似为\(m_i\cdot \frac{\sum_{ka}^{b}m_k}{i-(ab)/2}\)其实就是把分母全看成一样的。 显然分的段数越高效率越低答案越精确我这里是分了\(100\)段。 #includebits/stdc.h
using namespace std;void read(int x) {x0;int f1;char chgetchar();for(;!isdigit(ch);chgetchar()) if(ch-) f-f;for(;isdigit(ch);chgetchar()) xx*10ch-0;x*f;
}#define write(x) printf(%d\n,x)
#define ll long long
#define lf doubleconst int maxn 2e510;
const int T 100;int m[maxn],n;
ll sum[maxn];
lf a;int main() {read(n);scanf(%lf,a);for(int i1;in;i) read(m[i]),sum[i]sum[i-1]1ll*m[i];for(int i1;in;i) {int xfloor(a*(lf)i);double ans0; if(x1000) for(int k1;kx;k) ans(lf)m[i]*m[k]/(lf)(i-k);else {int tx/T;for(int k1;kx;kt) {int edmin(x,kt-1);int mid(ked)1;ans(lf)m[i]*(sum[ed]-sum[k-1])/(lf)(i-mid);}}printf(%.6lf\n,ans);}return 0;
} 转载于:https://www.cnblogs.com/hbyer/p/10351771.html
