第31次csp2023-9

第一题坐标变换（其一）

无脑做

#include<iostream>
using namespace std;
const int N=110,M=110;
int n,m;
int x[M],y[M];
int dx[N],dy[N]; 
void trans(int k){
    for(int i=0;i<n;i++){
        x[k]+=dx[i];
        y[k]+=dy[i];
    }
    return ;
}
int main(){
    //读入数据 
    scanf("%d%d",&n,&m);
    for(int i=0;i<n;i++) scanf("%d%d",&dx[i],&dy[i]);
    for(int i=0;i<m;i++) scanf("%d%d",&x[i],&y[i]);
    //执行操作
    for(int i=0;i<n;i++){
        trans(i);
    } 
    for(int i=0;i<m;i++) printf("%d %d\n",x[i],y[i]);
    
    return 0;
}

第二题坐标变换（其二）

暴力做法o(n^2) 会超时

暴力做法如下：

#include<iostream>
#include<cmath> 
using namespace std;
const int N=1e5+10;
int n,m;
double x,y;
int trans_type[N];
double trans_value[N];
void trans(int l,int r){
	for(int i=l;i<=r;i++){
		if(trans_type[i]==1){
			x=trans_value[i]*x;
			y=trans_value[i]*y;
		}
		else {
			double sita=trans_value[i];
			double temp1 =x*cos(sita)-y*sin(sita);
			double temp2 =x*sin(sita)+y*cos(sita);
			x=temp1;
			y=temp2;
		} 
	}
	return ; 
}
int main(){
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++) scanf("%d%lf",&trans_type[i],&trans_value[i]);
	while(m--){
		int l,r;
		scanf("%d%d%lf%lf",&l,&r,&x,&y);
		trans(l,r);
		printf("%.3f %.3f\n",x,y);
	} 
	return 0;
}

尝试结合前缀和思想进行优化，将每次操作看作对x乘上了某个系数，第l次操作到第r次操作就相当于Xr/Xl-1,但是跑样例不对，可能是坐标的伸长和旋转的混合运算在直角坐标系中不适用12212112

#include<iostream>
#include<cmath>
using namespace std;
const int N=1e6+10;
int n,m;
double changex[N],changey[N];
int main(){
	scanf("%d%d",&n,&m);
	changex[0]=1.0; changey[0]=1.0;
	for(int i=1;i<=n;i++){
		int type;double value;
		scanf("%d%lf",&type,&value);
        if(type==1) {
            changex[i]=changex[i-1]*value;
            changey[i]=changey[i-1]*value;
        }
        else{
            double temp1,temp2;
            temp1=changex[i-1]*cos(value)-changey[i-1]*sin(value);
            temp2=changex[i-1]*sin(value)+changey[i-1]*cos(value);
            changex[i]=temp1;
            changey[i]=temp2;
        }
	} 
	while(m--){ 
		int l,r;
		double x,y;
		scanf("%d%d%lf%lf",&l,&r,&x,&y);
		x*=changex[r]/changex[l-1];
		y*=changey[r]/changey[l-1];
		printf("%.3lf %.3lf\n",x,y);
	}
	
	return 0;
}

 改用极坐标系进行坐标变换即可，可以看到对坐标的伸长运算前缀和计算方法为len[r]/len[l-1],而对坐标的旋转运算计算方法为theta[r]-theta[l-1]。

#include<iostream>
#include<cmath>
using namespace std;
const int N=1e5+10;
int n,m;
double len[N],theta[N];
int main(){
	scanf("%d%d",&n,&m);
	len[0]=1.0;theta[0]=1.0;
	for(int i=1;i<=n;i++){
		int type;double value;
		scanf("%d%lf",&type,&value);
		if(type==1){
			len[i]=len[i-1]*value;
			theta[i]=theta[i-1];
		}
		else {
			len[i]=len[i-1];
			theta[i]=theta[i-1]+value;
		}
	}
	while(m--){
		int l,r;double x,y;
		scanf("%d%d%lf%lf",&l,&r,&x,&y);
		double tmplen=len[r]/len[l-1];
		double tmptheta=theta[r]-theta[l-1];
		x*=tmplen;
		y*=tmplen;
		double tmpx,tmpy; 	
		tmpx=x,tmpy=y;
		x=tmpx*cos(tmptheta)-tmpy*sin(tmptheta);
		y=tmpx*sin(tmptheta)+tmpy*cos(tmptheta);
		printf("%.3lf %.3lf\n",x,y);
	}
	return 0;
}