我是一个编程的初学者,这几天一直在写一个题目。
给定一个序列X = {a
1, a
2, …, a
n}, X中的一对元素(a
i,a
j) 称为逆对如果 i<j 且 a
i>a
j. 例如,序列{5, 3, 6, 4, 2, 3}具有10个逆对: (5,3), (5,4), (5,2), (5,3), (3,2), (6,4), (6,2), (6,3), (4,2), (4,3).
呵呵!
这个题目对很多人来说,是太简单了阿,可是我却花了5天才找到方法,写出来,我希望有向我一样的初学者,来看看,互相沟通一下,代码很简单,易懂,因为我也是初学者,呵呵!
此算法可以算到10的6次方 的 数。
本人英语学的不好,有很多变量都是用的汉语拼音 例如 " 苹果 " 就写成 pingguo 。
#include <stdlib.h>
#include <stdio.h>
#include <time.h>
#include <math.h>
#include <string>
#include <fstream>
#include <iostream>
#include <stdio.h>
#include <time.h>
#include <math.h>
#include <string>
#include <fstream>
#include <iostream>
using namespace std;
/*******************************************************************************************/
int i,j,a,b;
long int num;
int shuzu[10000000];
int shub[10000000];
int count[3];
/*********************************************************************************************/
//读取数组文件到shuzu[]中 函数名 number_in()
void number_in()
{
clock_t t1,t2;
//
string filename;
printf("/n");
cout<<"请属输入文件名:";
cin>>filename;
printf("/n");
/
t1=clock();
ifstream file(filename.c_str(),ios::in);
if(!file.fail())
{
file>>num;
cout<<"序列元素个数:"<<num<<endl;
for(long i=1;i<=num;i++)
{
file>>shuzu[i];
if(num<=10)cout<<shuzu[i]<<" ";
}
file.close();
printf("/n");
}
t2=clock();
double t=(double)(t2-t1)/CLK_TCK;
printf("读取数组文件时间 %f s/n/n",t);
}
/*********************************************************************************************/
//随即数组 suiji()
int suiji(void)
{
printf("/n随即生成数组:/n/n");
printf("数组元素个数(当元素个数少于等于10时显示数组):");
scanf("%d",&num);
printf("/n/n");
//
clock_t t1,t2;
t1=clock();
for(i=1; i<=num; i++)
{
shuzu[i]=rand() % 1000000;
/*******************************************************************************************/
int i,j,a,b;
long int num;
int shuzu[10000000];
int shub[10000000];
int count[3];
/*********************************************************************************************/
//读取数组文件到shuzu[]中 函数名 number_in()
void number_in()
{
clock_t t1,t2;
//
string filename;
printf("/n");
cout<<"请属输入文件名:";
cin>>filename;
printf("/n");
/
t1=clock();
ifstream file(filename.c_str(),ios::in);
if(!file.fail())
{
file>>num;
cout<<"序列元素个数:"<<num<<endl;
for(long i=1;i<=num;i++)
{
file>>shuzu[i];
if(num<=10)cout<<shuzu[i]<<" ";
}
file.close();
printf("/n");
}
t2=clock();
double t=(double)(t2-t1)/CLK_TCK;
printf("读取数组文件时间 %f s/n/n",t);
}
/*********************************************************************************************/
//随即数组 suiji()
int suiji(void)
{
printf("/n随即生成数组:/n/n");
printf("数组元素个数(当元素个数少于等于10时显示数组):");
scanf("%d",&num);
printf("/n/n");
//
clock_t t1,t2;
t1=clock();
for(i=1; i<=num; i++)
{
shuzu[i]=rand() % 1000000;
}
if(num<=20)
{
for(i=1;i<=num;i++)
{
printf("%6d",shuzu[i]);
}
printf("/n/n");
}
//
t2=clock();
double t=(double)(t2-t1)/CLK_TCK;
printf("生成数组的时间 %f s/n/n",t);
///
return 0;
}
if(num<=20)
{
for(i=1;i<=num;i++)
{
printf("%6d",shuzu[i]);
}
printf("/n/n");
}
//
t2=clock();
double t=(double)(t2-t1)/CLK_TCK;
printf("生成数组的时间 %f s/n/n",t);
///
return 0;
}
/********************************************************************************************************/
//快速排序 quk(int i,int j)
/*int quk(int i,int j)
{
if(i<j)
{
int p=i;int q=j;
int key[2];
key[0]=shuzu[i];
while(i<j)
{
while(i<j && shuzu[j]<=key[0])
{j=j-1;}
shuzu[i]=shuzu[j];
if(i<j)i++;
while(i<j && shuzu[i]>=key[0])
{i=i+1;}
shuzu[j]=shuzu[i];
if(i<j)j--;
if(i==j)
{
shuzu[j]=key[0];
}
}
int r=i;
quk(p,r);
quk(r+1,q);
}
//快速排序 quk(int i,int j)
/*int quk(int i,int j)
{
if(i<j)
{
int p=i;int q=j;
int key[2];
key[0]=shuzu[i];
while(i<j)
{
while(i<j && shuzu[j]<=key[0])
{j=j-1;}
shuzu[i]=shuzu[j];
if(i<j)i++;
while(i<j && shuzu[i]>=key[0])
{i=i+1;}
shuzu[j]=shuzu[i];
if(i<j)j--;
if(i==j)
{
shuzu[j]=key[0];
}
}
int r=i;
quk(p,r);
quk(r+1,q);
}
return 0;
}
}
/********************************************************************************************************/
int binggui(int i,int m,int n)
{
int x=i;
int y=n;
int j;
int k=i;
for(j=m+1;i<=m && j<=n;++k)
{
if(shuzu[i]>shuzu[j])shub[k]=shuzu[i++];
else shub[k]=shuzu[j++];
}
while(i<=m)
{
shub[k++]=shuzu[i++];
}
while(j<=n)
{
shub[k++]=shuzu[j++];
}
for(x=x;x<=y;x++)
{
shuzu[x]=shub[x];
}
return 0;
}
/*******************************************************************************************************/
//找逆对 nidui(int i,int j,int p,int q,int r)
int binggui(int i,int m,int n)
{
int x=i;
int y=n;
int j;
int k=i;
for(j=m+1;i<=m && j<=n;++k)
{
if(shuzu[i]>shuzu[j])shub[k]=shuzu[i++];
else shub[k]=shuzu[j++];
}
while(i<=m)
{
shub[k++]=shuzu[i++];
}
while(j<=n)
{
shub[k++]=shuzu[j++];
}
for(x=x;x<=y;x++)
{
shuzu[x]=shub[x];
}
return 0;
}
/*******************************************************************************************************/
//找逆对 nidui(int i,int j,int p,int q,int r)
int nidui(int i,int j,int p,int q,int r)
{
// quk(p,r);
// quk(r+1,q);
binggui(i,(i+r)/2,r);
binggui(j,(j+q)/2,q);
while(i<=r && j<=q)
{
if(shuzu[i]<=shuzu[j])j++;
else
{
count[0]=count[0]+q-j+1;
i++;
}
}
return 0;
}
/*********************************************************************************************************/
//递归算法 fenzhi(int p,int q)
int fenzhi(int p,int q)
{
//
if(p<q)
{
int r=(p+q)/2;
fenzhi(p,r);
fenzhi(r+1,q);
nidui(p,r+1,p,q,r);
// binggui();
}
return 0;
}
/*********************************************************************************************************/
//主函数
void main()
{
printf("——————————————————————————/n");
printf(" “ 找逆对” /n");
printf("——————————————————————————");
L1: printf("/n/n请选择 :");
printf("/n/n 1,随即生成数组。 2,读取已有数组文件。 3,退出。");
printf("/n/n 选择 : ");
scanf("%d",&a);
printf(" ......................................./n");
///
if(a!=3)
{
if(a==2)number_in();
else if(a==1)suiji();//随即数组
i=1;j=num;
/
clock_t t1,t2;
t1=clock();
/
fenzhi(i,j);
/
t2=clock();
double t=(double)(t2-t1)/CLK_TCK;
printf("逆对查找时间消耗 %f s/n/n",t);
/
printf("逆对个数是 %ld/n/n",count[0]);
printf("——————————————————————————");
goto L1;
L2: printf("/n");
}
{
// quk(p,r);
// quk(r+1,q);
binggui(i,(i+r)/2,r);
binggui(j,(j+q)/2,q);
while(i<=r && j<=q)
{
if(shuzu[i]<=shuzu[j])j++;
else
{
count[0]=count[0]+q-j+1;
i++;
}
}
return 0;
}
/*********************************************************************************************************/
//递归算法 fenzhi(int p,int q)
int fenzhi(int p,int q)
{
//
if(p<q)
{
int r=(p+q)/2;
fenzhi(p,r);
fenzhi(r+1,q);
nidui(p,r+1,p,q,r);
// binggui();
}
return 0;
}
/*********************************************************************************************************/
//主函数
void main()
{
printf("——————————————————————————/n");
printf(" “ 找逆对” /n");
printf("——————————————————————————");
L1: printf("/n/n请选择 :");
printf("/n/n 1,随即生成数组。 2,读取已有数组文件。 3,退出。");
printf("/n/n 选择 : ");
scanf("%d",&a);
printf(" ......................................./n");
///
if(a!=3)
{
if(a==2)number_in();
else if(a==1)suiji();//随即数组
i=1;j=num;
/
clock_t t1,t2;
t1=clock();
/
fenzhi(i,j);
/
t2=clock();
double t=(double)(t2-t1)/CLK_TCK;
printf("逆对查找时间消耗 %f s/n/n",t);
/
printf("逆对个数是 %ld/n/n",count[0]);
printf("——————————————————————————");
goto L1;
L2: printf("/n");
}
}