csp202309-2

题目

问题描述
对于平面直角坐标系上的坐标 (x,y)，小 P 定义了如下两种操作：

拉伸 k 倍：横坐标 x 变为 kx，纵坐标 y 变为 ky；

旋转 θ：将坐标 (x,y) 绕坐标原点 (0,0) 逆时针旋转 θ 弧度（0≤θ<2π）。易知旋转后的横坐标为 xcos⁡θ−ysin⁡θ，纵坐标为 xsin⁡θ+ycos⁡θ。

设定好了包含 n 个操作的序列 (t1,t2,⋯,tn) 后，小 P 又定义了如下查询：

i j x y：坐标 (x,y) 经过操作 ti,⋯,tj（1≤i≤j≤n）后的新坐标。
对于给定的操作序列，试计算 m 个查询的结果。

输入格式
从标准输入读入数据。

输入共 n+m+1 行。

输入的第一行包含空格分隔的两个正整数 n 和 m，分别表示操作和查询个数。

接下来 n 行依次输入 n 个操作，每行包含空格分隔的一个整数（操作类型）和一个实数（k 或 θ），形如 1 k（表示拉伸 k 倍）或 2 θ（表示旋转 θ）。

接下来 m 行依次输入 m 个查询，每行包含空格分隔的四个整数 i、j、x 和 y，含义如前文所述

输出格式
输出到标准输出中。

输出共
行，每行包含空格分隔的两个实数，表示对应查询的结果。

样例输入

10 5
2 0.59
2 4.956
1 0.997
1 1.364
1 1.242
1 0.82
2 2.824
1 0.716
2 0.178
2 4.094
1 6 -953188 -946637
1 9 969538 848081
4 7 -114758 522223
1 9 -535079 601597
8 8 159430 -511187

样例输出

-1858706.758 -83259.993
-1261428.46 201113.678
-75099.123 -738950.159
-119179.897 -789457.532
114151.88 -366009.892

80% 的测试数据满足：n,m<=1000.

ac代码

当直接用数组存储每个步骤的操作数时会出现只能过80%的用例情况
这时候要考虑换方法，类前缀和？

#include<iostream>
#include<math.h>
using namespace std;
double a[1000001][2];

void move(int i ,int j, double x,double y)
{
	double k ,t;
	k=a[j][0]/a[i-1][0];
	t=a[j][1]-a[i-1][1];
	double tmp=x;
	x=x*cos(t)-y*sin(t);
	y=tmp*sin(t)+y*cos(t);
	printf("%lf %lf\n",k*x,k*y);
}


int main()
{
	int n,m;
	cin>>n>>m;
	
	int start,end;
	double x,y;
	int type;
	double act;
	a[0][0]=1;
	a[0][1]=0;
	for(int i =1;i<=n;i++)
	{
		scanf("%d %lf",&type,&act);
		if(type==1)
		{
			a[i][0]=a[i-1][0]*act;
			a[i][1]=a[i-1][1];
		}
		else
		{
			a[i][1]=a[i-1][1]+act;	
			a[i][0]=a[i-1][0];
		}
	}
	for(int j =0;j<m;j++)
	{
		scanf("%d %d %lf %lf",&start,&end,&x,&y);
		move(start,end,x,y);
	}
	return 0;
}