#include <stdio.h> #include <time.h> #define SIZE 1000 #define CNT 10 typedef double array[SIZE][SIZE]; void init(array A, array B, array C, int n) { for(int i = 0; i < n; i++) { for(int j = 0; j < n; j++) { A[i][j] = B[i][j] = 1.0; C[i][j] = 0.0; } } } void de_init(array C, int n) { for(int i = 0; i < n; i++) { for(int j = 0; j < n; j++) { C[i][j] = 0.0; } } } //定义 i 行，j 列，k 索引 void ijk(array A, array B, array C, int n) { int i, j, k; double sum; for(i = 0; i < n; i++) { for(j = 0; j < n; j++) { sum = 0.0; for(k = 0; k < n; k++) { sum += A[i][k] * B[k][j]; } C[i][j] += sum; } } } //交换 i 和 j void jik(array A, array B, array C, int n) { int i, j, k; double sum; for(j = 0; j < n; j++) { for(i = 0; i < n; i++) { sum = 0.0; for(k = 0; k < n; k++) { sum += A[i][k] * B[k][j]; } C[i][j] += sum; } } } void ikj(array A, array B, array C, int n) { int i, j, k; double r; for(i = 0; i < n; i++) { for(k = 0; k < n; k++) { r = A[i][k]; for(j = 0; j < n; j++) { C[i][j] += r * B[k][j]; } } } } //交换 k 和 i void kij(array A, array B, array C, int n) { int i, j, k; double r; for(k = 0; k < n; k++) { for(i = 0; i < n; i++) { r = A[i][k]; for(j = 0; j < n; j++) { C[i][j] += r * B[k][j]; } } } } int main() { double t_ijk, t_ikj; clock_t clk_start, clk_end; array A, B, C; init(A, B, C, SIZE); clk_start = clock(); for(int i = 0; i < CNT; i++) { ijk(A, B, C, SIZE); de_init(C, SIZE); } clk_end = clock(); t_ijk = (double)(clk_end - clk_start) / CLOCKS_PER_SEC / 10.0; clk_start = clock(); for(int i = 0; i < CNT; i++) { ikj(A, B, C, SIZE); de_init(C, SIZE); } clk_end = clock(); t_ikj = (double)(clk_end - clk_start) / CLOCKS_PER_SEC / 10.0; printf("ijk:%f\t, ikj:%f\n", t_ijk, t_ikj); }修改此段代码让结果为ijk=3.536100，jik：2.44900

在函数 jik 中，将第二个 for 循环中的 i 和 j 交换即可：

void jik(array A, array B, array C, int n) {
    int i, j, k;
    double sum;
    for(j = 0; j < n; j++) {
        for(i = 0; i < n; i++) { // 交换 i 和 j
            sum = 0.0;
            for(k = 0; k < n; k++) {
                sum += A[i][k] * B[k][j];
            }
            C[i][j] += sum;
        }
    }
}

然后在 main 函数中输出结果：

printf("ijk:%f\t, jik:%f\n", t_ijk, t_jik);

相关问题

代码讲解：#include <stdio.h> #include <stdlib.h> #include <math.h> #define MAX_PAIRS 100 // 定义最大符合条件的数组长度 // 结构体，用于存储符合条件的两个数 typedef struct { int x; int y; } SquarePair; // 函数声明 int check_square_sum(int num, SquarePair *squares, int *cnt); // 检查是否可以表示成两个数的平方和 int main() { int num; SquarePair squares[MAX_PAIRS]; // 定义存储符合条件的两个数的结构体数组 int cnt = 0; // 初始化存储数量的变量 printf("请输入一个正整数："); scanf("%d", &num); // 输入一个正整数 if (check_square_sum(num, squares, &cnt)) { // 如果找到符合条件的两个数 printf("%d = %d^2 + %d^2\n", num, squares[cnt-1].x, squares[cnt-1].y); // 输出最后一个找到的两个数和原数 } else { // 如果找不到符合条件的两个数 printf("NO\n"); // 输出“NO” } return 0; } // 检查是否可以表示成两个数的平方和 int check_square_sum(int num, SquarePair *squares, int *cnt) { int i, j; int square_sum; for (i = 1; i <= sqrt(num); i++) { // 从1到num的平方根遍历一遍 square_sum = i * i; // 计算i的平方 for (j = i; j <= sqrt(num - square_sum); j++) { // 从i到num的平方根遍历一遍 if (i * i + j * j == num) { // 如果num可以表示成i和j的平方和 squares[*cnt].x = i; // 修改结构体中的x值 squares[*cnt].y = j; // 修改结构体中的y值 (*cnt)++; // 数量加1 if (*cnt == MAX_PAIRS) { // 如果达到数组最大长度，则退出循环 return 1; } } } } if (*cnt > 0) { // 如果找到了符合条件的两个数 return 1; } return 0; // 没有找到 }

这段代码是一个简单的平方和问题求解程序。以下是代码的详细讲解：

1. 第1至4行代码是头文件的引用，包括stdio.h、stdlib.h和math.h。这些头文件提供了一些常用的函数和宏定义，例如stdio.h提供了输入输出函数scanf和printf，而math.h提供了计算平方根的函数sqrt。

2. 第6至10行代码定义了一个结构体SquarePair，用于存储符合条件的两个数。结构体包含两个整型变量x和y，分别表示两个数。

3. 第13至22行代码是check_square_sum函数的定义。该函数的主要功能是检查一个正整数是否可以表示成两个数的平方和。函数接收三个参数：一个整数num，一个SquarePair类型的结构体数组squares和一个整型指针cnt。函数会遍历1到num的平方根，计算每个数的平方，并从该数到num的平方根遍历一遍，判断是否可以表示成两个数的平方和。如果找到符合条件的两个数，就将它们存储到squares数组中，并将cnt加1。如果squares数组已经存储了最大数量的符合条件的两个数，就退出循环。如果找到了符合条件的两个数，函数返回1，否则返回0。

4. 第25至36行代码是main函数的定义。该函数的主要功能是从标准输入中读取一个正整数，并调用check_square_sum函数来检查它是否可以表示成两个数的平方和。如果可以，就输出最后一个找到的两个数和原数，否则输出“NO”。在函数开头，定义了一个SquarePair类型的结构体数组squares，用于存储符合条件的两个数，同时初始化了一个整型变量cnt，用于存储squares数组中已经存储的符合条件的两个数的数量。

整个程序的运行逻辑比较简单，主要是调用check_square_sum函数来解决平方和问题。

#pragma GCC optimize ("O3") #pragma pack (16)//所有的存储都是以16个字节为单位的 #include <stdio.h> #include <stdlib.h> #include <string.h> #define tolower(c) (c>='A'&&c<='Z')?c-'A'+'a':c #define DATA 5200000 #define SIZE 1000005 int trie[4200000][26]; typedef struct node { int cnt; int logo; struct node *child[26]; } Node; Node *root; char str[35000000]; typedef struct word { char wor[85]; int cnt; } Word; Word w[300000]; struct node *creat() { Node *Root = (Node *)malloc(sizeof(Node)); Root->logo = 0; Root->cnt = 0; for (int i = 0; i < 26; i++) { Root->child[i] = NULL; } return Root; } void insert(Node *root, char *word, int flag) { struct node *leaf = root; for (int i = 0; word[i] != '\0'; i++) { int index = word[i] - 'a'; if (!leaf->child[index]) { leaf->child[index] = creat(); } leaf = leaf->child[index]; } if (leaf->logo != -1) leaf->logo = flag; leaf->cnt++; } int count = 0; void dfs(Node *leaf, char *word, int level) { if (leaf->logo == 1) { word[level] = '\0'; strcpy(w[count++].wor, word); w[count - 1].cnt = leaf->cnt; } for (int i = 0; i < 26; i++) { if (leaf->child[i]) { word[level] = i + 'a'; dfs(leaf->child[i], word, level + 1); } } } int cmp(const void *p1, const void *p2) { Word *v1, *v2; v1 = (Word *)p1; v2 = (Word *)p2; if (v1->cnt != v2->cnt) return v2->cnt - v1->cnt; else return strcmp(v1->wor, v2->wor); } int main(int argc, char *argv[]) { char s[1024]; int temp; int n, m;//读入n,m; //n = atoi(argv[1]); //m = atoi(argv[2]); scanf("%d%d", &n, &m); //读入stopwords中的元素,并令末序数组值为0，即该单词不计入 root = creat(); FILE *stopwords = fopen("stopwords.txt", "r"); while (fscanf(stopwords, "%s", s) != EOF) { insert(root, s, -1); } int cnt; FILE *article = fopen("article.txt", "r"); cnt = fread(str, sizeof(char), 35000000, article); char word[85]; int w_cnt = 0; for (int i = 0; i < cnt; i++) { char c = tolower(str[i]); if (c >= 'a' && c <= 'z') { word[w_cnt++] = c; } else { word[w_cnt] = '\0'; p = 0; w_cnt = 0; if (strlen(word) > 0) { insert(root, word, 1); } }//对article中的所有单词进行计数 } dfs(root, word, 0); qsort(w, count, sizeof(w[0]), cmp); printf("%s", w[0].cnt); return 0; }

这段代码是一个基于Trie树的单词统计程序，在给定的文章中统计出现频率最高的单词。程序从文件中读入stopwords（不计入统计的常用单词）和article（待统计的文章），使用Trie树统计文章中出现的单词频率，并按照出现频率从高到低输出结果。

具体实现中，程序使用了一个结构体Node表示Trie树中的节点，其中包含了该节点出现次数cnt和单词编号logo。同时，程序定义了一个结构体Word表示单词及其出现次数。程序中的insert函数用于向Trie树中插入单词，dfs函数用于遍历Trie树并统计单词出现次数，cmp函数用于按照出现频率从高到低排序输出结果。

注意，程序中的一些细节需要注意，比如将所有字母转换为小写字母，以及对于长度为0的单词不进行计数等。