菜 WA0
查看原帖
菜 WA0
54372
A_Đark_Horcrux楼主2021/10/18 22:57

开了long long

公式没打错

#include<cstdio>
#include<cmath> 
#define il inline
using namespace std;
typedef long long ll;
typedef double db;
const int N=1e6+10;
inline ll read()
{
	ll s=0; char c=getchar();
	while(c<'0'||c>'9') c=getchar();
	while(c>='0'&&c<='9') s=(s<<3)+(s<<1)+(c&15),c=getchar();
	return s;
}
ll n,m,t,l,r,i,a[N],tag[N];
db k,s[N],c[N];
il ll ls(ll x) {return x<<1;}
il ll rs(ll x) {return (x<<1)|1;}
il void pushup(ll x)
{
	s[x]=s[ls(x)]+s[rs(x)];
	c[x]=c[ls(x)]+c[rs(x)];
	return ;
}
void update(int now,db sinx,db cosx)
{
	db siny=s[now],cosy=c[now];
    s[now]=sinx*cosy+cosx*siny;
    c[now]=cosx*cosy-sinx*siny;
    return ;
}
void pushdown(ll now)
{
	if(!tag[now]) return ;
	update(ls(now),sin(tag[now]),cos(tag[now]));
	update(rs(now),sin(tag[now]),cos(tag[now]));
	tag[ls(now)]+=tag[now];
	tag[rs(now)]+=tag[now];
	tag[now]=0; return ;
}
void build(ll now,ll l,ll r)
{
	tag[now]=0;
	if(l==r)
	{
		s[now]=sin(a[l]);
        c[now]=cos(a[l]);
		return ;
	}
	ll mid=(l+r)>>1;
	build(ls(now),l,mid);
	build(rs(now),mid+1,r);
	pushup(now); return ;
}
void add(ll x,ll y,ll l,ll r,ll now,ll k)
{
	if(x<=l&&r<=y)
	{
		update(now,s[k],c[k]);
		tag[now]+=k; return ;
	}
	pushdown(now);
	ll mid=(l+r)>>1;
	if(x<=mid) add(x,y,l,mid,ls(now),k);
	if(y>mid) add(x,y,mid+1,r,rs(now),k);
	pushup(now); return ;
}
db ask(ll x,ll y,ll l,ll r,ll now)
{
	if(x<=l&&r<=y) return s[now];
	db cnt=0; ll mid=(l+r)>>1;
	pushdown(now);
	if(x<=mid) cnt+=ask(x,y,l,mid,ls(now));
	if(y>mid) cnt+=ask(x,y,mid+1,r,rs(now));
	return cnt;
}
int main()
{
	n=read();
	for(i=1;i<=n;i++) a[i]=read();
	build(1,1,n);
	m=read();
	while(m--)
	{
		t=read();
		if(t==1)
			l=read(),r=read(),k=read(),add(l,r,1,n,1,k);
		else
			l=read(),r=read(),printf("%.1lf\n",ask(l,r,1,n,1));
	}
	return 0;
}
2021/10/18 22:57
加载中...