用bitset写了一个高精度程序（待完善）

额外减小

就是模拟计算机的操作来进行运算。

打算做一个1024bit的虚拟操作系统。但是我不想用压位高精，这样就不像计算机了。我还是想用计算机最原始的运算方法来写。所以不要说这个代码常数怎么这么大。

代码如下。到时候慢慢完善。

目前做的是正整数的，应该没问题了，欢迎大家来hack.，但是使用补码应该可以很容易实现整数，浮点数的运算。

但是我现在没有时间做了。以后慢慢来。

1. #include <iostream>
   
2. #include <cstdio>
   
3. #include <typeinfo>
   
4. #include <bitset>
   
5. #include <string>
   
6. #include <utility>
   
7. using namespace std;
   
8. 
   
9. const int _bits=8;
   
10. 
    
11. // 准备写一个高精度类。用 std::bitset
    
12. 
    
13. inline int highest_true(bitset<_bits> x)
    
14. {
    
15.         for(int i=_bits-1;i>=0;i--)
    
16.         {
    
17.                 if(x[i])return i;
    
18.         }
    
19.         return -1;
    
20. }
    
21. 
    
22. inline bool operator>(bitset<_bits> x,bitset<_bits> y)
    
23. {
    
24.         for(int i=_bits-1;i>=0;i--)
    
25.         {
    
26.                 if(x[i]&&!y[i])
    
27.                 return true;
    
28.                 if(!x[i]&&y[i])
    
29.                 return false;
    
30.         }
    
31.         return false;
    
32. }
    
33. 
    
34. inline bool operator<(bitset<_bits> x,bitset<_bits> y)
    
35. {
    
36.         for(int i=_bits-1;i>=0;i--)
    
37.         {
    
38.                 if(x[i]&&!y[i])
    
39.                 return false;
    
40.                 if(!x[i]&&y[i])
    
41.                 return true;
    
42.         }
    
43.         return false;
    
44. }
    
45. 
    
46. inline bool operator>=(bitset<_bits> x,bitset<_bits> y)
    
47. {
    
48.         return !(x<y);
    
49. }
    
50. 
    
51. inline bool operator<=(bitset<_bits> x,bitset<_bits> y)
    
52. {
    
53.         return !(x>y);
    
54. }
    
55. 
    
56. inline bitset<_bits> operator+(bitset<_bits> x,bitset<_bits> y)
    
57. {
    
58.         bool c;
    
59.         bitset<_bits> ans;
    
60.         ans[0]=x[0]^y[0];
    
61.         c=x[0]&&y[0];
    
62.         for(int i=1;i<_bits;i++)
    
63.         {
    
64.                 ans[i]=(x[i]&&y[i]&&c)||(x[i]&&!y[i]&&!c)||(!x[i]&&y[i]&&!c)||(!x[i]&&!y[i]&&c);
    
65.                 c=(x[i]&&y[i]&&c)||(x[i]&&y[i]&&!c)||(!x[i]&&y[i]&&c)||(x[i]&&!y[i]&&c);
    
66.         }
    
67.         return ans;
    
68. }
    
69. 
    
70. inline bitset<_bits> operator-(bitset<_bits> x,bitset<_bits> y)
    
71. {
    
72.         return x+~y+(bitset<_bits>)1;
    
73. }
    
74. 
    
75. inline bitset<_bits> operator*(bitset<_bits> x,bitset<_bits> y)
    
76. {
    
77.         bitset<_bits> tmpx=x,tmpy=y,ans;
    
78.         while(tmpx.any())
    
79.         {
    
80.                 if(tmpx[0])
    
81.                 {
    
82.                         ans=ans+tmpy;
    
83.                 }
    
84.                 tmpx>>=1;
    
85.                 tmpy<<=1;
    
86.         }
    
87.         return ans;
    
88. }
    
89. 
    
90. inline bitset<_bits> operator/(bitset<_bits> x,bitset<_bits> y)
    
91. {
    
92.         int ht1=highest_true(x),ht2=highest_true(y);
    
93.         if(ht1<ht2)return (bitset<_bits>)0;
    
94.         bitset<_bits> tmpx=x,tmpy=y<<(ht1-ht2),ans;
    
95.         while(tmpx>=y)
    
96.         {
    
97.                 ans[0]=tmpx>=tmpy;
    
98.                 if(ans[0])
    
99.                 tmpx=tmpx-tmpy;
    
100.                
     
101.                 ans<<=1;
     
102.                 tmpy>>=1;
     
103.         }
     
104.         ans<<=(highest_true(tmpy)-ht2);
     
105.         return ans;
     
106.        
     
107. }
     
108. 
     
109. inline pair< bitset<_bits> , bitset<_bits> > divmod(bitset<_bits> x,bitset<_bits> y)
     
110. {
     
111.         int ht1=highest_true(x),ht2=highest_true(y);
     
112.         pair< bitset<_bits>,bitset<_bits> > ans_pair;
     
113.         if(ht1<ht2)
     
114.         {
     
115.                 ans_pair.first=0;
     
116.                 ans_pair.second=x;
     
117.                 return ans_pair;
     
118.         }
     
119.         bitset<_bits> tmpx=x,tmpy=y<<(ht1-ht2),div;
     
120.         while(1)
     
121.         {
     
122.                 cout<<tmpx<<endl<<tmpy<<endl<<div<<endl<<endl;
     
123.                
     
124.                 div[0]=(tmpx>=tmpy);
     
125.                
     
126.                
     
127.                 if(div[0])
     
128.                 tmpx=tmpx-tmpy;
     
129.                 if(tmpx<y)
     
130.                 break;
     
131.                 div<<=1;
     
132.                 tmpy>>=1;
     
133.         }
     
134.         div<<=(highest_true(tmpy)-ht2);
     
135.         cout<<highest_true(tmpy)<<' '<<ht2<<endl;
     
136.         ans_pair.first=div;
     
137.         ans_pair.second=tmpx;
     
138.         cout<<endl<<"div: "<<div<<" "<<(highest_true(tmpy)-ht2)<<endl;
     
139.         return ans_pair;
     
140. }
     
141. 
     
142. inline void in_int(bitset<_bits> &x)
     
143. {
     
144.         int a;
     
145.         cin>>a;
     
146.         x=a;
     
147. }
     
148. 
     
149. inline void out_int(bitset<_bits> &x)
     
150. {
     
151.         string str=x.to_string();
     
152.        
     
153.         int ans=0,l=str.length();
     
154.         for(int i=l-1,j=1;i>=0;i--,j*=2)
     
155.         {
     
156.                 if(str[i]=='1')ans+=j;
     
157.         }
     
158.         cout<<"output : "<<ans<<endl;
     
159. }
     
160. 
     
161. int main()
     
162. {
     
163.        
     
164. }

复制代码