#include "iostream.h"
#include "fstream.h"
#include "stdlib.h"
#include "conio.h"

void piksr2(int n, double arr[], double brr[])
{
	int i,j;
	double a,b;
    for (j = 2; j<=n; j++)
	{
        a = arr[j];
        b = brr[j];
        for (i = j - 1; i>=1; i--)
		{
            if (arr[i] <= a)
			{
				break;
			}
            arr[i + 1] = arr[i];
            brr[i + 1] = brr[i];
        }
        if (arr[i] > a)
		{
			i = 0;
		}
        arr[i + 1] = a;
        brr[i + 1] = b;
    }
}

void main()
{
    //program d7r2
    //driver for routine piksr2
    double a[101], b[101];
	int i,j;
	ifstream fin;
	fin.open("d:\\vc常用數(shù)值算法集\\data\\tarray.dat",ios::in|ios::nocreate);
	if (!fin)
	{
		cerr<<"could not open d:\\vc常用數(shù)值算法集\\data\\tarray.dat"<<endl;
		exit(1);
	}
	for (i=1; i<=100; i++)
	{
		fin>>a[i];
	}
	fin.close();
    //generate b-array
    for (i = 1; i<=100; i++)
	{
        b[i] = i;
    }
    //sort a and mix b
    piksr2(100, a, b);
    cout<<"After sorting A and mixing B, array A is:"<<endl;
    for (i = 1; i<=10; i++)
	{
        for (j = 1; j<=10; j++)
		{
			cout.width(7);
			cout.setf(ios::left);
            cout<<a[10 * (i - 1) + j];
        }
        cout<<endl;
    }
	cout<<endl;
    cout<<"...and array B is:"<<endl;
    for (i = 1; i<=10; i++)
	{
        for (j = 1; j<=10; j++)
		{
			cout.width(7);
			cout.setf(ios::left);
            cout<<b[10 * (i - 1) + j];
        }
        cout<<endl;
    }
    cout<<"press any key to continue..."<<endl;
    getch();
    //sort b and mix a
    piksr2(100, b, a);
    cout<<"After sorting A and mixing B, array A is:"<<endl;
    for (i = 1; i<=10; i++)
	{
        for (j = 1; j<=10; j++)
		{
			cout.width(7);
			cout.setf(ios::left);
            cout<<a[10 * (i - 1) + j];
        }
        cout<<endl;
    }
	cout<<endl;
    cout<<"...and array B is:"<<endl;
    for (i = 1; i<=10; i++)
	{
        for (j = 1; j<=10; j++)
		{
			cout.width(7);
			cout.setf(ios::left);
            cout<<b[10 * (i - 1) + j];
        }
        cout<<endl;
    }
}