HDU 4087 三维上的平移缩放旋转矩阵变化
题目大意:
就是根据它给的程序的要求,不断平移,缩放,旋转三维的点,最后计算出点的位置
这里主要是要列出三种转换方式的齐次矩阵描述
平移
translate tx ty tz
1 0 0 0
0 1 0 0
0 0 1 0
tx ty tz 1
缩放
scale a b c
a 0 0 0
0 b 0 0
0 0 c 0
0 0 0 1
绕任意轴(过原点)旋转(注意要把轴向量归一化,否则点在旋转轴上时有问题)
这里是以(x,y,z)向量指向我们人的方向逆时针旋转 d 的弧度
rotate x y z d
(1-cos(d))*x*x+cos(d) (1-cos(d))*x*y-sin(d)*z (1-cos(d))*x*z+sin(d)*y 0
(1-cos(d))*y*x+sin(d)*z (1-cos(d))*y*y+cos(d) (1-cos(d))*y*z-sin(d)*x 0
(1-cos(d))*z*x-sin(d)*y (1-cos(d))*z*y+sin(d)*x (1-cos(d))*z*z+cos(d) 0
0 0 0 1
然后这里因为循环的问题,所以用矩阵快速幂加速
而循环是可能出现嵌套的
我没用递归,而是判断当前循环属于第几个矩阵下的计算,每次进入一个新的循环都会给当前新的循环设置一个编号,以便找到它的矩阵
每次退出循环,都会将当前矩阵和循环外那个矩阵相乘即可
最后输出+eps,防止 -0.0 , -0.0 , -0.0 的情况发生
#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
#include <iostream>
using namespace std;
#define N 1005
#define eps 1e-6
const double PI = acos(-1.0);
char s[];
int rec[N]; struct Matrix{
double m[][];
Matrix operator*(const Matrix &t) const {
Matrix ans;
for(int i= ; i< ; i++){
for(int j= ; j< ; j++){
ans.m[i][j] = ;
for(int k= ; k< ; k++){
ans.m[i][j]+=m[i][k]*t.m[k][j];
}
}
}
return ans;
}
Matrix(){}
Matrix(double p[][])
{
for(int i= ; i< ; i++)
for(int j= ; j< ; j++)
m[i][j] = p[i][j];
}
void init(){
memset( m , , sizeof(m));
for(int i= ; i< ; i++) m[i][i] = ;
}
void out(){
for(int i= ; i< ; i++){
for(int j= ; j< ; j++)
cout<<m[i][j]<<" ";
cout<<endl;
}
}
}mat[] , tmp; Matrix q_pow(Matrix a , int k)
{
Matrix ret;
ret.init();
while(k)
{
if(k&) ret = ret*a;
a = a*a;
k>>=;
}
return ret;
}
//当向量方向指向自己的时候逆时针旋转A的弧度
void Rotate(Matrix &p , double a , double b , double c , double A)
{
double len = sqrt(a*a+b*b+c*c);
double x = a/len , y = b/len , z = c/len;
double sine = sin(A) , cosine = cos(A);
double m[][] = {
{cosine+(-cosine)*x*x, x*y*(-cosine)-z*sine, x*z*(-cosine)+y*sine, },
{y*x*(-cosine)+z*sine, cosine+y*y*(-cosine), y*z*(-cosine)-x*sine, },
{z*x*(-cosine)-y*sine, z*y*(-cosine)+x*sine, cosine+z*z*(-cosine), },
{ , , , }
};
p = Matrix(m);
} void Trans(Matrix &p , double a , double b , double c)
{
p.init();
p.m[][]=a , p.m[][]=b , p.m[][]=c;
} void Scale(Matrix &p , double a , double b , double c)
{
p.init();
p.m[][] = a , p.m[][] = b , p.m[][] = c;
} int main()
{
// freopen("in.txt" , "r" , stdin);
int n;
double a,b,c,d; while(scanf("%d" , &n) , n){
mat[].init();
int cnt = ;//repeat次数
while(scanf("%s" , s)){
if(s[]=='e' && cnt==) break;
if(s[]=='e'){
// cout<<"id: "<<cnt<<": "<<endl;
// mat[cnt].out();
mat[cnt] = q_pow(mat[cnt] , rec[cnt]);
// cout<<"id: "<<cnt<<": "<<endl;
// mat[cnt].out();
mat[cnt-] = mat[cnt-]*mat[cnt];
// cout<<"id: "<<cnt<<": "<<endl;
// mat[cnt-1].out();
cnt--;
}
if(s[]=='r' && s[]=='e') {
scanf("%d" , &rec[++cnt]);
mat[cnt].init();
}
else if(s[]=='t'){
scanf("%lf%lf%lf" , &a , &b , &c);
Trans(tmp , a , b , c);
// tmp.out();
mat[cnt] = mat[cnt]*tmp;
}
else if(s[]=='s'){
scanf("%lf%lf%lf" , &a , &b , &c);
Scale(tmp , a , b , c);
mat[cnt] = mat[cnt]*tmp;
}
else if(s[]=='r' && s[]=='o'){
scanf("%lf%lf%lf%lf" , &a , &b , &c , &d);
Rotate(tmp , a , b , c , -d/*PI);
mat[cnt] = mat[cnt]*tmp;
}
}
// mat[0].out();
for(int i= ; i<n ; i++){
double x , y , z , xx , yy , zz;
scanf("%lf%lf%lf" , &x , &y , &z);
xx = x*mat[].m[][]+y*mat[].m[][]+z*mat[].m[][]+mat[].m[][];
yy = x*mat[].m[][]+y*mat[].m[][]+z*mat[].m[][]+mat[].m[][];
zz = x*mat[].m[][]+y*mat[].m[][]+z*mat[].m[][]+mat[].m[][];
printf("%.2f %.2f %.2f\n" , xx+eps , yy+eps , zz+eps);
}
puts("");
}
return ;
}
HDU 4087 三维上的平移缩放旋转矩阵变化的更多相关文章
- 【OpenGL(SharpGL)】支持任意相机可平移缩放的轨迹球实现
[OpenGL(SharpGL)]支持任意相机可平移缩放的轨迹球 (本文PDF版在这里.) 在3D程序中,轨迹球(ArcBall)可以让你只用鼠标来控制模型(旋转),便于观察.在这里(http://w ...
- WPF中的平移缩放和矩阵变换(TranslateTransform、ScaleTransform、MatrixTransform)
在WPF中的平移缩放都是通过RenderTransform这个类来实现这些效果的,在这个类中,除了平移和缩放还有旋转.扭曲变换.矩阵变换.这些都差不多的,都是坐标的变换. 这里我就先简单弄个平移和缩放 ...
- 一步一步手写GIS开源项目-(2)地图平移缩放实现
系列文章目录 一步一步手写GIS开源项目-(1)500行代码实现基础GIS展示功能 一步一步手写GIS开源项目-(2)地图平移缩放实现 项目github地址:https://github.com/Hu ...
- WPF/Silverlight中图形的平移,缩放,旋转,倾斜变换演示
原文:WPF/Silverlight中图形的平移,缩放,旋转,倾斜变换演示 为方便描述, 这里仅以正方形来做演示, 其他图形从略. 运行时效果图:XAML代码:// Transform.XAML< ...
- CSS 实现背景图尺寸不随浏览器缩放而变化
<!-- Author:博客园小dee --> 一些网站的首页背景图尺寸不随浏览器缩放而变化,例如百度个人版的首页,缩放后背景图的尺寸并不改变: 再比如花瓣网( http://www.hu ...
- CSS实现背景图尺寸不随浏览器缩放而变化
方法一. 把图片作为background,方法二使用img标签.同时要有一张足够大尺寸的图片. 方法一. 把图片作为background 有几个CSS的属性要提一下:background-size:c ...
- HDU 1253 三维数组的图上找最短路
题目大意: 从三维空间的(0,0,0)出发到(a-1,b-1,c-1),每移动一个都要时间加一,计算最短时间 根据六个方向,开个bfs,像spfa那样计算最短路径就行了,但是要1200多ms,也不知道 ...
- Graphics平移缩放旋转(转载)+点睛
点睛:可以进行多次旋转和平移,也就是平移以后再平移,旋转以后再旋转,有时候一次达不到要求,如果你想一次调整完美的话很麻烦,所以最好多次,上代码 private void btnTranslate_Cl ...
- C#图片上传服务器缩放存储功能
项目需求上需要用户上传头像.本来是用第三方的插件的.但是很多都是收费的. 不过,我需要花这钱么?是不?所以网络上搜集了些资料.. 果然.这个世界 大神是很多很多的. 用了个免费插件. <scri ...
随机推荐
- python 脚本传递参数
python查找指定字符 #!/usr/bin/env python import sys import re f = open("log.txt", "rb" ...
- 转!!java泛型概念(泛型类,接口,方法)
一. 泛型概念的提出(为什么需要泛型)? 首先,我们看下下面这段简短的代码: 1 public class GenericTest { 2 3 public static void main(Stri ...
- 图解javascript中this指向
JavaScript 是一种脚本语言,支持函数式编程.闭包.基于原型的继承等高级功能.JavaScript一开始看起来感觉会很容易入门,但是随着使用的深入,你会发JavaScript其实很难掌握,有些 ...
- 抛弃vboot不格盘用grub4dos+firadisk安装Ghost版XP到VHD,轻松RAMOS!
http://bbs.wuyou.net/forum.php?mod=viewthread&tid=363198&extra=抛弃vboot不格盘用grub4dos+firadisk安 ...
- [saiku] 访问saiku首页的时候前后台处理流程
这篇文章讲述:项目启动后,首次访问SAIKU的登录页,前后台分别做了什么处理 (1) 访问的到底是什么页面? 浏览器输入:localhost:8080 啪一回车 根据web访问的尿性,访问的是 ind ...
- Qt之QSS(黑色炫酷)
简述 Qt助手中有关于各种部件的QSS详细讲解,资源很丰富,请参考:Qt Style Sheets Examples. 黑色炫酷 - 一款漂亮的QSS风格. 之前博客中分享了很多关于Qt的样式效果,几 ...
- (07)odoo扩展API
* 打开XML-RPC 连接 >>> import xmlrpclib >>> srv, db = 'http://localhost:8069', ' ...
- 对Web标准的理解
Web标准是一系列标准的组合,包括结构化语言标准(XHTML/XML/HTML).表现标准语言(CSS)及行为标准语言(JavaScript).具体表现为: (1)Web标准规范要求,标签书写必须闭合 ...
- Java 之 I/O 系列 02 ——序列化(一)
Java 之 I/O 系列 目录 Java 之 I/O 系列 01 ——基础 Java 之 I/O 系列 02 ——序列化(一) Java 之 I/O 系列 02 ——序列化(二) 一 序列化概述 序 ...
- ERP开发分享 1 数据库表设计
这是我的ERP设计经验分享系列,今天讲的是数据库的表设计(1),主要阐述: 1.单字段的主键:2.使用int32作为主键类型:3.使用版本字段处理乐观锁定:4.生效字段标明是否允许“被使用”:5.锁定 ...