OpenJudge | 计算鞍点

总时间限制: 1000ms 内存限制: 65536kB

描述

给定一个5*5的矩阵,每行只有一个最大值,每列只有一个最小值,寻找这个矩阵的鞍点。
鞍点指的是矩阵中的一个元素,它是所在行的最大值,并且是所在列的最小值。
例如:在下面的例子中(第4行第1列的元素就是鞍点,值为8 )。
11 3 5 6 9
12 4 7 8 10
10 5 6 9 11
8 6 4 7 2
15 10 11 20 25

输入

输入包含一个5行5列的矩阵

输出

如果存在鞍点,输出鞍点所在的行、列及其值,如果不存在,输出"not found"

样例输入

11 3 5 6 9
12 4 7 8 10
10 5 6 9 11
8  6 4 7 2
15 10 11 20 25

样例输出

4 1 8

思路

  1. 读入矩阵
  2. 找出每行的最大值的位置,将其存入数组line
  3. 找出每列的最小值的位置,将其存入数组row
  4. 遍历linerow,找出鞍点

图例

总时间限制: 1000ms 内存限制: 65536kB

描述

给定一个5*5的矩阵,每行只有一个最大值,每列只有一个最小值,寻找这个矩阵的鞍点。
鞍点指的是矩阵中的一个元素,它是所在行的最大值,并且是所在列的最小值。
例如:在下面的例子中(第4行第1列的元素就是鞍点,值为8 )。
11 3 5 6 9
12 4 7 8 10
10 5 6 9 11
8 6 4 7 2
15 10 11 20 25

输入

输入包含一个5行5列的矩阵

输出

如果存在鞍点,输出鞍点所在的行、列及其值,如果不存在,输出"not found"

样例输入

11 3 5 6 9
12 4 7 8 10
10 5 6 9 11
8  6 4 7 2
15 10 11 20 25

样例输出

4 1 8

思路

  1. 读入矩阵
  2. 找出每行的最大值的位置,将其存入数组line
  3. 找出每列的最小值的位置,将其存入数组row
  4. 遍历linerow,找出鞍点

图例

pk7bT41.png

Code

C++ STL

#include <bits/stdc++.h>
using namespace std;

int main() {
 array<array<int, 5>, 5> a;
 array<pair<int, int>, 5> line, row; // (first, second)
 for(int i = 0; i < 5; i++) {
  for(int j = 0; j < 5; j++) cin >> a[i][j];
 }
 for(int i = 0; i < 5; i++) { // (i, j)
  line[i].first = i, line[i].second = 0;
  for(int j = 0; j < 5; j++) {
   if(a[line[i].first][line[i].second] < a[i][j]) {
    line[i].first = i;
    line[i].second = j;
   }
  }
 }
 for(int i = 0; i < 5; i++) { // (j, i)
  row[i].first = 0, row[i].second = i;
  for(int j = 0; j < 5; j++) {
   if(a[row[i].first][row[i].second] > a[j][i]) {
    row[i].first = j;
    row[i].second = i;
   }
  }
 }
 int flag = 0;
 for(int i = 0; i < 5; i++) {
  for(int j = 0; j < 5; j++) {
   if(row[i].first == line[j].first && line[j].second == row[i].second) {
    cout << row[i].first+1 << " " << row[i].second+1 << " " << a[row[i].first][row[i].second] << endl;
    flag = 1;
   }
  }
 }
 if(!flag) cout << "not found";
}

C

#include <stdio.h>

int main() {
 int a[5][5];
 int line[5][2], row[5][2]; // (first, second)
 for(int i = 0; i < 5; i++) {
  for(int j = 0; j < 5; j++) scanf("%d", &a[i][j]);
 }
 for(int i = 0; i < 5; i++) { // (i, j)
  line[i][0] = i, line[i][1] = 0;
  for(int j = 0; j < 5; j++) {
   if(a[line[i][0]][line[i][1]] < a[i][j]) {
    line[i][0] = i;
    line[i][1] = j;
   }
  }
 }
 for(int i = 0; i < 5; i++) { // (j, i)
  row[i][0] = 0, row[i][1] = i;
  for(int j = 0; j < 5; j++) {
   if(a[row[i][0]][row[i][1]] > a[j][i]) {
    row[i][0] = j;
    row[i][1] = i;
   }
  }
 }
 int flag = 0;
 for(int i = 0; i < 5; i++) {
  for(int j = 0; j < 5; j++) {
   if(row[i][0] == line[j][0] && line[j][1] == row[i][1]) {
    printf("%d %d %d\n", row[i][0]+1, row[i][1]+1, a[row[i][0]][row[i][1]]);
    flag = 1;
   }
  }
 }
 if(!flag) printf("not found");
}

Code

C++ STL

#include <bits/stdc++.h>
using namespace std;

int main() {
 array<array<int, 5>, 5> a;
 array<pair<int, int>, 5> line, row; // (first, second)
 for(int i = 0; i < 5; i++) {
  for(int j = 0; j < 5; j++) cin >> a[i][j];
 }
 for(int i = 0; i < 5; i++) { // (i, j)
  line[i].first = i, line[i].second = 0;
  for(int j = 0; j < 5; j++) {
   if(a[line[i].first][line[i].second] < a[i][j]) {
    line[i].first = i;
    line[i].second = j;
   }
  }
 }
 for(int i = 0; i < 5; i++) { // (j, i)
  row[i].first = 0, row[i].second = i;
  for(int j = 0; j < 5; j++) {
   if(a[row[i].first][row[i].second] > a[j][i]) {
    row[i].first = j;
    row[i].second = i;
   }
  }
 }
 int flag = 0;
 for(int i = 0; i < 5; i++) {
  for(int j = 0; j < 5; j++) {
   if(row[i].first == line[j].first && line[j].second == row[i].second) {
    cout << row[i].first+1 << " " << row[i].second+1 << " " << a[row[i].first][row[i].second] << endl;
    flag = 1;
   }
  }
 }
 if(!flag) cout << "not found";
}

C

#include <stdio.h>

int main() {
 int a[5][5];
 int line[5][2], row[5][2]; // (first, second)
 for(int i = 0; i < 5; i++) {
  for(int j = 0; j < 5; j++) scanf("%d", &a[i][j]);
 }
 for(int i = 0; i < 5; i++) { // (i, j)
  line[i][0] = i, line[i][1] = 0;
  for(int j = 0; j < 5; j++) {
   if(a[line[i][0]][line[i][1]] < a[i][j]) {
    line[i][0] = i;
    line[i][1] = j;
   }
  }
 }
 for(int i = 0; i < 5; i++) { // (j, i)
  row[i][0] = 0, row[i][1] = i;
  for(int j = 0; j < 5; j++) {
   if(a[row[i][0]][row[i][1]] > a[j][i]) {
    row[i][0] = j;
    row[i][1] = i;
   }
  }
 }
 int flag = 0;
 for(int i = 0; i < 5; i++) {
  for(int j = 0; j < 5; j++) {
   if(row[i][0] == line[j][0] && line[j][1] == row[i][1]) {
    printf("%d %d %d\n", row[i][0]+1, row[i][1]+1, a[row[i][0]][row[i][1]]);
    flag = 1;
   }
  }
 }
 if(!flag) printf("not found");
}
  • 7
    点赞
  • 4
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
能带理论中的鞍点是指在布里渊区中,能带的极值点,它们不是最大值也不是最小值,而是驻点。计算布里渊区中的鞍点,需要进行以下步骤: 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`分别是能带曲线的极值点、最小值点和鞍点的索引。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

Mryan2005

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值