csp 坐标变换(其二) C语言

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

拉伸 k 倍:横坐标 x 变为 kx,纵坐标 y 变为 ky;

旋转 a:将坐标 (x,y) 绕坐标原点 (0,0) 逆时针旋转 a弧度(0≤a<2*PI)。易知旋转后的横坐标为 xcos⁡(a)−ysin⁡(a),纵坐标为 xsin⁡(a)+ycos⁡(a)。

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

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

样例输入

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
解析:这个题真的真的就是错的很隐形,老是找不出来哪里有问题,最后居然是循环里的问题,笑死啊啊啊啊

错误代码大赏:
 

#include<stdio.h>
#include<math.h>

int main()
{
	int n,m;
	
	scanf("%d %d",&n,&m);
	int dx[n];
	double dy[n];
	for(int i=1;i<=n;i++)
	{
		scanf("%d %lf",&dx[i],&dy[i]);
	}
	int a[m],b[m],c[m],d[m];
	for(int j=1;j<=m;j++)
	{
		scanf("%d %d %d %d",&a[j],&b[j],&c[j],&d[j]);
	}
	for(int i=1;i<=m;i++)
	{
		double x=c[i], y=d[i];
		double zong,heng;
		for(int j=a[i];j<=b[i];j++)
		{
			if(dx[j]==1)
			{
				zong = x*dy[j];
				heng =y*dy[j];
			}
			else if(dx[j]==2)
			{
				double h=dy[j];
				double new_x = x*cos(h)-y*sin(h);
				double new_y = x*sin(h)+y*cos(h);
				zong = new_x;
				heng = new_y;
			}
		}
		printf("%.3f %.3f\n",zong,heng);
	}
	return 0;
}

错误解析:这个就是明显的定义了heng和zong导致那个循环只进行了一步,也就是说对于这个坐标的操作只进行了一次,它没办法在第一次操作后坐标的基础上进行第二次

正确代码:
 

#include<stdio.h>
#include<math.h>

int main()
{
	int n,m;
	scanf("%d %d",&n,&m);
	int dx[n];
	double dy[n];
	for(int i=1;i<=n;i++)
	{
		scanf("%d %lf",&dx[i],&dy[i]);
	}
	int a[m],b[m],c[m],d[m];
	for(int j=1;j<=m;j++)
	{
		scanf("%d %d %d %d",&a[j],&b[j],&c[j],&d[j]);
	}
	for(int i=1;i<=m;i++)
	{
		double x=c[i], y=d[i];
		for(int j=a[i];j<=b[i];j++)
		{
			if(dx[j]==1)
			{
				x*=dy[j];
				y*=dy[j];
			}
			else if(dx[j]==2)
			{
				double h=dy[j];
				double new_x = x*cos(h)-y*sin(h);
				double new_y = x*sin(h)+y*cos(h);
				x = new_x;
				y = new_y;
			}
		}
		printf("%.3f %.3f\n",x,y);
	}
	return 0;
}

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
根据引用[1]中的描述,本题中操作共有n个,坐标共有m个,且n和m均小于等于100000。因此,暴力算法的时间复杂度为O(n*m),会浪费宝贵的时间。我们需要寻找更高效的算法。 根据题目描述,我们需要进行坐标变换。具体来说,我们需要将一个矩阵中的所有点按照一定规则进行变换。这个规则可以通过观察样例来得到。 假设原始矩阵中的一个点的坐标为(x, y),则变换后的坐标为(x', y'),其中: x' = a1 * x + b1 * y + c1 y' = a2 * x + b2 * y + c2 其中a1, b1, c1, a2, b2, c2是给定的常数。 我们可以将这个变换看作是一个线性变换,即将原始矩阵中的每个点映射到一个新的矩阵中。因此,我们可以使用矩阵乘法来实现这个变换。 具体来说,我们可以将原始矩阵中的每个点表示为一个列向量(x, y, 1),然后将这个列向量与一个3x3的矩阵M相乘,得到一个新的列向量(x', y', w')。其中w'表示一个常数,我们可以将其忽略,因为我们只需要x'和y'。 因此,我们可以将矩阵M表示为: a1 b1 c1 a2 b2 c2 0 1 然后,我们可以将原始矩阵中的每个点表示为一个列向量(x, y, 1),并将这个列向量与矩阵M相乘,得到一个新的列向量(x', y', w')。最后,我们可以将x'和y'作为变换后的坐标。 下面是一个Python实现的例子: ```python import numpy as np # 读入数据 n, m = map(int, input().split()) a = np.zeros((n, 3)) for i in range(n): a[i] = list(map(int, input().split())) # 构造矩阵M a1, b1, c1, a2, b2, c2 = map(int, input().split()) M = np.array([[a1, b1, c1], [a2, b2, c2], [0, 0, 1]]) # 进行坐标变换 b = np.ones((m, 3)) for i in range(m): b[i, :2] = list(map(int, input().split())) c = np.dot(b, M.T) for i in range(m): print(int(c[i, 0]), int(c[i, 1])) ```

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值