C++爱心代码

1. #include <stdio.h>
   
2. #include <math.h>
   
3. #include <windows.h>
   
4. #include <tchar.h>
   
5. 
   
6. float f(float x, float y, float z) {
   
7. float a = x * x + 9.0f / 4.0f * y * y + z * z - 1;
   
8. return a * a * a - x * x * z * z * z - 9.0f / 80.0f * y * y * z * z * z;
   
9. }
   
10. 
    
11. float h(float x, float z) {
    
12.   for (float y = 1.0f; y >= 0.0f; y -= 0.001f)
    
13. if (f(x, y, z) <= 0.0f)
    
14. return y;
    
15. return 0.0f;
    
16. }
    
17. 
    
18. int main() {
    
19. HANDLE o = GetStdHandle(STD_OUTPUT_HANDLE);
    
20. _TCHAR buffer[25][80] = { _T(' ') };
    
21. _TCHAR ramp[] = _T("@#$%&*~~");
    
22. 
    
23. for (float t = 0.0f;; t += 0.1f) {
    
24.         int sy = 0;
    
25.         float s = sinf(t);
    
26.         float a = s * s * s * s * 0.2f;
    
27.         for (float z = 1.3f; z > -1.2f; z -= 0.1f) {
    
28.         _TCHAR* p = &buffer[sy++][0];
    
29.         float tz = z * (1.2f - a);
    
30.         for (float x = -1.5f; x < 1.5f; x += 0.05f) {
    
31.         float tx = x * (1.2f + a);
    
32.         float v = f(tx, 0.0f, tz);
    
33.         if (v <= 0.0f) {
    
34.         float y0 = h(tx, tz);
    
35.         float ny = 0.01f;
    
36.         float nx = h(tx + ny, tz) - y0;
    
37.         float nz = h(tx, tz + ny) - y0;
    
38.         float nd = 1.0f / sqrtf(nx * nx + ny * ny + nz * nz);
    
39.         float d = (nx + ny - nz) * nd * 0.5f + 0.5f;
    
40.         *p++ = ramp[(int)(d * 5.0f)];
    
41. }
    
42. else
    
43.         *p++ = ' ';
    
44. }
    
45. }
    
46. 
    
47.         for (sy = 0; sy < 25; sy++) {
    
48.         COORD coord = { 0, sy };
    
49.         SetConsoleCursorPosition(o, coord);
    
50.         WriteConsole(o, buffer[sy], 79, NULL, 0);
    
51. }
    
52. Sleep(33);
    
53. }
    
54. }
    

复制代码