转置一维阵列

Transpose 1 Dimensional Array

本文关键字:阵列 一维 转置      更新时间:2023-10-16

所以我有一个具有N个值的一维数组,其中N是一个完美的正方形。我把这个一维数组想象成二维数组(尽管它不是(。例如,值为int Array = { 0,1,2,3,4,5,6,7,8 } 的数组

那是

int *Array = new int [9];                                                                                                                                                                                                    
for ( int i = 0 ; i < 9 ; i ++ )
         Array[i] = i; // For example

这被打印为

0 1 2
3 4 5 
6 7 8

所以,我想交换一维数组中的位置,这样我就得到了它的转置,。。。

例如。。。

0 3 6
1 4 7
2 5 8

这基本上是相同的一维数组,但值被交换,因此数组现在是int Array = {0,3,6,1,4,7,2,5,8}

如果我将其缩放为1024*1024维的数组,逻辑会如何?

使用n = sqrt(N),您可以尝试一些简单的东西,比如:

for(int i = 0; i < n; ++i)
    for(int j = i+1; j < n; ++j)
        std::swap(Array[n*i + j], Array[n*j + i]);

转置运算对不包括矩阵对角线的上三角形或下三角形执行swap(v[y][x],v[x][y])(假设为上三角形(。

在C一维矢量vc中,v[y][x]对应于vc[y*n+x]。所以你想做vc[y*n+x] = vc[x*n+y]

要交换的元素是x > y所针对的元素。

你最终会做:

for(int y = 0; y < n; ++y)
    for(int x = y+1; x < n; ++x)
        swap(vc[x*n + y], vc[y*n + x]);

你本可以自己想出来的。。。

#include <iostream>
#include <cmath>
using namespace std;
int xyToIndex(const int x, const int y, const int size){
    return x + y * size;
}
int main(){
    int a[] = { 0,1,2,3,4,5,6,7,8 };
    const int size = sqrt(sizeof(a)/sizeof(int));
    //print only
    for(int x = 0;x < size; ++x){
        for(int y = 0; y < size; ++y)
            cout << a[xyToIndex(x,y,size)] << " ";;
        cout << endl;
    }
    //make array
    int b[size*size];
    int index = 0;
    for(int x = 0;x < size; ++x)
        for(int y = 0; y < size; ++y)
             b[index++] = a[xyToIndex(x,y,size)];
    for(int i = 0; i< size * size ; ++i){
        cout << b[i] << " ";
    }
}

您可以交换矩阵中的值,也可以交换后面函数中的解释。

例如,您可以打印一个(j,i(并打印转换姿势,而不是打印(i,j(。

话虽如此,你想做什么检查?如果你看看LAPACK和BLAS,它们的例程采用控制算法的标志来正常解释它们或将它们转换为转置。

不使用交换函数。len是数组的长度。

int i,j;
  N = sqrt(len);    
  int temp[len];
    for(i=0;i<N;i++)
    {   for(j=0;j<N;j++)
        {
            temp[j+(i*N)] = a[(j*N)+i];
        }
    }
static unsigned other(unsigned dim0, unsigned dim1, unsigned index)
{
#if 0
unsigned x0,x1;
x0 = index % dim0 ;
x1 = index / dim0 ;
return x0 * dim1 + x1;
#else
unsigned mod,val;
mod = dim0 * dim1 -1;
val = (index==mod) ? mod: (dim1*index) % mod;
return val;
#endif
}

上面的函数返回索引的"其他"索引(转置矩阵中的索引:=,x和y交换(。dim0和dim1是矩阵的"水平"answers"垂直"大小。#ifdeffed-out部分是天真的实现。在您的情况下,您可以使用初始化(或分析(一维阵列

for (i=0; i < 9; i++)
    arr[i] = other(3, 3, i);