Moorhsum Round #1 获取名额 乱搞

38 篇文章 0 订阅
9 篇文章 0 订阅

Description


给定长为n的序列a[],m次询问(l,r,x)求这样一个柿子 1 − ∏ i = l r ( 1 − a i x ) 1-\prod\limits_{i=l}^{r}{\left(1-a_ix\right)} 1i=lr(1aix)
n ≤ 6 × 1 0 5 a i , x ≤ 1 0 9 n\le6\times10^5\\a_i,x\le10^9 n6×105ai,x109

Solution


正解太强了%%%%
首先四十分是送的。。
然后有七十分是我们把x当作未知数做出这个多项式,然后等价于多点求值问题
一个并不能过的响法就是将询问区间内的a排序,这样就收敛得比较快可以满足精度了

考虑优化上面那个不成熟的响法,我们对柿子两边ln变成加法最后exp回去。最终答案可以分成两部分,分别是比较大直接算的,和其余的用ln的和做

考虑把 ln ⁡ ( 1 − a i x ) \ln\left(1-\frac{a_i}{x}\right) ln(1xai)泰勒展开取前M项系数,那么我们只需要带进当前的x就能算出指数了
然后我们根据精度需要每次找区间内最大的数字直接乘,剩余的就用exp算

Code


#include <stdio.h>
#include <string.h>
#include <algorithm>
#include <math.h>
#define rep(i,st,ed) for (int i=st;i<=ed;++i)

const double eps=1e-8;
const int N=1000005;

int max[21][N],lg[N];

double s[21][N],a[N],rec[21],x,ans;

int read() {
	int x=0,v=1; char ch=getchar();
	for (;ch<'0'||ch>'9';v=(ch=='-')?(-1):(v),ch=getchar());
	for (;ch<='9'&&ch>='0';x=x*10+ch-'0',ch=getchar());
	return x*v;
}

int query(int l,int r) {
	if (l==r) return l;
	int llgg=lg[r-l+1];
	l=max[llgg][l],r=max[llgg][r-(1<<llgg)+1];
	return (a[l]>a[r]?l:r);
}

void solve(int l,int r) {
	if (r<l||ans<eps) return ;
	int pos=query(l,r);
	if (a[pos]*x>0.5) {
		rep(i,1,20) rec[i]-=s[i][pos]-s[i][pos-1];
		ans=ans*(1-a[pos]*x);
		solve(l,pos-1);
		solve(pos+1,r);
	}
}

int main(void) {
	freopen("data.in","r",stdin);
	freopen("myp.out","w",stdout);
	int n=read(),m=read(); double mx=0;
	rep(i,1,n) a[i]=read(),mx=std:: max(mx,a[i]);
	rep(i,1,n) {
		if (i>1) lg[i]=lg[i>>1]+1;
		max[0][i]=i; a[i]/=mx;
		double tmp=a[i];
		rep(j,1,20) s[j][i]=s[j][i-1]-tmp/j,tmp*=a[i];
	}
	for (int i=1;i<=n;++i) printf("%d\n", lg[i]); puts("");
	rep(j,1,lg[n]) rep(i,1,n-(1<<j)+1) {
		int l=max[j-1][i],r=max[j-1][i+(1<<j-1)];
		if (a[l]>a[r]) max[j][i]=l;
		else max[j][i]=r;
	}
	for (;m--;) {
		int l=read(),r=read(); x=read();
		ans=1; x=mx/x;
		rep(i,1,20) rec[i]=s[i][r]-s[i][l-1];
		solve(l,r);
		double tot=0,tmp=1;
		rep(i,1,20) tmp=tmp*x,tot+=rec[i]*tmp;
		ans=ans*exp(tot);
		printf("%.10lf\n", 1-ans);
	}
	return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值