[已解决]高阶矩阵转置算法(百万阶的矩阵)

LordHdy

需要从文件读取矩阵，获取行数和列数，并进行转置。可能不是方阵！有没有算法的空间复杂度和时间复杂度比较好的算法啊

最佳答案

月排行榜 / 总排行榜

风扫地

2018-11-9 12:00:09

思路1.先用MATLAB  textscan出来,然后transpose一下，再写进新的文件，要数据的话就从新文件里面读就可以了，我这个主意你觉得怎么样；
思路2.先获得数据的规模，申请空间，读一个然后算对应位置(应该就交换一下下标)，存入对应空间。再读一个，再存一个，直到读完，存完。

跳转到最佳答案楼层

风扫地

思路1.先用MATLAB  textscan出来,然后transpose一下，再写进新的文件，要数据的话就从新文件里面读就可以了，我这个主意你觉得怎么样；
思路2.先获得数据的规模，申请空间，读一个然后算对应位置(应该就交换一下下标)，存入对应空间。再读一个，再存一个，直到读完，存完。

Charles未晞

矩阵转置还有考虑是不是方阵的么。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。

LordHdy

风扫地 发表于 2018-11-9 12:00
思路1.先用MATLAB  textscan出来,然后transpose一下，再写进新的文件，要数据的话就从新文件里面读就可以了 ...

不能用matlab,老师要求自己用C++写

LordHdy

1. #include <iostream>
   
2. #include <fstream>
   
3. #include <stdlib.h>
   
4. #include <time.h>
   
5. 
   
6. using namespace std;
   
7. 
   
8. void MatrixTranspos(char* infilename, char* outfilename, unsigned long M, unsigned long N, int line)
   
9. {// 原矩阵文件名,  转置矩阵文件名, 行数, 列数, 一次读取行数
   
10. 
    
11.     ifstream in;
    
12.     ofstream out;
    
13. 
    
14.     int buffsize = N*line*3+line*2;               // 缓冲区大小 列数*要读取的行数*3(一个数字占三个字符的位置)+要读取的行数(换行符)
    
15.     char *Matrix = new char[buffsize];            // 分配缓冲区内存
    
16. 
    
17.     in.open(infilename, ios::in|ios::binary);     // 打开原矩阵文件
    
18.     out.open(outfilename, ios::out|ios::app);              // 打开转置矩阵文件
    
19.     if(!in || !out)
    
20.     {
    
21.         cout << "File open error!" << endl;
    
22.         return ;
    
23.     }
    
24. 
    
25.     in.read((char*)Matrix, buffsize);                  // 按块读取数据
    
26.     in.seekg(0, in.beg);
    
27. //    cout.write((char*)Matrix, buffsize);             // 输出缓冲区内数据
    
28. int k = 0;
    
29.     for(int i = 0; i < N; i++)                         // 转置
    
30.     {
    
31.         for(int j = 0; j < line; j++)
    
32.         {
    
33.             out.write((char*)Matrix+3*N*j+2*j, 3);
    
34.             cout.write((char*)Matrix+3*N*j+2*j, 3);
    
35.         }
    
36.         out << "  " << ++k << '\n';
    
37.         cout << endl;
    
38.         Matrix = Matrix+3;
    
39.     }
    
40. 
    
41.     delete[] Matrix;
    
42.     in.close();
    
43.     out.close();
    
44. 
    
45. }
    
46. 
    
47. int main()
    
48. {
    
49.     char *infile = "matrix.txt";
    
50.     char *outfile = "transpos.txt";
    
51. 
    
52. //    // 创建大矩阵文件
    
53. //    ofstream out;
    
54. //    out.open("matrix.txt", ios::out);
    
55. //
    
56. //    srand((int)time(0));
    
57. //    for(int i = 0; i < 3000; i++)
    
58. //    {
    
59. //        for(int j = 0; j <2000; j++)
    
60. //        {
    
61. //            out << setw(3) << (rand()%100);  // 一个数字占三个字节的位置
    
62. //        }
    
63. //        out << endl;
    
64. //    }
    
65. //    out.close();
    
66. 
    
67.     MatrixTranspos(infile, outfile, 3000, 2000, 10);
    
68. 
    
69. 
    
70.     return 0;
    
71. }
    

复制代码

以上是我写的矩阵转置的部分代码，只是用来测试，但是有问题，控制台输出的部分矩阵转置是正确的，但是输出到文件中的矩阵少了几行，有没有哪位可以帮忙看一下怎么回事，该怎么修改呀？
主函数中注释部分是生成原矩阵的代码。

LordHdy

1. #include <iostream>
   
2. #include <fstream>
   
3. #include <stdlib.h>
   
4. #include <time.h>
   
5. 
   
6. using namespace std;
   
7. 
   
8. void MatrixTranspos(char* infilename, char* outfilename, unsigned long M, unsigned long N, int line)
   
9. {// 原矩阵文件名,  转置矩阵文件名, 行数, 列数, 一次读取行数
   
10. 
    
11.     ifstream in;
    
12.     ofstream out;
    
13. 
    
14.     int buffsize = N*line*3+line*2;               // 缓冲区大小 列数*要读取的行数*3(一个数字占三个字符的位置)+要读取的行数(换行符)
    
15.     char *Matrix = new char[buffsize];            // 分配缓冲区内存
    
16.     char *ptr = Matrix;
    
17. 
    
18.     in.open(infilename, ios::in|ios::binary);     // 打开原矩阵文件
    
19.     out.open(outfilename, ios::out|ios::app);              // 打开转置矩阵文件
    
20.     if(!in || !out)
    
21.     {
    
22.         cout << "File open error!" << endl;
    
23.         return ;
    
24.     }
    
25. 
    
26.     in.read((char*)Matrix, buffsize);                  // 按块读取数据
    
27.     in.seekg(0, in.beg);
    
28. //    cout.write((char*)Matrix, buffsize);             // 输出缓冲区内数据
    
29. 
    
30.     for(int i = 0; i < N; i++)                         // 转置
    
31.     {
    
32.         for(int j = 0; j < line; j++)
    
33.         {
    
34.             //out.write((char*)Matrix+3*N*j+2*j, 3);
    
35.             //cout << ptr[3*i+3*N*j+2*j] << ptr[3*i+3*N*j+2*j+1] << ptr[3*i+3*N*j+2*j+2];
    
36.             out << ptr[3*i+3*N*j+2*j] << ptr[3*i+3*N*j+2*j+1] << ptr[3*i+3*N*j+2*j+2];
    
37.             //ptr = ptr+3*N*j+2*j;
    
38.             //cout.write((char*)Matrix+3*N*j+2*j, 3);
    
39.         }
    
40.         out <<'\n';
    
41.         //cout << endl;
    
42.         //cout << "   " << ++k << '\n';
    
43.         //Matrix = Matrix+3;
    
44.         //ptr = Matrix;
    
45.     }
    
46. 
    
47.     delete[] Matrix;
    
48.     in.close();
    
49.     out.close();
    
50. 
    
51. }
    
52. 
    
53. int main()
    
54. {
    
55.     char *infile = "matrix.txt";
    
56.     char *outfile = "transpos.txt";
    
57. 
    
58. //    // 创建大矩阵文件
    
59. //    ofstream out;
    
60. //    out.open("matrix.txt", ios::out);
    
61. //
    
62. //    srand((int)time(0));
    
63. //    for(int i = 0; i < 3000; i++)
    
64. //    {
    
65. //        for(int j = 0; j <2000; j++)
    
66. //        {
    
67. //            out << setw(3) << (rand()%100);  // 一个数字占三个字节的位置
    
68. //        }
    
69. //        out << endl;
    
70. //    }
    
71. //    out.close();
    
72. 
    
73.     MatrixTranspos(infile, outfile, 3000, 2000, 10);
    
74. 
    
75. 
    
76.     return 0;
    
77. }
    

复制代码

这样输出就可以了。
谢谢大家的回复！