所以,实际上我需要的是在排序后保留旧数组的索引。例如,如果我输入[2,4,1,5,7,9,6],那么输出就是[2,0,1,3,6,4,5]。我已经使用了qsort,如果没有重复的元素,它会运行得很好。
如果有重复的元素,有时第一个重复的元素放在最后。例如,如果输入是[5,4,6,5,2,1,3],我想要输出的是[5,4,6,1,0,3,2]。因此,5的索引0放在有索引3的5之前。但是,使用qsort有时使输出[5,4,6,1,3,0,2]。
你能帮我修一下吗?或者我应该创建自己的排序函数?你能帮我创作一下吗?
这是我的代码:
#include <stdlib.h>
int* sortidx(double *X,int n)
{
int *idx,i,j;
int cmp(const void *a,const void *b)
{
return X[*(int*)a]>=X[*(int*)b]?1:-1;
}
idx=(int*)calloc(n,sizeof(int));
for(i=0;i<n;i++)
{
idx[i]=i;
}
qsort(idx,n,sizeof(int),cmp);
return idx;
}发布于 2017-08-13 14:40:33
如果一个元素的值更大,或者它的值相等,并且它的索引更大,您希望一个元素比另一个元素更大。(这就是稳定排序算法背后的想法。)
在本例中,您知道要比较的元素的索引,因此可以轻松地将其添加到比较标准中:
int cmp(const void *a, const void *b)
{
return X[*(int*)a] > X[*(int*)b] ||
(X[*(int*)a] == X[*(int*)b] && *(int*)a > *(int*)b)
?1:-1;
}或者,更容易理解和更正确(因为没有证明a和b是不同的):
int cmp(const void *a, const void *b)
{
int idxa = *(const int*)a, idxb = *(const int*)b;
if (X[idxa] > X[idxb]) return 1;
if (X[idxa] < X[idxb]) return -1;
return idxa - idxb;
}引用参数X的嵌套函数的使用是gcc扩展,可能不适用于其他编译器。标准C库的Gnu实现还包含函数qsort_r,该函数可用于将X传递给比较例程,但编写函数的一种更可移植的方法是使用指针数组而不是索引数组:
int idxcmp(const void *a,const void *b)
{
double *ap = *(double *const*)a, *bp = *(double *const*)b;
if (*ap > *bp) return 1;
if (*ap < *bp) return -1;
return ap - bp;
}
double** sortidx(double *X, size_t n)
{
double **idx = calloc(n, sizeof(double*));
for (size_t i=0; i<n; ++i) idx[i] = X + i;
qsort(idx, n, sizeof(idx[0]), idxcmp);
return idx;
}(如果确实希望返回索引,则可以在返回之前将指针转换为索引。)
https://stackoverflow.com/questions/45661361
复制相似问题