#include <stdio.h> #include<stdlib.h> #include <math.h> const int max_iterations = 128; const float stop_threshold = 0.01f; const float grad_step = 0.01f; const float clip_far = 10.0f; const float PI = 3.14159265359f; const float PI2 = 6.28318530718f; const float DEG_TO_RAD = PI / 180.0f; typedef struct { float x, y; } vec2; typedef struct { float x, y, z; } vec3; typedef struct { float m[9]; } mat3; const vec3 light_pos = { 20.0f, 60.0f, 20.0f }; float min(float a, float b) { return a < b ? a : b; } float max(float a, float b) { return a > b ? a : b; } float clamp(float f, float a, float b) { return max(min(f, b), a); } vec2 make2(float x, float y) { vec2 r = { x, y }; return r; } vec2 add2(vec2 a, vec2 b) { vec2 r = { a.x + b.x, a.y + b.y }; return r; } vec2 sub2(vec2 a, vec2 b) { vec2 r = { a.x - b.x, a.y - b.y }; return r; } float dot2(vec2 a, vec2 b) { return a.x * b.x + a.y * b.y; } float length2(vec2 v) { return sqrt(dot2(v, v)); } vec3 make3(float x, float y, float z) { vec3 r = { x, y, z }; return r; } vec3 add3(vec3 a, vec3 b) { vec3 r = { a.x + b.x, a.y + b.y, a.z + b.z }; return r; } vec3 sub3(vec3 a, vec3 b) { vec3 r = { a.x - b.x, a.y - b.y, a.z - b.z }; return r; } vec3 mul3(vec3 a, vec3 b) { vec3 r = { a.x * b.x, a.y * b.y, a.z * b.z }; return r; } vec3 scale3(vec3 v, float s) { vec3 r = { v.x * s, v.y * s, v.z *s }; return r; } float dot3(vec3 a, vec3 b) { return a.x * b.x + a.y * b.y + a.z * b.z; } float length3(vec3 v) { return sqrt(dot3(v, v)); } vec3 normalize3(vec3 v) { return scale3(v, 1.0f / length3(v)); } vec3 mul(mat3 m, vec3 v) { return make3( m.m[0] * v.x + m.m[3] * v.y + m.m[6] * v.z, m.m[1] * v.x + m.m[4] * v.y + m.m[7] * v.z, m.m[2] * v.x + m.m[5] * v.y + m.m[8] * v.z); } mat3 rotationXY(float x, float y) { vec2 c = { cos(x), cos(y) }, s = { sin(x), sin(y) }; mat3 m = { c.y , 0.0f, -s.y, s.y * s.x, c.x, c.y * s.x, s.y * c.x, -s.x, c.y * c.x }; return m; } float opI(float d1, float d2) { return max(d1, d2); } float opU(float d1, float d2) { return min(d1, d2); } float opS(float d1, float d2) { return max(-d1, d2); } float sdPetal(vec3 p, float s) { p = add3(mul3(p, make3(0.8f, 1.5f, 0.8f)), make3(0.1f, 0.0f, 0.0f)); vec2 q = make2(length2(make2(p.x, p.z)), p.y); float lower = length2(q) - 1.0f; lower = opS(length2(q) - 0.97f, lower); lower = opI(lower, q.y); float upper = length2(sub2(q, make2(s, 0.0f))) + 1.0f - s; upper = opS(upper, length2(sub2(q, make2(s, 0.0f))) + 0.97f - s); upper = opI(upper, -q.y); upper = opI(upper, q.x - 2.0f); float region = length3(sub3(p, make3(1.0f, 0.0f, 0.0f))) - 1.0f; return opI(opU(upper, lower), region); } float map(vec3 p) { float d = 1000.0f, s = 2.0f; mat3 r = rotationXY(0.1f, PI2 * 0.618034f); r.m[0] *= 1.08f; r.m[1] *= 1.08f; r.m[2] *= 1.08f; r.m[3] *= 0.995f; r.m[4] *= 0.995f; r.m[5] *= 0.995f; r.m[6] *= 1.08f; r.m[7] *= 1.08f; r.m[8] *= 1.08f; for (int i = 0; i < 21; i++) { d = opU(d, sdPetal(p, s)); p = mul(r, p); p = add3(p, make3(0.0, -0.02, 0.0)); s *= 1.05f; } return d; } vec3 gradient(vec3 pos) { const vec3 dx = { grad_step, 0.0, 0.0 }; const vec3 dy = { 0.0, grad_step, 0.0 }; const vec3 dz = { 0.0, 0.0, grad_step }; return normalize3(make3( map(add3(pos, dx)) - map(sub3(pos, dx)), map(add3(pos, dy)) - map(sub3(pos, dy)), map(add3(pos, dz)) - map(sub3(pos, dz)))); } float ray_marching(vec3 origin, vec3 dir, float start, float end) { float depth = start; for (int i = 0; i < max_iterations; i++) { float dist = map(add3(origin, scale3(dir, depth))); if (dist < stop_threshold) return depth; depth += dist * 0.3; if ( depth >= end) return end; } return end; } float shading(vec3 v, vec3 n, vec3 eye) { vec3 ev = normalize3(sub3(v, eye)); vec3 vl = normalize3(sub3(light_pos, v)); float diffuse = dot3(vl, n) * 0.5f + 0.5f; vec3 h = normalize3(sub3(vl, ev)); float rim = pow(1.0f - max(-dot3(n, ev), 0.0f), 2.0f) * 0.15f; float ao = clamp(v.y * 0.5f + 0.5f, 0.0f, 1.0f); return (diffuse + rim) * ao; } vec3 ray_dir(float fov, vec2 pos) { vec3 r = { pos.x, pos.y, -tan((90.0f - fov * 0.5f) * DEG_TO_RAD) }; return normalize3(r); } float f(vec2 fragCoord) { vec3 dir = ray_dir(45.0f, fragCoord); vec3 eye = { 0.0f, 0.0f, 4.5f }; //mat3 rot = rotationXY(-1.0f, 1.0f);mat3 rot = rotationXY(-1.0f, 1.0f); // mat3 rot = rotationXY(-0.6f, 1.0f); dir = mul(rot, dir); eye = mul(rot, eye); float depth = ray_marching(eye, dir, 0.0f, clip_far); vec3 pos = add3(eye, scale3(dir, depth)); if (depth >= clip_far) return 0.0f; else return shading(pos, gradient(pos), eye); } int main() {system("color 0C"); // puts(" 33[91m"); for (int y = 5; y < 50; y++) { for (int x = 0; x < 140; x++) { //putchar(" .,-:;+=*pier#@"[(int)(f(make2((x / 160.0f - 0.5f) * 2.0f, (y / 80.0f - 0.5f) * -2.0f)) * 12.0f)]);putchar(" .,-:;+=*#@"[(int)(f(make2((x / 140.0f - 0.5f) * 2.0f, (y / 50.0f - 0.5f) * -2.0f)) * 12.0f)]);if(y==46&&x==50){puts("一个闲得蛋疼的程序员所写! 三月十四日凌晨n");}} putchar('n'); }system("pause");return 0; } //#include <stdlib.h> //#include<fstream> //#include "svpng.h" // //#define C(x, y) (img + ((y) * 6 + 3) * a + (x) * 6 + 3) //#define W(o) c = C(x, y) + o; c[0] = 255; c[1] = c[2] = d % 512 < 256 ? d % 512 : 511 - d % 512 //#define V(x, y) (x) >= 0 && (x) < w && (y) >= 0 && (y) < h && !*C(x, y) && m(x, y) // //static const int w = 256, h = 256, a = w * 2 + 1, b = h * 2 + 1; //static unsigned char img[a * b * 3], *c; // //int m(int x, int y) { // float u = 2.4f * x / w - 1.2f, v = -2.4f * y / h + 1.3f, w = u * u + v * v - 1; // return w * w * w - u * u * v * v * v < 0; //} // //void f(int x, int y, int d) { // W(0); // for (int j = rand() % 24 * 4, i = 0; i < 4; i++) // switch ("ABCDABDCACBDACDBADBCADCBBACDBADCBCADBCDABDACBDCACABDCADBCBADCBDACDABCDBADABCDACBDBACDBCADCABDCBA"[i + j]) { // case 'A': if (V(x - 1, y)) { W( -3); f(x - 1, y, d + 1); } break; // case 'B': if (V(x + 1, y)) { W( 3); f(x + 1, y, d + 1); } break; // case 'C': if (V(x, y - 1)) { W(a * -3); f(x, y - 1, d + 1); } break; // case 'D': if (V(x, y + 1)) { W(a * 3); f(x, y + 1, d + 1); } break; // } //} // //int main() { // f(w / 2, h / 2, 0); // svpng(fopen("heart.png", "wb"), a, b, img, 0); //}
相关知识
情人节表白代码:浪漫玫瑰花,俘获少女心
c语言表白程序源码玫瑰花,小伙子用C语言写出绽放的玫瑰花,成功表白C++代码女神!...
c语言玫瑰花代码devc++
python仿抖音表白神器资源
C语言情人节玫瑰花代码
c语言制作超炫流星雨表白
C#桌面玫瑰花程序情人节玫瑰花礼物资源
vbs玫瑰花表白代码带音乐?简单的病毒编程代码
计算机编程玫瑰花,c语言实现玫瑰花的方法
comblockengine windows
网址: C语言程序windows窗口显示玫瑰花,表白,装B神器 https://m.huajiangbk.com/newsview958629.html
上一篇: 520最实用的两个Python表 |
下一篇: 表白神器玫瑰花 |