信息学奥赛一本通 1122:计算鞍点 | OpenJudge NOI 1.8 05:计算鞍点

【题目链接】

ybt 1122:计算鞍点
OpenJudge NOI 1.8 05:计算鞍点

【题目考点】

1. 二维数组
2. 求最大最小值

【思路及题解代码】

证明:证明:在一个n行n列的矩阵中,每行只有一个最大值,每列只有一个最小值。证明鞍点数量最多为1。

鞍点指的是矩阵中的一个元素,它是所在行的最大值,并且是所在列的最小值。
首先:鞍点数量可以为0
例:

1 2
2 1

鞍点数量可以为1:
例:

1

接下来使用反证法证明鞍点数量不可能大于等于2。
假设n行n列的矩阵中有大于等于2个鞍点
设其中两个鞍点分别为a[x1][y1], a[x2][y2]
由于a[x1][y1]是x1行的最大值,a[x2][y2]是x2行的最大值,且每行只有一个最大值,因此x1不等于x2。同理y1不等于y2。
由于a[x1][y1]是鞍点,那么一定存在a[x1][y2]和a[x2][y1]满足:
(1) a[x1][y2] < a[x1][y1]
(2) a[x2][y1] > a[x1][y1]
由于a[x2][y2]是鞍点,那么一定存在a[x1][y2]和a[x2][y1]满足:
(3) a[x2][y1] < a[x2][y2]
(4) a[x1][y2] > a[x2][y2]
根据(1)和(4)可得a[x2][y2] < a[x1][y2] < a[x1][y1]
根据(2)和(3)可得a[x1][y1] < a[x2][y1] < a[x2][y2]
a[x1][y1]和a[x2][y2]大大小关系出现矛盾,假设不成立,原命题得证。

解法1:

遍历各行,先找到这一行的最大值,假设最大值在第m_j列,然后判断该值是不是第m_j列的最小值,如果是,那么该位置就是鞍点。

#include<bits/stdc++.h>
using namespace std;
int main()
{
	bool isFound = false;//是否找到鞍点
    int a[6][6];
    for(int i = 1; i <= 5; ++i)
    	for(int j = 1; j <= 5; ++j)
    		cin >> a[i][j];
	for(int i = 1; i <= 5; ++i)//找到每行的最大值,而后判断该值是否是这一列的最小值
	{
		int m_j = 1;//m_j:最大值的纵坐标
		for(int j = 1; j <= 5; ++j)//遍历第i行,j是列号,找到最大值
			if(a[i][j] > a[i][m_j])
				m_j = j;//此时a[i][m_j]是第i行的最大值
		bool isAndian = true;//判断a[i][m_j]是否是第m_j列的最小值,如果是,那么a[i][m_j]就是鞍点
		for(int k = 1; k <= 5; ++k)//遍历第m_j列 k是行号
			if(a[k][m_j] < a[i][m_j])//只要在第m_j列找到一个小于a[i][m_j]的元素,那么该点就不是鞍点
			{
				isAndian = false;
				break;
			}
		if(isAndian)
		{
			cout << i << ' ' << m_j << ' ' << a[i][m_j] << endl;
			isFound = true;
			break;
		}
	}
	if(isFound == false)
		cout << "not found" << endl;
    return 0;
}

解法2:

设整型二维数组x,初始化为0。寻找各行的最大值的位置,若该位置是(i,j),那么将x[i][j]++。寻找各行最小值的位置,若位置是(i,j),将x[i][j]++。遍历数组x,如果发现某个位置x[i][j] == 2,那么该位置就是鞍点。否则就没有鞍点。

#include<bits/stdc++.h>
using namespace std;
int main()
{
    int a[6][6], x[6][6] = {};
    for(int i = 1; i <= 5; ++i)
    	for(int j = 1; j <= 5; ++j)
    		cin>>a[i][j];
	for(int i = 1; i <= 5; ++i)//找到每行的最大值
	{
		int m_j = 1;//m_j:最大值的纵坐标
		for(int j = 1; j <= 5; ++j)//遍历第i行,j是列号,找到最大值
			if(a[i][j] > a[i][m_j])
				m_j = j;
		x[i][m_j]++;
	}
	for(int j = 1; j <= 5; ++j)//找到每列的最小值
	{
		int m_i = 1;//m_i:最小值的横坐标
		for(int i = 1; i <= 5; ++i)//遍历第j列,i是列号,找最小值
			if(a[i][j] < a[m_i][j])
				m_i = i;
		x[m_i][j]++;
	}
    for(int i = 1; i <= 5; ++i)
        for(int j = 1; j <= 5; ++j)
            if(x[i][j] == 2)//找到鞍点
            {
                cout << i << ' ' << j << ' ' << a[i][j];
                return 0;
            }
    cout << "not found";
    return 0;
}
  • 5
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 2
    评论
能带理论中的鞍点是指在布里渊区中,能带的极值点,它们不是最大值也不是最小值,而是驻点。计算布里渊区中的鞍点,需要进行以下步骤: 1. 构造晶格,并计算出其倒格子。 2. 计算出布里渊区的边界。 3. 在布里渊区内选取一些离散的点,计算出它们对应的能量值。 4. 对能量值进行插值,得到一个平滑的能带曲线。 5. 对能带曲线求导,得到能带曲线的斜率。 6. 找到能带曲线的极值点和鞍点,极值点是斜率为零的点,而鞍点是斜率为零且具有负曲率的点。 MATLAB中可以使用插值函数(如interp1)来插值计算能带曲线,可以使用差分函数(如diff)来计算能带曲线的斜率,可以使用find函数结合逻辑运算来找出能带曲线的极值点和鞍点。例如: ```matlab % 构造晶格和计算倒格子 ... % 计算布里渊区的边界 ... % 在布里渊区内选取一些离散的点,计算出它们对应的能量值 ... % 对能量值进行插值,得到一个平滑的能带曲线 xq = linspace(min(k), max(k), 100); eq = interp1(k, e, xq); % 对能带曲线求导,得到能带曲线的斜率 dq = diff(eq); % 找到能带曲线的极值点和鞍点 max_idx = find(dq(1:end-1) > 0 & dq(2:end) < 0) + 1; min_idx = find(dq(1:end-1) < 0 & dq(2:end) > 0) + 1; saddle_idx = find(dq(1:end-1) < 0 & dq(2:end) < 0) + 1; ``` 其中,`k`是离散的布里渊区点的波矢,`e`是对应的能量值。最后,`max_idx`、`min_idx`和`saddle_idx`分别是能带曲线的极值点、最小值点和鞍点的索引。

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值