主頁 > 後端開發 > bfprt-線性查找演算法-(topK問題)

bfprt-線性查找演算法-(topK問題)

2020-09-13 20:35:04 後端開發

bfprt-線性查找演算法
百度百科: bfprt
bfprt 流程:
   呼叫 selectPivot 得到 劃分軸.
   呼叫 partition 按斬訓分軸進行區間劃分(小于劃分軸區域,等于劃分軸區域,大于劃分軸區域),
   根據 k 所在的區域進行處理:
	  在等于區域: 直接回傳等于區間任意元素即可.
	  在小于區域: 遞回呼叫 bfprt 處理 小于區域的第k大或小的元素.
	  在大于區域: 遞回呼叫 bfprt 處理 大于區域的第k大或小的元素.

selectPivot 流程:
   陣列區域中每5個連續元素成為一組, 最后不足5個元素的也成為一組,
   每組進行組內排序, 并把 每組的中位數 組成中位數陣列,
   繼續遞回呼叫自身 或 bfprt 處理中位數陣列,
   直至中位數陣列只有一個元素時回傳(即是中位陣列的中位數).

partition 流程:
   根據劃分軸進行劃分: 小于的放左邊,等于的放中間,大于的放右邊,
   最后回傳等于區域的左右邊界(閉區間).




創建3個檔案:bfprt.h、bfprtArray.c、bfprtArrayTest.c,




bfprtArray.h
#ifndef BFPRT_H_
#define BFPRT_H_

// 功能: 比較兩個資料.
// 引數: a(資料a), b(資料b).
// 回傳: a>b回傳正數, a<b回傳負數, 否則回傳0.
// 注意: 其實回傳值可以是任意int型別值, 但是為了統一規范使用所以才強制回傳值.
typedef int( CompareFunc )( const void *a, const void *b );

// 功能: 陣列某一閉區間內的第k小或大的數.
// 引數: a(陣列首地址), left(左閉邊界), right(右閉邊界), k(第k小或大), cmp(比較函式).
// 回傳: 陣列指定區間的第k小或大的元素值.
// 注意: 當 a=NULL 或 cmp=NULL 或 區間范圍表示錯誤 時, 將錯誤退出程式.
extern int bfprtArray( int32_t a[], int32_t left, int32_t right, int32_t k, CompareFunc *cmp );

#endif

bfprtArray.c
#include <stdio.h>
#include <stdlib.h>
#include <stdint.h>
#include "bfprtArray.h"

// 功能: 列印錯誤資訊后就錯誤退出程式.
// 引數: expression(錯誤判斷運算式), message(需列印的錯誤資訊).
// 回傳: 無.
// 注意: 當運算式 expression 為真時, 才觸發.
#define ERROR_EXIT( expression, message )                                    \
if( (expression) ) {                                                         \
	fprintf( stderr, "\nerror location: file = %s, func = %s, line = %d.\n", \
	                       __FILE__, __func__, __LINE__ );                   \
	fprintf( stderr, "error  message: %s%s.\n\a",                            \
	                       (message) != NULL ? (message) : __func__,         \
		                   (message) != NULL ? "" : " function error" );     \
	exit( EXIT_FAILURE );                                                    \
}

// 功能: 陣列的兩個元素進行互換.
// 引數: a(陣列首地址), i(陣列下標), j(陣列下標).
// 回傳: 無.
// 注意: 沒有副作用的宏, C99標準: #define f ({...}).
#define SWAP( a, i, j ) ({            \
	int32_t i1i1i = (i), j1j1j = (j); \
	int32_t t1t1t = *((a) + i1i1i);   \
	*((a) + i1i1i) = *((a) + j1j1j);  \
	*((a) + j1j1j) = t1t1t;           \
})

#define MAX( a, b ) ((a) >= (b) ? (a) : (b))
#define MIN( a, b ) ((a) <= (b) ? (a) : (b))
  • 版本一: 每次 selectPivot 程序動態生成新的中位數陣列.
    // 功能: 陣列某一閉區間進行排序.
    // 引數: a(陣列首地址), left(左閉邊界), right(右閉邊界), cmp(比較函式).
    // 回傳: 無.
    static void insertSort( int32_t a[], int32_t left, int32_t right, CompareFunc *cmp ) {
    	int32_t i = 0, j = 0, t = 0;
    
    	for( int32_t i = left; i <= right; ++i ) {
    		t = a[i];
    		for( j = i - 1; j >= left && cmp( &t, &a[j] ) < 0; --j ) {
    			a[j + 1] = a[j];
    		}
    		a[j + 1] = t;
    	}
    }
    
    // 功能: 陣列某一閉區間的bfprt劃分軸.
    // 引數: a(陣列首地址), left(左閉邊界), right(右閉邊界), cmp(比較函式).
    // 回傳: bfprt劃分軸.
    // 注意: 每5個數為一組進行排序,并得到一個中位數陣列,
    //       中位數陣列繼續bfprt演算法程序, 直至陣列中位數陣列大小為1.
    static int32_t selectPivot( int32_t a[], int32_t left, int32_t right, CompareFunc *cmp ) {
    	int32_t n = right - left + 1;
    	int32_t i = 0, sortSize = 0, median = 0;
    	int32_t *medianA = NULL, mi = 0;
    
    	if( n == 1 ) {
    		return a[left];
    	}
    
    	medianA = malloc( sizeof(*medianA) * (n / 5 + MIN( n % 5, 1 )) );
    
    	// 每5個數為一組進行插入排序, 并把每組的上中位數構成一個中位數陣列.
    	for( i = left; i <= right; i += 5 ) {
    		sortSize = MIN( 5, right - i + 1 );
    		insertSort( a, i, i + sortSize - 1, cmp );
    		medianA[mi++] = a[i + (sortSize - 1) / 2];
    	}
    
    	median = mi > 1 ? bfprtArray( medianA, 0, mi - 1, (mi + 1) / 2, cmp ) : medianA[0];
    	free( medianA );
    
    	return median;
    }
    
    // 功能: 陣列某一閉區間按照指定的劃分軸進行區間劃分.
    // 引數: a(陣列首地址), left(左閉邊界), right(右閉邊界),
    //       pivot(劃分軸), boundary(傳入傳出引數), cmp(比較函式).
    // 回傳: 無, 實際是使用了傳入傳出引數 boundary 進行結果回傳.
    // 注意: boundary[0] 表示等于區域的左閉邊界, boundary[1] 表示等于區域的右閉邊界.
    static void partition( int32_t a[], int32_t left, int32_t right, int32_t pivot, int32_t boundary[2], CompareFunc *cmp ) {
    	int32_t l = left - 1, c = left, r = right + 1;
    
    	while( c < r ) {
    		int32_t cmpRes = cmp( &a[c], &pivot );
    		if( cmpRes == 0 ) {
    			++c;
    		} else if( cmpRes < 0 ) {
    			SWAP( a, c++, ++l );
    		} else {
    			SWAP( a, c, --r );
    		}
    	}
    	boundary[0] = l + 1; // 等于區域的左閉邊界.
    	boundary[1] = r - 1; // 等于區域的右閉邊界.
    }
    
    // 功能: 陣列某一閉區間內的第k小或大的數.
    // 引數: a(陣列首地址), left(左閉邊界), right(右閉邊界), k(第k小或大).
    // 回傳: 陣列指定區間的第k小或大的元素值.
    // 注意: 當 a=NULL 或 cmp=NULL 或 區間范圍表示錯誤 時, 將錯誤退出程式.
    int32_t bfprtArray( int32_t a[], int32_t left, int32_t right, int32_t k, CompareFunc *cmp ) {
    	int32_t pivot = 0, boundary[2] = {-1,-1};
    
    	ERROR_EXIT( a == NULL || cmp == NULL, "NullPointerException" );
    	ERROR_EXIT( left < 0 || left > right, "IndexOutOfBoundsException" );
    	ERROR_EXIT( k < 1, "ParemeterError: k" );
    
    	pivot = selectPivot( a, left, right, cmp );        // 選取劃分軸.
    	partition( a, left, right, pivot, boundary, cmp ); // 根據劃分軸進行區間劃分.
    
    	if( k - 1 < boundary[0] ) {
    		return bfprtArray( a, left, boundary[0] - 1, k, cmp );
    	} else if( k - 1 > boundary[1] ) {
    		return bfprtArray( a, boundary[1] + 1, right, k, cmp );
    	}
    
    	return a[boundary[0]];
    }
    
  • 版本二: 優化上一版本, 優化掉中位數陣列, 并且在 selectPivot 程序中, 不在遞回呼叫bfprt求中位數陣列的中位數, 而是遞回呼叫本身.
    // 功能: 陣列某一閉區間進行排序.
    // 引數: a(陣列首地址), left(左閉邊界), right(右閉邊界), cmp(比較函式).
    // 回傳: 無.
    static void insertSort( int32_t a[], int32_t left, int32_t right, CompareFunc *cmp ) {
    	int32_t i = 0, j = 0, t = 0;
    
    	for( int32_t i = left; i <= right; ++i ) {
    		t = a[i];
    		for( j = i - 1; j >= left && cmp( &t, &a[j] ) < 0; --j ) {
    			a[j + 1] = a[j];
    		}
    		a[j + 1] = t;
    	}
    }
    
    // 功能: 陣列某一閉區間的bfprt劃分軸.
    // 引數: a(陣列首地址), left(左閉邊界), right(右閉邊界), cmp(比較函式).
    // 回傳: bfprt劃分軸.
    // 注意: 每5個數為一組進行排序,并得到一個中位數陣列,
    //       中位數陣列繼續bfprt演算法程序, 直至陣列中位數陣列大小為1.
    static int32_t selectPivot( int32_t a[], int32_t left, int32_t right, CompareFunc *cmp ) {
    	int32_t i = 0, sortSize = 0, mi = 0;
    
    	if( left == right) {
    		return a[left];
    	}
    	// 每5個數為一組進行排序, 并把每組的上中位數構成一個中位數陣列.
    	for( i = left; i <= right; i += 5 ) {
    		sortSize = MIN( 5, right - i + 1 );
    		insertSort( a, i, i + sortSize - 1, cmp );
    		// 將每組上中位數交換至陣列前部分.
    		SWAP( a, left + mi++, i + (sortSize - 1) / 2 );
    	}
    
    	return  mi > 1 ? selectPivot( a, left, left + mi - 1, cmp ) : a[left];
    }
    
    // 功能: 陣列某一閉區間按照指定的劃分軸進行區間劃分.
    // 引數: a(陣列首地址), left(左閉邊界), right(右閉邊界),
    //       pivot(劃分軸), boundary(傳入傳出引數), cmp(比較函式).
    // 回傳: 無, 實際是使用了傳入傳出引數 boundary 進行結果回傳.
    // 注意: boundary[0] 表示等于區域的左閉邊界, boundary[1] 表示等于區域的右閉邊界.
    static void partition( int32_t a[], int32_t left, int32_t right, int32_t pivot, int32_t boundary[2], CompareFunc *cmp ) {
    	int32_t l = left - 1, c = left, r = right + 1;
    
    	while( c < r ) {
    		int32_t cmpRes = cmp( &a[c], &pivot );
    		if( cmpRes == 0 ) {
    			++c;
    		} else if( cmpRes < 0 ) {
    			SWAP( a, c++, ++l );
    		} else {
    			SWAP( a, c, --r );
    		}
    	}
    	boundary[0] = l + 1; // 等于區域的左閉邊界.
    	boundary[1] = r - 1; // 等于區域的右閉邊界.
    }
    
    // 功能: 陣列某一閉區間內的第k小或大的數.
    // 引數: a(陣列首地址), left(左閉邊界), right(右閉邊界), k(第k小或大).
    // 回傳: 陣列指定區間的第k小或大的元素值.
    // 注意: 當 a=NULL 或 cmp=NULL 或 區間范圍表示錯誤 時, 將錯誤退出程式.
    int32_t bfprtArray( int32_t a[], int32_t left, int32_t right, int32_t k, CompareFunc *cmp ) {
    	int32_t pivot = 0, boundary[2] = {-1,-1};
    
    	ERROR_EXIT( a == NULL || cmp == NULL, "NullPointerException" );
    	ERROR_EXIT( left < 0 || left > right, "IndexOutOfBoundsException" );
    	ERROR_EXIT( k < 1, "ParemeterError: k" );
    
    	pivot = selectPivot( a, left, right, cmp );        // 選取劃分軸.
    	partition( a, left, right, pivot, boundary, cmp ); // 根據劃分軸進行區間劃分.
    
    	if( k - 1 < boundary[0] ) {
    		return bfprtArray( a, left, boundary[0] - 1, k, cmp );
    	} else if( k - 1 > boundary[1] ) {
    		return bfprtArray( a, boundary[1] + 1, right, k, cmp );
    	}
    
    	return a[boundary[0]];
    }
    
  • 版本三: 優化上一版本, 使 selectPivot 回傳下標索引而不是元素值, 而該 partition 程序將會與快排的 partition 一毛一樣.
    // 功能: 陣列某一閉區間進行排序.
    // 引數: a(陣列首地址), left(左閉邊界), right(右閉邊界), cmp(比較函式).
    // 回傳: 無.
    static void insertSort( int32_t a[], int32_t left, int32_t right, CompareFunc *cmp ) {
    	int32_t i = 0, j = 0, t = 0;
    
    	for( int32_t i = left; i <= right; ++i ) {
    		t = a[i];
    		for( j = i - 1; j >= left && cmp( &t, &a[j] ) < 0; --j ) {
    			a[j + 1] = a[j];
    		}
    		a[j + 1] = t;
    	}
    }
    
    // 功能: 陣列某一閉區間的bfprt劃分軸.
    // 引數: a(陣列首地址), left(左閉邊界), right(右閉邊界), cmp(比較函式).
    // 回傳: bfprt劃分軸.
    // 注意: 每5個數為一組進行排序,并得到一個中位數陣列,
    //       中位數陣列繼續bfprt演算法程序, 直至陣列中位數陣列大小為1.
    static int32_t selectPivot( int32_t a[], int32_t left, int32_t right, CompareFunc *cmp ) {
    	int32_t i = 0, sortSize = 0, mi = 0;
    
    	if( left == right ) {
    		return left;
    	}
    	// 每5個數為一組進行排序, 并把每組的上中位數構成一個中位數陣列.
    	for( i = left; i <= right; i += 5 ) {
    		sortSize = MIN( 5, right - i + 1 );
    		insertSort( a, i, i + sortSize - 1, cmp );
    		// 將每組上中位數交換至陣列前部分.
    		SWAP( a, left + mi++, i + (sortSize - 1) / 2 );
    	}
    
    	return  mi > 1 ? selectPivot( a, left, left + mi - 1, cmp ) : left;
    }
    
    // 功能: 陣列某一閉區間按照指定的劃分軸進行區間劃分.
    // 引數: a(陣列首地址), left(左閉邊界), right(右閉邊界),
    //       pivot(劃分軸), boundary(傳入傳出引數), cmp(比較函式).
    // 回傳: 無, 實際是使用了傳入傳出引數 boundary 進行結果回傳.
    // 注意: boundary[0] 表示等于區域的左閉邊界, boundary[1] 表示等于區域的右閉邊界.
    static void partition( int32_t a[], int32_t left, int32_t right, int32_t boundary[2], CompareFunc *cmp ) {
    	int32_t l = left - 1, c = left, r = right;
    
    	while( c < r ) {
    		int32_t cmpRes = cmp( &a[c], &a[right] );
    		if( cmpRes == 0 ) {
    			++c;
    		} else if( cmpRes < 0 ) {
    			SWAP( a, c++, ++l );
    		} else {
    			SWAP( a, c, --r );
    		}
    	}
    	SWAP( a, r, right );
    	boundary[0] = l + 1; // 等于區域的左閉邊界.
    	boundary[1] = r;     // 等于區域的右閉邊界.
    }
    
    // 功能: 陣列某一閉區間內的第k小或大的數.
    // 引數: a(陣列首地址), left(左閉邊界), right(右閉邊界), k(第k小或大).
    // 回傳: 陣列指定區間的第k小或大的元素值.
    // 注意: 當 a=NULL 或 cmp=NULL 或 區間范圍表示錯誤 時, 將錯誤退出程式.
    int32_t bfprtArray( int32_t a[], int32_t left, int32_t right, int32_t k, CompareFunc *cmp ) {
    	int32_t pivot = 0, boundary[2] = {-1,-1};
    
    	ERROR_EXIT( a == NULL || cmp == NULL, "NullPointerException" );
    	ERROR_EXIT( left < 0 || left > right, "IndexOutOfBoundsException" );
    	ERROR_EXIT( k < 1, "ParemeterError: k" );
    
    	pivot = selectPivot( a, left, right, cmp ); // 選取劃分軸.
    	SWAP( a, pivot, right );                    // 劃分軸位置與最后位置進行交換.
    	partition( a, left, right, boundary, cmp ); // 根據劃分軸進行區間劃分.
    
    	if( k - 1 < boundary[0] ) {
    		return bfprtArray( a, left, boundary[0] - 1, k, cmp );
    	} else if( k - 1 > boundary[1] ) {
    		return bfprtArray( a, boundary[1] + 1, right, k, cmp );
    	}
    
    	return a[boundary[0]];
    }
    
  • 版本四: 修改上一版本, 使 insertSort 回傳排序區間中位數的下標索引, 而 selectPivot 不必遞回處理直至中位數陣列只有一個元素, 而是當不多于5個元素時,直接排序回傳.
    // 功能: 陣列某一閉區間進行排序.
    // 引數: a(陣列首地址), left(左閉邊界), right(右閉邊界), cmp(比較函式).
    // 回傳: 無.
    static int32_t insertSort( int32_t a[], int32_t left, int32_t right, CompareFunc *cmp ) {
    	int32_t i = 0, j = 0, t = 0;
    
    	for( int32_t i = left; i <= right; ++i ) {
    		t = a[i];
    		for( j = i - 1; j >= left && cmp( &t, &a[j] ) < 0; --j ) {
    			a[j + 1] = a[j];
    		}
    		a[j + 1] = t;
    	}
    
    	return left + (right - left) / 2;
    }
    
    // 功能: 陣列某一閉區間的bfprt劃分軸.
    // 引數: a(陣列首地址), left(左閉邊界), right(右閉邊界), cmp(比較函式).
    // 回傳: bfprt劃分軸.
    // 注意: 每5個數為一組進行排序,并得到一個中位數陣列,
    //       中位數陣列繼續bfprt演算法程序, 直至陣列中位數陣列大小為1.
    static int32_t selectPivot( int32_t a[], int32_t left, int32_t right, CompareFunc *cmp ) {
    	int32_t i = 0, sortSize = 0, mi = 0;
    
    	if( right - left + 1 <= 5 ) {
    		return insertSort( a, left, right, cmp );
    	}
    	// 每5個數為一組進行排序, 并把每組的上中位數構成一個中位數陣列.
    	for( i = left; i <= right; i += 5 ) {
    		sortSize = MIN( 5, right - i + 1 );
    		insertSort( a, i, i + sortSize - 1, cmp );
    		// 將每組上中位數交換至陣列前部分.
    		SWAP( a, left + mi++, i + (sortSize - 1) / 2 );
    	}
    	if( mi <= 5 ) {
    		return insertSort( a, left, left + mi - 1, cmp );
    	}
    
    	return selectPivot( a, left, left + mi - 1, cmp );
    }
    
    // 功能: 陣列某一閉區間按照指定的劃分軸進行區間劃分.
    // 引數: a(陣列首地址), left(左閉邊界), right(右閉邊界),
    //       pivot(劃分軸), boundary(傳入傳出引數), cmp(比較函式).
    // 回傳: 無, 實際是使用了傳入傳出引數 boundary 進行結果回傳.
    // 注意: boundary[0] 表示等于區域的左閉邊界, boundary[1] 表示等于區域的右閉邊界.
    static void partition( int32_t a[], int32_t left, int32_t right, int32_t boundary[2], CompareFunc *cmp ) {
    	int32_t l = left - 1, c = left, r = right;
    
    	while( c < r ) {
    		int32_t cmpRes = cmp( &a[c], &a[right] );
    		if( cmpRes == 0 ) {
    			++c;
    		} else if( cmpRes < 0 ) {
    			SWAP( a, c++, ++l );
    		} else {
    			SWAP( a, c, --r );
    		}
    	}
    	SWAP( a, r, right );
    	boundary[0] = l + 1; // 等于區域的左閉邊界.
    	boundary[1] = r;     // 等于區域的右閉邊界.
    }
    
    // 功能: 陣列某一閉區間內的第k小或大的數.
    // 引數: a(陣列首地址), left(左閉邊界), right(右閉邊界), k(第k小或大).
    // 回傳: 陣列指定區間的第k小或大的元素值.
    // 注意: 當 a=NULL 或 cmp=NULL 或 區間范圍表示錯誤 時, 將錯誤退出程式.
    int32_t bfprtArray( int32_t a[], int32_t left, int32_t right, int32_t k, CompareFunc *cmp ) {
    	int32_t pivot = 0, boundary[2] = {-1,-1};
    
    	ERROR_EXIT( a == NULL || cmp == NULL, "NullPointerException" );
    	ERROR_EXIT( left < 0 || left > right, "IndexOutOfBoundsException" );
    	ERROR_EXIT( k < 1, "ParemeterError: k" );
    
    	pivot = selectPivot( a, left, right, cmp ); // 選取劃分軸.
    	SWAP( a, pivot, right );                    // 劃分軸位置與最后位置進行交換.
    	partition( a, left, right, boundary, cmp ); // 根據劃分軸進行區間劃分.
    
    	if( k - 1 < boundary[0] ) {
    		return bfprtArray( a, left, boundary[0] - 1, k, cmp );
    	} else if( k - 1 > boundary[1] ) {
    		return bfprtArray( a, boundary[1] + 1, right, k, cmp );
    	}
    
    	return a[boundary[0]];
    }
    

bfprtArrayTest.c
實作對數器.
#include <stdio.h>
#include <stdlib.h>
#include <stdint.h>
#include <time.h>
#include "bfprtArray.c"


#define MAX( a, b ) ((a) >= (b) ? (a) : (b))


// 功能: 對陣列的某一區間內的元素值進行隨機化.
// 引數: a(陣列首地址), left(左閉區間), right(右閉區間), v(最大隨機值).
// 回傳: 無.
// 注意: 當v是正數/負數/零時,隨機值的區間分別為[-v, v]/[]/[].
static void randomArray( int32_t a[], int32_t left, int32_t right, int32_t v ) {
	v -= v != INT32_MAX ? 0 : 1;
	while( left <= right ) {
		a[left++] = rand() % (v + 1) - rand() % (v + 1);
	}
}

// 功能: 將源陣列的某一區間內的元素值拷貝至目的陣列中.
// 引數: a(源陣列首地址), left(左閉區間), right(右閉區間), b(目的陣列首地址).
// 回傳: 無.
// 注意: 無.
static void cloneArray( const int32_t a[], int32_t left, int32_t right, int32_t b[] ) {
	while( left <= right ) {
		b[left] = a[left];
		++left;
	}
}

// 功能: 將陣列的某一區間內的元素值送入到檔案流中.
// 引數: a(陣列首地址), left(左閉區間), right(右閉區間), fp(檔案流指標).
// 回傳: 無.
// 注意: 無.
static void printArray( const int32_t a[], int32_t left, int32_t right, FILE *fp ) {
	fprintf( fp, "[" );
	while( left <= right ) {
		fprintf( fp, "%5d%s", a[left], left != right ? ", " : "" );
		++left;
	}
	fprintf( fp, "]\n" );
}

// 功能: 資料比較.
// 引數: a(資料1), b(資料2).
// 回傳: a<b回傳負數, a>b回傳正數, 否則回傳0.
// 注意: 無.
static int cmp( const void *a, const void *b ) {
	#if 1
	return *(int32_t *) a - *(int32_t *) b;
	#else
	return *(int32_t *) b - *(int32_t *) a;
	#endif
}

// 功能: 絕對正確的方法.
// 引數: a(陣列首地址), left(左閉區間), right(右閉區間).
// 回傳: 無.
// 注意: 使用穩定的庫函式方法 或 非常笨的方法(如暴力回圈,暴力遞回搜索).
static int correct( int32_t a[], int32_t left, int32_t right, int32_t k ) {
	qsort( &a[left], right - left + 1, sizeof(*a), cmp );

	return a[k - 1];
}

// 功能: 待測驗的方法.
// 引數: a(陣列首地址), left(左閉區間), right(右閉區間).
// 回傳: 無.
// 注意: 使用技巧性的高效的正確性未知的方法.
static int test( int32_t a[], int32_t left, int32_t right, int32_t k ) {
	return bfprtArray( a, left, right, k, cmp );
}

int main( int argc, char *argv[] ) {
	const int32_t times = 7654321; // 測驗的總測驗.
	const int32_t maxSize = 123;   // 隨機的最大長度.
	int32_t size = 0, left = 0, right = 0;
	int32_t *a = NULL, *t = NULL, *c = NULL;
	int32_t i = 0;

	srand( time( NULL ) );
	a = malloc( sizeof(*a) * maxSize );
	c = malloc( sizeof(*c) * maxSize );
	t = malloc( sizeof(*t) * maxSize );

	for( i = 0; i < times; ++i ) {
		size = rand() % maxSize + 1; // 隨機化陣列長度, 取值區間(0,maxSize].

		randomArray( a, 0, size - 1, 321 );
		cloneArray( a, 0, size - 1, c );
		cloneArray( a, 0, size - 1, t );

		left = 0;
		right = size - 1;

		int k = rand() % size + 1; // 取值區間(0,size].

		int cr = correct( c, left, right, k ); // 絕對正確的方法.
		int tr = test( t, left, right, k );    // 待測驗的方法.
		if( cr != tr ) {
			fprintf( stderr, "source  :" );
			printArray( a, 0, size - 1, stderr );
			fprintf( stderr, "correct : %d\n", cr );
			fprintf( stderr, "test    : %d\n", tr );
			fprintf( stderr, " 本次測驗不通過!\n\a" );
			exit( EXIT_FAILURE );
		}
	}

	free( t );
	free( c );
	free( a );

	printf( "總共 %d 次測驗且全部通過!\n", times );

	return EXIT_SUCCESS;
}

bfprtArrayTest.sh
# !/bin/bash

for(( i = 1; i <= 21; ++i )) do
	printf "%02d" ${i}
	echo -n ______________
	./bfprtTest
done

參考: https://my.oschina.net/dfbb/blog/393346
參考: https://en.academic.ru/dic.nsf/enwiki/301561

重點參考人物: 左神.



轉載請註明出處,本文鏈接:https://www.uj5u.com/houduan/27319.html

標籤:C

上一篇:關于atmega32u4的CDC模擬串口程式

下一篇:解決c語言編譯運行速度慢的問題

標籤雲
其他(157675) Python(38076) JavaScript(25376) Java(17977) C(15215) 區塊鏈(8255) C#(7972) AI(7469) 爪哇(7425) MySQL(7132) html(6777) 基礎類(6313) sql(6102) 熊猫(6058) PHP(5869) 数组(5741) R(5409) Linux(5327) 反应(5209) 腳本語言(PerlPython)(5129) 非技術區(4971) Android(4554) 数据框(4311) css(4259) 节点.js(4032) C語言(3288) json(3245) 列表(3129) 扑(3119) C++語言(3117) 安卓(2998) 打字稿(2995) VBA(2789) Java相關(2746) 疑難問題(2699) 细绳(2522) 單片機工控(2479) iOS(2429) ASP.NET(2402) MongoDB(2323) 麻木的(2285) 正则表达式(2254) 字典(2211) 循环(2198) 迅速(2185) 擅长(2169) 镖(2155) 功能(1967) .NET技术(1958) Web開發(1951) python-3.x(1918) HtmlCss(1915) 弹簧靴(1913) C++(1909) xml(1889) PostgreSQL(1872) .NETCore(1853) 谷歌表格(1846) Unity3D(1843) for循环(1842)

熱門瀏覽
  • 【C++】Microsoft C++、C 和匯編程式檔案

    ......

    uj5u.com 2020-09-10 00:57:23 more
  • 例外宣告

    相比于斷言適用于排除邏輯上不可能存在的狀態,例外通常是用于邏輯上可能發生的錯誤。 例外宣告 Item 1:當函式不可能拋出例外或不能接受拋出例外時,使用noexcept 理由 如果不打算拋出例外的話,程式就會認為無法處理這種錯誤,并且應當盡早終止,如此可以有效地阻止例外的傳播與擴散。 示例 //不可 ......

    uj5u.com 2020-09-10 00:57:27 more
  • Codeforces 1400E Clear the Multiset(貪心 + 分治)

    鏈接:https://codeforces.com/problemset/problem/1400/E 來源:Codeforces 思路:給你一個陣列,現在你可以進行兩種操作,操作1:將一段沒有 0 的區間進行減一的操作,操作2:將 i 位置上的元素歸零。最終問:將這個陣列的全部元素歸零后操作的最少 ......

    uj5u.com 2020-09-10 00:57:30 more
  • UVA11610 【Reverse Prime】

    本人看到此題沒有翻譯,就附帶了一個自己的翻譯版本 思考 這一題,它的第一個要求是找出所有 $7$ 位反向質數及其質因數的個數。 我們應該需要質數篩篩選1~$10^{7}$的所有數,這里就不慢慢介紹了。但是,重讀題,我們突然發現反向質數都是 $7$ 位,而將它反過來后的數字卻是 $6$ 位數,這就說明 ......

    uj5u.com 2020-09-10 00:57:36 more
  • 統計區間素數數量

    1 #pragma GCC optimize(2) 2 #include <bits/stdc++.h> 3 using namespace std; 4 bool isprime[1000000010]; 5 vector<int> prime; 6 inline int getlist(int ......

    uj5u.com 2020-09-10 00:57:47 more
  • C/C++編程筆記:C++中的 const 變數詳解,教你正確認識const用法

    1、C中的const 1、區域const變數存放在堆疊區中,會分配記憶體(也就是說可以通過地址間接修改變數的值)。測驗代碼如下: 運行結果: 2、全域const變數存放在只讀資料段(不能通過地址修改,會發生寫入錯誤), 默認為外部聯編,可以給其他源檔案使用(需要用extern關鍵字修飾) 運行結果: ......

    uj5u.com 2020-09-10 00:58:04 more
  • 【C++犯錯記錄】VS2019 MFC添加資源不懂如何修改資源宏ID

    1. 首先在資源視圖中,添加資源 2. 點擊新添加的資源,復制自動生成的ID 3. 在解決方案資源管理器中找到Resource.h檔案,編輯,使用整個專案搜索和替換的方式快速替換 宏宣告 4. Ctrl+Shift+F 全域搜索,點擊查找全部,然后逐個替換 5. 為什么使用搜索替換而不使用屬性視窗直 ......

    uj5u.com 2020-09-10 00:59:11 more
  • 【C++犯錯記錄】VS2019 MFC不懂的批量添加資源

    1. 打開資源頭檔案Resource.h,在其中預先定義好宏 ID(不清楚其實ID值應該設定多少,可以先新建一個相同的資源項,再在這個資源的ID值的基礎上遞增即可) 2. 在資源視圖中選中專案資源,按F7編輯資源檔案,按 ID 型別 相對路徑的形式添加 資源。(別忘了先把檔案拷貝到專案中的res檔案 ......

    uj5u.com 2020-09-10 01:00:19 more
  • C/C++編程筆記:關于C++的參考型別,專供新手入門使用

    今天要講的是C++中我最喜歡的一個用法——參考,也叫別名。 參考就是給一個變數名取一個變數名,方便我們間接地使用這個變數。我們可以給一個變數創建N個參考,這N + 1個變數共享了同一塊記憶體區域。(參考型別的變數會占用記憶體空間,占用的記憶體空間的大小和指標型別的大小是相同的。雖然參考是一個物件的別名,但 ......

    uj5u.com 2020-09-10 01:00:22 more
  • 【C/C++編程筆記】從頭開始學習C ++:初學者完整指南

    眾所周知,C ++的學習曲線陡峭,但是花時間學習這種語言將為您的職業帶來奇跡,并使您與其他開發人員區分開。您會更輕松地學習新語言,形成真正的解決問題的技能,并在編程的基礎上打下堅實的基礎。 C ++將幫助您養成良好的編程習慣(即清晰一致的編碼風格,在撰寫代碼時注釋代碼,并限制類內部的可見性),并且由 ......

    uj5u.com 2020-09-10 01:00:41 more
最新发布
  • Rust中的智能指標:Box<T> Rc<T> Arc<T> Cell<T> RefCell<T> Weak

    Rust中的智能指標是什么 智能指標(smart pointers)是一類資料結構,是擁有資料所有權和額外功能的指標。是指標的進一步發展 指標(pointer)是一個包含記憶體地址的變數的通用概念。這個地址參考,或 ” 指向”(points at)一些其 他資料 。參考以 & 符號為標志并借用了他們所 ......

    uj5u.com 2023-04-20 07:24:10 more
  • Java的值傳遞和參考傳遞

    值傳遞不會改變本身,參考傳遞(如果傳遞的值需要實體化到堆里)如果發生修改了會改變本身。 1.基本資料型別都是值傳遞 package com.example.basic; public class Test { public static void main(String[] args) { int ......

    uj5u.com 2023-04-20 07:24:04 more
  • [2]SpinalHDL教程——Scala簡單入門

    第一個 Scala 程式 shell里面輸入 $ scala scala> 1 + 1 res0: Int = 2 scala> println("Hello World!") Hello World! 檔案形式 object HelloWorld { /* 這是我的第一個 Scala 程式 * 以 ......

    uj5u.com 2023-04-20 07:23:58 more
  • 理解函式指標和回呼函式

    理解 函式指標 指向函式的指標。比如: 理解函式指標的偽代碼 void (*p)(int type, char *data); // 定義一個函式指標p void func(int type, char *data); // 宣告一個函式func p = func; // 將指標p指向函式func ......

    uj5u.com 2023-04-20 07:23:52 more
  • Django筆記二十五之資料庫函式之日期函式

    本文首發于公眾號:Hunter后端 原文鏈接:Django筆記二十五之資料庫函式之日期函式 日期函式主要介紹兩個大類,Extract() 和 Trunc() Extract() 函式作用是提取日期,比如我們可以提取一個日期欄位的年份,月份,日等資料 Trunc() 的作用則是截取,比如 2022-0 ......

    uj5u.com 2023-04-20 07:23:45 more
  • 一天吃透JVM面試八股文

    什么是JVM? JVM,全稱Java Virtual Machine(Java虛擬機),是通過在實際的計算機上仿真模擬各種計算機功能來實作的。由一套位元組碼指令集、一組暫存器、一個堆疊、一個垃圾回收堆和一個存盤方法域等組成。JVM屏蔽了與作業系統平臺相關的資訊,使得Java程式只需要生成在Java虛擬機 ......

    uj5u.com 2023-04-20 07:23:31 more
  • 使用Java接入小程式訂閱訊息!

    更新完微信服務號的模板訊息之后,我又趕緊把微信小程式的訂閱訊息給實作了!之前我一直以為微信小程式也是要企業才能申請,沒想到小程式個人就能申請。 訊息推送平臺🔥推送下發【郵件】【短信】【微信服務號】【微信小程式】【企業微信】【釘釘】等訊息型別。 https://gitee.com/zhongfuch ......

    uj5u.com 2023-04-20 07:22:59 more
  • java -- 緩沖流、轉換流、序列化流

    緩沖流 緩沖流, 也叫高效流, 按照資料型別分類: 位元組緩沖流:BufferedInputStream,BufferedOutputStream 字符緩沖流:BufferedReader,BufferedWriter 緩沖流的基本原理,是在創建流物件時,會創建一個內置的默認大小的緩沖區陣列,通過緩沖 ......

    uj5u.com 2023-04-20 07:22:49 more
  • Java-SpringBoot-Range請求頭設定實作視頻分段傳輸

    老實說,人太懶了,現在基本都不喜歡寫筆記了,但是網上有關Range請求頭的文章都太水了 下面是抄的一段StackOverflow的代碼...自己大修改過的,寫的注釋挺全的,應該直接看得懂,就不解釋了 寫的不好...只是希望能給視頻網站開發的新手一點點幫助吧. 業務場景:視頻分段傳輸、視頻多段傳輸(理 ......

    uj5u.com 2023-04-20 07:22:42 more
  • Windows 10開發教程_編程入門自學教程_菜鳥教程-免費教程分享

    教程簡介 Windows 10開發入門教程 - 從簡單的步驟了解Windows 10開發,從基本到高級概念,包括簡介,UWP,第一個應用程式,商店,XAML控制元件,資料系結,XAML性能,自適應設計,自適應UI,自適應代碼,檔案管理,SQLite資料庫,應用程式到應用程式通信,應用程式本地化,應用程式 ......

    uj5u.com 2023-04-20 07:22:35 more