搜索
查看: 7|回复: 1

精益求精

[复制链接]
 楼主| 发表于 2020-9-17 15:23:44 | 显示全部楼层 |阅读模式
16啊哈币

  1. #include <stdio.h>
  2. #include <windows.h>
  3. #include <math.h>
  4. #include <graphics.h>
  5. #include <mmsystem.h>
  6. #pragma comment(lib,"winmm.lib")
  7. //定义全局变量
  8. int rosesize = 500;
  9. int h = -250;
  10. //定义结构体
  11. struct DOT {
  12.   double x;
  13.   double y;
  14.   double z;
  15.   double r;
  16.   double g;
  17. };

  18. bool calc(double a,double b,double c,DOT &d) {
  19.   double j, n, o, w, z;

  20.   if (c > 60) { //花柄
  21.     d.x = sin(a * 7)*(13 + 5 / (0.2 + pow(b * 4, 4))) - sin(b) * 50;
  22.     d.y = b*rosesize + 50;
  23.     d.z = 625 + cos(a * 7)*(13 + 5 / (0.2 + pow(b * 4, 4))) + b * 400;
  24.     d.r = a * 1 - b / 2;
  25.     d.g = a;
  26.     return true;
  27.   }
  28.   double A = a * 2 - 1;
  29.   double B = b * 2 - 1;
  30.   if (A*A + B*B < 1) { //绿叶
  31.     if (c > 37) {
  32.       j = (int(c) & 1);
  33.       n = j ? 6 : 4;
  34.       o = 0.5 / (a + 0.01) + cos(b * 125) * 3 - a * 300;
  35.       w = b * h;
  36.       d.x = o * cos(n) + w * sin(n) + j * 610 - 390;
  37.       d.y = o * sin(n) - w * cos(n) + 550 - j * 350;
  38.       d.z = 1180 + cos(B + A) * 99 - j * 300;
  39.       d.r = 0.4 - a * 0.1 + pow(1 - B * B, -h * 6) * 0.15 - a * b * 0.4 + cos(a + b) / 5 +
  40.             pow(cos((o * (a + 1) + (B > 0 ? w : -w)) / 25), 30) * 0.1 * (1 - B * B);
  41.       d.g = o / 1000 + 0.7 - o * w * 0.000003;
  42.       return true;
  43.     }

  44.     if (c > 32) { //花萼
  45.       c = c*1.16 - 0.15;
  46.       o = a * 45 - 20;
  47.       w = b * b * h;
  48.       z = o * sin(c) + w * cos(c) + 620;
  49.       d.x = o * cos(c) - w * sin(c);
  50.       d.y = 28 + cos(B * 0.5) * 99 - b * b * b * 60 - z / 2 - h;
  51.       d.z = z;
  52.       d.r = (b * b * 0.3 + pow((1 - (A * A)), 7)*0.15 + 0.3)*b;
  53.       d.g = b * 0.7;
  54.       return true;
  55.     }

  56. //花朵
  57.     o = A*(2 - b)*(80 - c * 2);
  58.     w = 99 - cos(A) * 120 - cos(b)*(-h - c * 4.9) + cos(pow(1 - b, 7)) * 50 + c * 2;
  59.     z = o * sin(c) + w * cos(c) + 700;
  60.     d.x = o * cos(c) - w * sin(c);
  61.     d.y = B * 99 - cos(pow(b, 7)) * 50 - c / 3 - z / 1.35 + 450;
  62.     d.z = z;
  63.     d.r = (1 - b / 1.2) * 0.9 + a * 0.1;
  64.     d.g = pow((1 - b), 20) / 4 + 0.05;
  65.     return true;
  66.   }

  67.   return false;
  68. }

  69. int main()
  70. {
  71.   short ?*zBuffer;
  72.   int x, y, z, zBufferIndex;
  73.   DOT dot;
  74.   initgraph(640,480);//创建绘图窗口

  75.   /*添加音乐*/
  76.   mciSendString("open BadApple.mp3 alias BadApple", NULL, 0, NULL);
  77.   mciSendString("play BAdApple", NULL, 0, NULL);

  78.   setbkcolor(YELLOW);
  79.   cleardevice();
  80.   zBuffer = new short[rosesize * rosesize];
  81.   memset(zBuffer, 0, sizeof(short) * rosesize * rosesize);
  82.   for (int j = 0; j < 2000; j++)// && !_kbhit()
  83.   {
  84.     for (int i = 0; i < 10000; i++)
  85.       if (calc(double(rand()) / RAND_MAX, double(rand()) / RAND_MAX, rand() % 46 / 0.74, dot))
  86.       {
  87.         z = int(dot.z + 0.5);
  88.         x = int(dot.x*rosesize / z - h + 0.5);
  89.         y = int(dot.y*rosesize / z - h + 0.5);
  90.         if (y > rosesize) continue;
  91.         zBufferIndex = y*rosesize + x;
  92.         if (!zBuffer[zBufferIndex] || zBuffer[zBufferIndex] > z)
  93.         {
  94.           zBuffer[zBufferIndex] = z;
  95. //画点
  96.           int r = ~int((dot.r*h));
  97.           ?
  98.           if (r < 0)r = 0;
  99.           ?
  100.           if (r > 255)r = 255;

  101.           int g = ~int((dot.g*h));
  102.           if (g < 0)g = 0;
  103.           if (g > 255)g = 255;

  104.           int b = ~int((dot.r*dot.r*-80));
  105.           if (b < 0)b = 0;
  106.           if (b > 255)b = 255;

  107.           putpixel(x + 50, y - 20, RGB(r, g, b));
  108.         }
  109.       }
  110.     Sleep(100);
  111.   }
  112.   delete[]zBuffer;
  113.   getchar();
  114.   closegraph();
  115.   return 0;
  116. }
复制代码

 楼主| 发表于 2020-9-17 15:24:24 | 显示全部楼层
就这样就这样就这样就这样就这样就这样就这样就这样就这样就这样
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

广播台
特别关注
快速回复 返回顶部 返回列表