poj 3083

  1 #include <iostream>
  2 #include <limits>
  3 #include <queue>
  4 using namespace std;
  5 
  6 typedef pair<int, int> P;
  7 
  8 char maze[41][41]; 
  9 int dist[41][41];
 10 int spDirX[4] = {1, 0, -1, 0};
 11 int spDirY[4] = {0, 1, 0, -1};
 12 int n;
 13 int w, h;
 14 int startX, startY;
 15 int endX, endY;
 16 const int INF = numeric_limits<int>::max();
 17 int sp, lo, ro;
 18 int dx[4] = {
 19         0, 1, 0, -1
 20 };
 21 int dy[4] = {
 22         -1, 0, 1, 0    
 23 };
 24 void leftOrder(int x, int y, int sum, int dir)
 25 {
 26     //ÖÕÖ¹ 
 27     if(x == endX && y == endY)
 28     {
 29         lo = sum;
 30         return ;
 31     }
 32     //µÚÒ»¸öÔªËØÈ·¶¨·½Ïò
 33     if(sum == 1)
 34     {
 35         for(int i = 0; i < 4; ++i)
 36         {
 37             int nx = x + dx[i];
 38             int ny = y + dy[i];
 39             if(nx >= 0 && nx < w && ny >= 0 && ny < h && maze[ny][nx] == '.')
 40             {
 41                     leftOrder(nx, ny , sum+1, i);
 42                     break;
 43             }
 44             
 45         }
 46     }
 47     else
 48     {
 49         //Ò»°ãÇé¿ö 
 50         dir = (dir+3) % 4;
 51         for(int i = 0; i < 4; ++i)
 52         {
 53             int nx = x + dx[(dir+i)%4];
 54             int ny = y + dy[(dir+i)%4];
 55             if(nx >= 0 && nx < w && ny >= 0 && ny < h && maze[ny][nx] != '#')
 56             {
 57                 leftOrder(nx, ny , sum+1, (dir+i)%4);
 58                 break;
 59             }        
 60         }    
 61     }
 62     
 63 }
 64 
 65 void rightOrder(int x, int y, int sum, int dir)
 66 {
 67         //ÖÕÖ¹ 
 68     if(x == endX && y == endY)
 69     {
 70         ro = sum;
 71         return ;
 72     }
 73     //µÚÒ»¸öÔªËØÈ·¶¨·½Ïò
 74     if(sum == 1)
 75     {
 76         for(int i = 0; i < 4; ++i)
 77         {
 78             int nx = x + dx[i];
 79             int ny = y + dy[i];
 80             if(nx >= 0 && nx < w && ny >= 0 && ny < h && maze[ny][nx] == '.')
 81             {
 82                     rightOrder(nx, ny , sum+1, i);
 83                     break;
 84             }
 85             
 86         }
 87     }
 88     else
 89     {
 90         //Ò»°ãÇé¿ö 
 91         dir = (dir+1) % 4;
 92         for(int i = 0; i < 4; ++i)
 93         {
 94             int nx = x + dx[(dir-i+4)%4];
 95             int ny = y + dy[(dir-i+4)%4];
 96             if(nx >= 0 && nx < w && ny >= 0 && ny < h && maze[ny][nx] != '#')
 97             {
 98                 rightOrder(nx, ny , sum+1, (dir-i+4)%4);
 99                 break;
100             }        
101         }    
102     }
103     
104 } 
105 
106 int shortPath()
107 {
108     queue<P> que; 
109     dist[startY][startX] = 1;
110     que.push(P(startX, startY));
111     while(que.size())
112     {
113         P temp = que.front();
114         que.pop();
115         if(temp.first == endX && temp.second == endY)
116             break;
117         for(int i = 0; i < 4; ++i)
118         {
119             int nx = temp.first + spDirX[i];
120             int ny = temp.second + spDirY[i];
121             if(nx >= 0 && ny >= 0 && nx < w && ny < h && dist[ny][nx] == INF && maze[ny][nx] != '#')
122             {
123                 que.push(P(nx, ny));
124                 dist[ny][nx] = dist[temp.second][temp.first] + 1;
125             }
126         }        
127     }
128     sp = dist[endY][endX];
129 }
130 
131 
132 int main()
133 {
134     cin >> n;
135     for(int i = 0; i < n; ++i)
136     {
137         cin >> w >> h;
138         for(int j = 0; j < h; ++j)
139         {
140             for(int k = 0; k < w; ++k)
141             {
142                 cin >> maze[j][k];
143                 dist[j][k] = INF;
144                 if(maze[j][k] == 'S')
145                 {
146                     startX = k;
147                     startY = j;
148                 }
149                 if(maze[j][k] == 'E')
150                 {
151                     endX = k;
152                     endY = j;
153                 }    
154             }
155         }
156         leftOrder(startX, startY, 1, 0);
157         rightOrder(startX, startY, 1, 0);
158         shortPath();    
159         cout << lo << " " << ro << " " << sp << endl;    
160     }
161     return 0;
162 } 
View Code

 

转载于:https://www.cnblogs.com/shuanghong/p/4075237.html

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
### 支持向量机非线性回归通用MATLAB程序解析 #### 一、概述 本文将详细介绍一个基于MATLAB的支持向量机(SVM)非线性回归的通用程序。该程序采用支持向量机方法来实现数据的非线性回归,并通过不同的核函数设置来适应不同类型的数据分布。此外,该程序还提供了数据预处理的方法,使得用户能够更加方便地应用此程序解决实际问题。 #### 二、核心功能与原理 ##### 1. 支持向量机(SVM) 支持向量机是一种监督学习模型,主要用于分类和回归分析。对于非线性回归任务,SVM通过引入核技巧(kernel trick)将原始低维空间中的非线性问题转换为高维空间中的线性问题,从而实现有效的非线性建模。 ##### 2. 核函数 核函数的选择直接影响到模型的性能。本程序内置了三种常用的核函数: - **线性核函数**:`K(x, y) = x'y` - **多项式核函数**:`K(x, y) = (x'y + 1)^d` - **径向基函数(RBF)**:`K(x, y) = exp(-γ|x - y|^2)` 其中RBF核函数被广泛应用于非线性问题中,因为它可以处理非常复杂的非线性关系。本程序默认使用的是RBF核函数,参数`D`用于控制高斯核函数的宽度。 ##### 3. 数据预处理 虽然程序本身没有直接涉及数据预处理的过程,但在实际应用中,对数据进行适当的预处理是非常重要的。常见的预处理步骤包括归一化、缺失值处理等。 ##### 4. 模型参数 - **Epsilon**: ε-insensitive loss function的ε值,控制回归带宽。 - **C**: 松弛变量的惩罚系数,控制模型复杂度与过拟合的风险之间的平衡。 #### 三、程序实现细节 ##### 1. 函数输入与输出 - **输入**: - `X`: 输入特征矩阵,维度为(n, l),其中n是特征数量,l是样本数量。 - `Y`: 目标值向量,长度为l。 - `Epsilon`: 回归带宽。 - `C`: 松弛变量的惩罚系数。 - `D`: RBF核函数的参数。 - **输出**: - `Alpha1`: 正的拉格朗日乘子向量。 - `Alpha2`: 负的拉格朗日乘子向量。 - `Alpha`: 拉格朗日乘子向量。 - `Flag`: 标记向量,表示每个样本的类型。 - `B`: 偏置项。 ##### 2. 核心代码解析 程序首先计算所有样本间的核矩阵`K`,然后构建二次规划问题并求解得到拉格朗日乘子向量。根据拉格朗日乘子的值确定支持向量,并计算偏置项`B`。 - **核矩阵计算**:采用RBF核函数,通过`exp(-(sum((xi-xj).^2)/D))`计算任意两个样本之间的相似度。 - **二次规划**:构建目标函数和约束条件,使用`quadprog`函数求解最小化问题。 - **支持向量识别**:根据拉格朗日乘子的大小判断每个样本是否为支持向量,并据此计算偏置项`B`。 #### 四、程序扩展与优化 - **多核函数支持**:可以通过增加更多的核函数选项,提高程序的灵活性。 - **自动调参**:实现参数自动选择的功能,例如通过交叉验证选择最优的`Epsilon`和`C`值。 - **并行计算**:利用MATLAB的并行计算工具箱加速计算过程,特别是当样本量很大时。 #### 五、应用场景 该程序适用于需要进行非线性回归预测的场景,如经济预测、天气预报等领域。通过调整核函数和参数,可以有效应对各种类型的非线性问题。 ### 总结 本程序提供了一个支持向量机非线性回归的完整实现框架,通过灵活的核函数设置和参数调整,能够有效地处理非线性问题。对于需要进行回归预测的应用场景,这是一个非常实用且强大的工具。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值