题目链接：
https://jzoj.net/senior/#main/show/6087
题目：
题解：
- 只需要统计$\prod_{i=l}^r (1-\frac{a_i}{x})$
- =$exp(\sum_{i=l}^r ln(1-\frac{a_i}{x}))(x>a_i)$
- 我们可以把$ln(1-x)|x<1|$泰勒展开，得到$-\sum_{i=1}^{∞}\frac{x^i}{i}=0-\frac{x}{1}-\frac{x^2}{2}-\frac{x^3}{3}-...$
- 那么里面化简得到$-\sum_{i=1}^{r} \sum_{j=1}^∞\frac{a_i^j}{x^j*j})=-\sum_{i=1}^∞\frac{1}{i*x^i}\sum_{j=l}^{r}a_j^i$
- 实测$∞$取$50$就可以，预处理前缀和来快速计算$\sum_{j=l}^{r}a_j^i$
- 注意把$x=x/max(a_i)$,$a_i=a_i/max(a_i)$，防止爆$double$的范围
- 但是发现精度会有很大问题
- 当$\frac{a_i}{x}$较大时，$ln$里面的东西会比较接近$0$，很难保证精度
- 为$\frac{a_i}{x}$设置一个$lim$，当$\frac{a_i}{x}>lim$时直接算，再递归左右区间
- $lim$取$0.5$是可以过的
- 寻找区间内最大的$\frac{a_i}{x}$用$RMQ$
- 参考自大米饼大佬的博客
代码：
#include<bits/stdc++.h>
using namespace std;
typedef double db; const int N=6e5+;
const int M=;
int n,q;
int bin[],lg[N],f[N][];
db res;
db a[N],s[M+][N];
inline char gc()
{
static char*p1,*p2,s[];
if(p1==p2)p2=(p1=s)+fread(s,,,stdin);
return(p1==p2)?EOF:*p1++;
}
inline int read()
{
char ch=gc();int s=;
while (ch<''||ch>'') ch=gc();
while (ch>=''&&ch<='') {s=(s<<)+(s<<)+ch-'';ch=gc();}
return s;
}
int Max(int x,int y) {return a[x]>a[y]?x:y;}
void pre()
{
for (int i=bin[]=;i<=;bin[i]=bin[i-]<<,i++);
lg[]=-;
for (int i=;i<=n;i++) f[i][]=i,lg[i]=lg[i>>]+;
for (int i=;i<=;i++)
for (int j=;j+bin[i]-<=n;j++)
f[j][i]=Max(f[j][i-],f[j+bin[i-]][i-]);
}
int ask(int l,int r)
{
int t=lg[r-l+];
return Max(f[l][t],f[r-bin[t]+][t]);
}
db x;
db cal(int l,int r)
{
db re=,t=;
for (int i=;i<=M;i++)
{
t*=x;
re-=(s[i][r]-s[i][l-])/(i*t);
}
return exp(re);
}
void solve(int l,int r)
{
int mid=ask(l,r);
if (a[mid]/x<0.5) {res*=cal(l,r);return;}
res*=-a[mid]/x;
if (l<mid) solve(l,mid-);
if (r>mid) solve(mid+,r);
}
int main()
{
freopen("orz.in","r",stdin);
freopen("orz.out","w",stdout);
n=read();q=read();
db mx=;
for (int i=;i<=n;i++) a[i]=read(),mx=max(mx,a[i]);
for (int i=;i<=n;i++) a[i]/=mx;
pre();
for (int i=;i<=n;i++)
{
db t=;
for (int j=;j<=M;j++)
{
t*=a[i];
s[j][i]=s[j][i-]+t;
}
}
while (q--)
{
int l=read(),r=read();
x=read()/mx;
res=;solve(l,r);
printf("%.10f\n",-res);
}
return ;
}