sicily 6415. linear correlation

Description

向量u=(a1, a2, ...., an)和向量v=(b1, b2, ....., bn)为Rⁿ中的两个向量，判断它们是否线性相关。

Input

输入包括多组测试用例，每组测试用例包括三行。

第一行是一个数n(2<=n<=10)，表示向量的维数。

第二行是向量u的n个元素。

第三行是向量v的n个元素。

当输入的n为0时结束。

Output

 对于每组测试用例输出一行。

如果u和v线性相关，则输出Yes，否则输出No。

这道题的测试用例似乎很糟糕，一开始写了N个漏洞多多的程序居然都给AC（包括判定一个含零向量的向量组线性无关，判定一个每项都是另一向量不同倍的向量为线性相关，某一项有零就判断出错之类的漏洞程序），还需要自己在检查过程中不断修改，最后成品是这个样子，就自己设计测试的用例（包括零向量，含零非零向量，小数，负数）来说暂时没发现有什么大问题

懒得再写一遍题解了，直接存一下提交作业时写的英文函数注释

/*
 * algorithm:
 *  1. create an array to store each scaler that transform each entry of v into
 *     corresponding entry of u. if one or two of the entires is 0, then the 
 *     scaler would be stored as 0;
 *  2. bubble sort the array of scalers;
 *  3. if the smallest nonzero scaler equals to the greatst one, then this 
 *     nonzero scaler could transform the entire v into u, thus u and v are
 *     linearly correlated. otherwise they are not linearly correlated.
 */

#include<stdio.h>
#include<math.h>

int correlation( double u[], int sizeu, double v[], int sizev );
void bubbleSort ( double array[], int size );

int main()
{
    int n, i = 0, j = 0;
    double u[10];
    double v[10];
    
    while ( scanf("%d", &n ) && ( n != 0 ) )
    {
        for ( i = 0; i < n; i++ )
        {
            scanf("%lf", &u[i]); 
        }
        
        for ( i = 0; i < n; i++ )
        {
            scanf("%lf", &v[i]);
        }
        
        if ( correlation(u, n, v, n) == 1 )
        {
            printf( "%s\n", "Yes" );
        }
        else
        {
            printf( "%s\n", "No" );
        } 
    }
    
    return 0;
}
 
int correlation( double u[], int sizeu, double v[], int sizev )
{
    int i;
    int n = sizeu;
    double scaler[n];
    
    for ( i = 0; i < n; i++ )
    {
        if ( fabs( u[i] ) < 0.000001 || fabs( v[i] ) < 0.000001 )
        {
            scaler[i] = 0.0;
        }
        else
        {
            scaler[i] = u[i] / v[i];
        }
    }
    
    bubbleSort( scaler, n );
    
    i = 0;
    while ( fabs(scaler[i]) < 0.000001 && i < n )
    {
        i++;
    }
    
    if ( fabs( scaler[n-1] - scaler[i] ) < 0.000001 )
    {
        return 1;
    }
    else
    {
        return 0;
    }
}

 
void bubbleSort ( double array[], int size )
{
    int i, j;
    double temp;
    
    for ( i = 0; i < size; i++ )
    {
        for( j = 0; j < size - 1 - i; j++ )
        {
            if ( array[j] > array[j+1] )
            {
                temp = array[j];
                array[j] = array[j+1];
                array[j+1] = temp;
            }
        }
    }
}