清华大学06年复试机试题小解

第一题：

 

计算n<=1000的阶乘。

 

由于计算机中能表达的最大整数就是long long ，只有64位，当n超过13时就会溢出。

 

因此对于小规模处理时，使用整形数组便可以。模拟手工乘法。

 

代码如下：

 

#include <iostream> #include <cstdlib> using namespace std; const int N = 100000; typedef struct{ int d[N];//用数组存储数据 int last;//数据长度 }bignum; void mutiply(bignum *ans,int n){ //把每位都乘以n for(int i = 0;i <= ans->last;i++) ans->d[i] *= n; //处理进位 int j; for(j = 0;(j <=ans->last)||(ans->d[j] != 0);j++){ int tmp = ans->d[j]; ans->d[j] = ans->d[j] % 10; ans->d[j+1] += (tmp/10); } ans->last = j-1; } int main(){ int n; bignum a; while(cin>>n){ //初始化 a.last = 0; memset(a.d ,0,sizeof(a.d)); a.d[0] = 1; for(int i = 1;i <= n;i++) mutiply(&a,i); for(int i = a.last;i >= 0;i--) cout<<a.d[i]; cout<<endl; } return 1; }

 

第二题：

 

是典型的求最大子向量问题，在编程叽珠第八章中详细讨论了这个问题，有兴趣的可以好好研究一下：

 

暂时做了三种不同效率的算法，供对比

 

 

#include<iostream> #include<fstream> #include<cstdio> using namespace std; const int N = 100000; int a[N]; int n; //o(n^2) __int64 fun1(int &start,int &end) { __int64 sum = 0; __int64 max = a[0]; for(int i = 0; i < n;i++){ sum = 0; for(int j = i;j < n;j++){ sum += a[j]; if(sum > max){ max = sum; start = i; end = j; } } } return max; } //o(nlogn),分治法 __int64 fun2(int left,int right) { if(left > right) return 0; if(left == right) return a[left]; int m = (left+right)/2; __int64 lmax = a[m]; __int64 rmax = a[m+1]; __int64 sum = 0; __int64 max = 0; for(int i = m;i >= left;i--){ sum += a[i]; if(sum > lmax) lmax = sum; } sum = 0; for(int j = m+1;j <= right;j++){ sum += a[j]; if(sum > rmax) rmax = sum; } max = lmax+rmax; __int64 ltmp = fun2(left,m); __int64 rtmp = fun2(m+1,right); if(max < ltmp) max = ltmp; if(max < rtmp) max = rtmp; return max; } //o(n),DP __int64 fun(int &start,int &end) { int bstart,bend; //用于记录当前临时的子向量的起始位置 __int64 b = 0; //记录当前子向量和 __int64 sum = a[0]; bstart = bend = 0; for(int i = 0;i < n;i++){ if(b>0){ //若之前的最大子向量和为正,可继续加 b += a[i]; bend = i; }else{ //若之前的最大子向量和已为负,则需要重新开始记录 b = a[i]; bstart = i; bend = i; } if(b > sum){ //找到更优的解 start = bstart; end = bend; sum = b; } } return sum; } int main() { fstream f("t.txt"); memset(a,0,sizeof(a)); int start,end; while(f>>n){ for(int i = 0;i < n;i++) f>>a[i]; printf("%I64d/n",fun1(start,end)); cout<<"start = "<<start<<" end = "<<end<<endl; printf("%I64d/n",fun2(0,n-1)); printf("%I64d/n",fun(start,end)); cout<<"start = "<<start<<" end = "<<end<<endl; } return 0; }  

第三题与上海交大04年第三题完全相同，详情见：

 

http://blog.csdn.net/HereChaos/archive/2010/03/11/5369833.aspx