Tweetable Mathematical Art 代码框架

最新推荐文章于 2022-12-31 23:15:35 发布
最新推荐文章于 2022-12-31 23:15:35 发布
阅读量189 收藏 1
点赞数
原文链接：https://my.oschina.net/2bit/blog/498193
版权

为什么80%的码农都做不了架构师？>>>   hot3.png

昨晚在微信公众号看到 《有没有一段代码，让你觉得人类的智慧也可以璀璨无比》里面介绍了一个叫《Tweetable Mathematical Art》的代码游戏，看完我激动不已，这些人都太NB了。

我迫不及待要试一下，可惜该代码游戏原来的框架生成的图片是ppm格式的，ppm格式的图片在Windwos上打开比较麻烦，所以我索性写了一个生成bmp格式的代码框架。

在这个代码框架下，如果有新idea，只需在body.hpp中加入一个新类，然后用REG注册，就可以生成新的图片了。

这个代码框架支持gcc和msvc.

pic.hpp

#ifndef _PIC_HPP_
#define _PIC_HPP_

#include <cstdio>
#include <cstring>

class pic {
public:
  pic() {}
  virtual ~pic() {}

  //  ----x--------->
  //  |
  //  y
  //  |
  //  v
  virtual void set(unsigned int x, unsigned int y, unsigned char r, unsigned char g, unsigned char b) = 0;

  bool save(const char * fname) {
    char * new_name = new char[strlen(fname) + strlen(ext) + 2];
    sprintf(new_name, "%s.%s", fname, ext);

    FILE * f = fopen(new_name, "wb");
    delete [] new_name;

    if(f == NULL)
      return false;

    fwrite(p, 1, sz, f);
    fclose(f);
    return true;
  }

protected:
  unsigned char * p;
  unsigned int sz;
  unsigned int w;
  unsigned int h;
  char ext[8];
};

class ppm : public pic {
public:
  ppm(unsigned int width, unsigned int height) {
    w = width & 0xFFF;
    h = height & 0xFFF;

    char hdr[64];
    sprintf(hdr, "P6\n%u %u\n255\n", w, h);

    hdr_offset = strlen(hdr);
    sz = hdr_offset + (w * 3 * h);

    p = new unsigned char[sz];
    memcpy(p, hdr, hdr_offset);
    memset(p + hdr_offset, 0xFF, sz - hdr_offset);

    strcpy(ext, "ppm");
  }

  ~ppm() { delete [] p; }

  void set(unsigned int x, unsigned int y, unsigned char r, unsigned char g, unsigned char b) {
    if(x < w && y < h) {
      unsigned int offset = hdr_offset + (w * 3 * y) + (x * 3);
      p[offset] = r;
      p[offset + 1] = g;
      p[offset + 2] = b;
    }
  }

private:
  unsigned int hdr_offset;
};

class bmp : public pic {
public:
  bmp(unsigned int width, unsigned int height) {
    w = width & 0xFFF;
    h = height & 0xFFF;
    w3 = (w * 3 + 3) & ~3;
    sz = 54 + (w3 * h);

    p = new unsigned char[sz];
    memset(p, 0, 54);
    memset(p + 54, 0xFF, sz - 54);

    p[0] = 'B';
    p[1] = 'M';
    p[2] = sz;
    p[3] = sz >> 8;
    p[4] = sz >> 16;
    p[5] = sz >> 24;
    p[10] = 54;
    p[14] = 40;
    p[18] = w;
    p[19] = w >> 8;
    p[22] = h;
    p[23] = h >> 8;
    p[26] = 1;
    p[28] = 24;
    p[34] = sz - 54;
    p[35] = (sz - 54) >> 8;
    p[36] = (sz - 54) >> 16;
    p[37] = (sz - 54) >> 24;

    strcpy(ext, "bmp");
  }

  ~bmp() { delete [] p; }

  void set(unsigned int x, unsigned int y, unsigned char r, unsigned char g, unsigned char b) {
    if(x < w && y < h) {
      unsigned int offset = 54 + w3 * (h - 1 - y) + x * 3;
      p[offset] = b;
      p[offset + 1] = g;
      p[offset + 2] = r;
    }
  }

private:
  unsigned int w3;
};

#endif

core.hpp

#ifndef _CORE_HPP_
#define _CORE_HPP_

#include <cstdio>
#include <cstdlib>
#include <cstring>

struct core {
  core():DIM(1024) {}
  virtual ~core() {}

  virtual unsigned char RD(int i, int j) = 0;
  virtual unsigned char GR(int i, int j) = 0;
  virtual unsigned char BL(int i, int j) = 0;

  unsigned int DIM;
};

typedef core * (*CTOR)();

struct name_ctor {
  char name[64];
  CTOR ctor;
};

static name_ctor nc[128];
static unsigned int nc_num = 0;

static int reg_nc(const char * name, CTOR ctor){
  if(strlen(name) >= sizeof(nc[0].name) || nc_num >= sizeof(nc)/sizeof(nc[0])) {
    printf("Error at %s line %d\n", __FILE__, __LINE__);
    exit(1);
  }

  strcpy(nc[nc_num].name, name);
  nc[nc_num].ctor = ctor;
  nc_num++;

  return 0;
}

#define REG(NAME) \
static core * NAME##_ctor_() { return new NAME; } \
static int NAME##_reg_ = reg_nc(#NAME, NAME##_ctor_);

#endif

body.cc

#ifndef _BODY_HPP_
#define _BODY_HPP_

#include "core.hpp"
#include <cmath>

#define _sq(x) ((x)*(x)) // square
#define _cb(x) abs((x)*(x)*(x)) // absolute value of cube
#define _cr(x) (unsigned char)(pow((x),1./3)) // cube root

struct colorful : public core {
  unsigned char RD(int i,int j){
    return (char)(_sq(cos(atan2(j-512.,i-512)/2))*255);
  }
  unsigned char GR(int i,int j){
    return (char)(_sq(cos(atan2(j-512.,i-512)/2-2*acos(-1.)/3))*255);
  }
  unsigned char BL(int i,int j){
    return (char)(_sq(cos(atan2(j-512.,i-512)/2+2*acos(-1.)/3))*255);
  }
};
REG(colorful);

struct flat : public core {
  unsigned char RD(int i,int j){
    float s=3./(j+99);
    return (int((i+DIM)*s+j*s)%2+int((DIM*2-i)*s+j*s)%2)*127;
  }
  unsigned char GR(int i,int j){
    float s=3./(j+99);
    return (int((i+DIM)*s+j*s)%2+int((DIM*2-i)*s+j*s)%2)*127;
  }
  unsigned char BL(int i,int j){
    float s=3./(j+99);
    return (int((i+DIM)*s+j*s)%2+int((DIM*2-i)*s+j*s)%2)*127;
  }
};
REG(flat);

struct ripple : public core {
  unsigned char RD(int i,int j){
    float s=3./(j+99);
    float y=(j+sin((i*i+_sq(j-700)*5)/100./DIM)*35)*s;
    return (int((i+DIM)*s+y)%2+int((DIM*2-i)*s+y)%2)*127;
  }
  unsigned char GR(int i,int j){
    float s=3./(j+99);
    float y=(j+sin((i*i+_sq(j-700)*5)/100./DIM)*35)*s;
    return (int((i+DIM)*s+y)%2+int((DIM*2-i)*s+y)%2)*127;
  }
  unsigned char BL(int i,int j){
    float s=3./(j+99);
    float y=(j+sin((i*i+_sq(j-700)*5)/100./DIM)*35)*s;
    return (int((i+DIM)*s+y)%2+int((DIM*2-i)*s+y)%2)*127;
  }
};
REG(ripple);

struct cripple : public core {
  unsigned char RD(int i,int j){
    float s=3./(j+99);
    float y=(j+sin((i*i+_sq(j-700)*5)/100./DIM)*35)*s;
    return (int((i+DIM)*s+y)%2+int((DIM*2-i)*s+y)%2)*127;
  }
  unsigned char GR(int i,int j){
    float s=3./(j+99);
    float y=(j+sin((i*i+_sq(j-700)*5)/100./DIM)*35)*s;
    return (int(5*((i+DIM)*s+y))%2+int(5*((DIM*2-i)*s+y))%2)*127;
  }
  unsigned char BL(int i,int j){
    float s=3./(j+99);
    float y=(j+sin((i*i+_sq(j-700)*5)/100./DIM)*35)*s;
    return (int(29*((i+DIM)*s+y))%2+int(29*((DIM*2-i)*s+y))%2)*127;
  }
};
REG(cripple);

struct randx: public core {
#define r(n)(rand()%n)
  unsigned char RD(int i,int j){
    static char c[1024][1024];return!c[i][j]?c[i][j]=!r(999)?r(256):RD((i+r(2))%1024,(j+r(2))%1024):c[i][j];
  }
  unsigned char GR(int i,int j){
    static char c[1024][1024];return!c[i][j]?c[i][j]=!r(999)?r(256):GR((i+r(2))%1024,(j+r(2))%1024):c[i][j];
  }
  unsigned char BL(int i,int j){
    static char c[1024][1024];return!c[i][j]?c[i][j]=!r(999)?r(256):BL((i+r(2))%1024,(j+r(2))%1024):c[i][j];
  }
};
REG(randx);

struct omg1: public core {
  unsigned char RD(int i,int j){
    float x=0,y=0;int k;for(k=0;k++<256;){float a=x*x-y*y+(i-768.0)/512;y=2*x*y+(j-512.0)/512;x=a;if(x*x+y*y>4)break;}return k>31?256:k*8;
  }
  unsigned char GR(int i,int j){
    float x=0,y=0;int k;for(k=0;k++<256;){float a=x*x-y*y+(i-768.0)/512;y=2*x*y+(j-512.0)/512;x=a;if(x*x+y*y>4)break;}return k>63?256:k*4;
  }
  unsigned char BL(int i,int j){
    float x=0,y=0;int k;for(k=0;k++<256;){float a=x*x-y*y+(i-768.0)/512;y=2*x*y+(j-512.0)/512;x=a;if(x*x+y*y>4)break;}return k;
  }
};
REG(omg1);

struct omg2: public core {
  unsigned char RD(int i,int j){
    double a=0,b=0,c,d,n=0;
    while((c=a*a)+(d=b*b)<4&&n++<880)
    {b=2*a*b+j*8e-9-.645411;a=c-d+i*8e-9+.356888;}
    return 255*pow((n-80)/800,3.);
  }
  unsigned char GR(int i,int j){
    double a=0,b=0,c,d,n=0;
    while((c=a*a)+(d=b*b)<4&&n++<880)
    {b=2*a*b+j*8e-9-.645411;a=c-d+i*8e-9+.356888;}
    return 255*pow((n-80)/800,.7);
  }
  unsigned char BL(int i,int j){
    double a=0,b=0,c,d,n=0;
    while((c=a*a)+(d=b*b)<4&&n++<880)
    {b=2*a*b+j*8e-9-.645411;a=c-d+i*8e-9+.356888;}
    return 255*pow((n-80)/800,.5);
  }
};
REG(omg2);

struct size800 : public core {
  size800() { DIM = 640; }

  unsigned char RD(int i,int j){
    return (unsigned short)sqrt((double)(_sq(i-DIM/2)*_sq(j-DIM/2))*2.0);
  }
  unsigned char GR(int i,int j){
    return (unsigned short)sqrt((double)((_sq(i-DIM/2)|_sq(j-DIM/2))*(_sq(i-DIM/2)&_sq(j-DIM/2))));
  }
  unsigned char BL(int i,int j){
    return (unsigned short)sqrt((double)(_sq(i-DIM/2)&_sq(j-DIM/2))*2.0);
  }
};
REG(size800);

#endif

main.cc

/*
http://www.zhihu.com/question/30262900
http://www.matrix67.com/blog/archives/6039
http://codegolf.stackexchange.com/questions/35569/tweetable-mathematical-art
*/

#include "body.hpp"
#include "pic.hpp"

int main(int argc, char * argv[]) {
  if(argc < 2) {
    printf("Usage: %s <-all|name> [bmp|ppm]\n", argv[0]);
    printf("name list:\n");
    for(unsigned int i = 0; i < nc_num; i++) {
      printf("\t%s\n", nc[i].name);
    }
    return 1;
  }

  for(unsigned int i = 0; i < nc_num; i++) {
    if(strcmp(argv[1], nc[i].name)==0 || strcmp(argv[1], "-all")==0) {
      core * c = nc[i].ctor();
      pic * p;

      if(argc > 2 && strcmp(argv[2],"ppm") == 0) {
        p = new ppm(c->DIM, c->DIM);
      } else {
        p = new bmp(c->DIM, c->DIM);
      }

      for(unsigned int y = 0; y < c->DIM; y++)
        for(unsigned int x = 0; x < c->DIM; x++)
          p->set(x, y, c->RD(x,y), c->GR(x,y), c->BL(x,y));

      p->save(nc[i].name);

      delete p;
      delete c;
    }
  }

  return 0;
}

Makefile

ifndef T
A=-all
endif

run :: math_art
	./math_art $A

math_art : pic.hpp core.hpp body.hpp main.cc
	g++ -Wall main.cc -lm -o math_art

clean ::
	rm -f math_art *.bmp *.ppm

build.bat

@echo off
call "C:\Program Files\Microsoft Visual Studio 10.0\VC\vcvarsall.bat" x86
del *.bmp 2>nul
cl main.cc & main.exe -all & del main.obj main.exe
pause


214049_EWtc_947735.bmp

214051_09vU_947735.bmp


转载于:https://my.oschina.net/2bit/blog/498193

weixin_33786077
关注
  • 0
    点赞
  • 踩
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
StackExchange 上的 Tweetable Mathematical Art 比赛优秀解法合集
01-10
StackExchange 上的 Tweetable Mathematical Art 比赛优秀解法合集，用C/C++生成1024x1024图像，图像格式ppm，压缩包中有转JPG后的图像。
Tweetable Mathematical Art
MCM的专栏
11-02 1028
http://weibo.com/p/1001603759188458255826
Tweetable Mathematical Art 太极阴阳图
xshbx的专栏
04-24 3696
从Matrix67 blog里面看到了Tweetable Mathematical Art。 使用C++实现3个函数RD(i, j), GR(i, j), BL(i, j)分别返回像素(i, j)的R/G/B值，最终生成1幅1024*1024的图片。要求是每个函数代码不超过140字符，即Tweetable。 来个太极八卦图，娱乐一下。
Tweetable Mathematical Art 习作
weixin_30701521的博客
09-03 113
从Matrix67 blog里面看到了Tweetable Mathematical Art比赛，感觉挺有意思。 这个比赛的方法很简单，实现3个函数RD(i, j), GR(i, j), BL(i, j)分别返回像素(i, j)的R/G/B值，最终生成1幅1024*1024的图片，看谁的图片被赞最多。要求是每个函数代码不超过140字符，即Tweetable。 手痒写了一个：放大镜效果，...
发现个有趣的东西：Tweetable Mathematical Art（用三个140字符以内的函数生成一个1024尺寸的图片）
主要记录了工作之余胡乱学习的知识，大部分是入门级别的内容。 主要是为了自己未来有用到时能快速上手，相当于是个备忘录。 当然，也很开心其他人可以从中找到对自己有用的信息。^_^
11-16 809
发现 我是在看《构建之法》这本书时，看到作者提到这个： 好厉害！用三段140字符以内的代码生成一张1024×1024的图片_IT新闻_博客园 这是2014年一个人在 Code Golf Stack Exchange (a question and answer site for programming puzzle enthusiasts and code golfers) 发起的编程挑战：math - Tweetable Mathematical Art - Code Golf Stack Exchang
Tweetable Mathematical Art 算法制图程序
11-07
"Tweetable Mathematical Art 算法制图程序"是一种利用编程技术创作艺术图形的方法，它将数学原理与美学结合，通过算法生成独特的视觉效果。在这个项目中，用户可以制作28种不同的数学艺术图形，为编程爱好者和数学...
tweetable-polyglot-png:将最多3MB的数据打包到可发推的PNG多语言文件中
03-17
将最多3MB的数据打包到可发布Twitter消息的PNG多语言文件中。如何？ Twitter会从PNG上载中剥离不必要的数据，但是，如果整个图像文件符合避免重新编码的要求，它们就不会从IDAT块内的DEFLATE流中剥离尾随数据。...
Spring.Boot.in.Action.2015.12.pdf
09-14
ary 9Oish-character tweetable Boot application to an in-depth analysis of Boots Actuator in chapter 7, which enables a host of auto-magical operational features required for any production application...
Kyle McCormick 在 StackExchange 上发起了一个叫做 Tweetable Mathematical Art 的比赛，参赛者需要用 C++ 语言编写 RD 、 GR 、 BL ...
weixin_35756892的博客
12-31 144
这题目要求用 Kotlin 语言实现一个 Android 应用，使用 ImageView 显示由三个函数生成的图像。这三个函数分别为 RD、GR 和 BL，每个函数都接收两个整型参数 i 和 j，要求返回一个 0 到 255 之间的整数，表示位于 (i,j) 的像素点的颜色值。举个例子，如果 RD(0, 0) 和 GR(0, 0) 返回的都是 0，但 BL(0, 0) 返回的是 255，那么图像的...
数学之美，C++之美，代码之美
OliverkingLi的博客
09-25 691
用三段 140 字符以内的代码生成一张 1024×1024 的图片 Kyle McCormick 在 StackExchange 上发起了一个叫做 Tweetable Mathematical Art 的比赛，参赛者需要用三条推这么长的代码来生成一张图片。具体地说，参赛者需要用 C++ 语言编写 RD 、 GR 、 BL 三个函数，每个函数都不能超过 140 个字符。每个函数都会接到 i 和
【cool】三段 140 字以内代码生成一张图
nickwong & nickwang
09-12 1105
Tweetable Mathematical Art比赛，用三段 140 字以内代码生成一张1024*1024图，感觉挺cool。 自己模(chao)仿(xi)了几张，感觉后三张不忍直视。 5张1024*1024图片的代码如下： /*picture1 unsigned char RD(int i,int j){ return i&&j?((i%j)|(
I do. 【附：C++用不多于一条推博字符数的代码能作出怎样的图像？】
stereohomology
05-28 1300
以上为原创主贴 以下为转帖及少量翻译。 很多C++程序员在进行的一场比赛： 用不多于一条微博的字符数的代码能作出怎样的数学图像？ http://codegolf.stackexchange.com/questions/35569/tweetable-mathematical-art 悄悄问一声：看到有人画美国国旗，但那个难度似乎高些，画的不够完整； 我国国旗简单明快，谁能把
趣味C++
火雨(Nick)
09-06 2044
用三段 140 字符以内的代码生成一张 1024×1024 的图片                 Kyle McCormick 在 StackExchange 上发起了一个叫做TweetableMathematical Art 的比赛，参赛者需要用三条推这么长的代码来生成一张图片。具体地说，参赛者需要用 C++ 语言编写 RD 、 GR 、 BL 三个函数，每个函数都不能超过 140
百行内代码生成图片
musks的博客
05-25 828
看到一个流弊的精简的代码生成各种图片 好像是之前一个比赛吧，先上官网 https://codegolf.stackexchange.com/questions/35569/tweetable-mathematical-art 上面已经有生成不少漂亮的图片，先上两张拔过来的图片 我还是看看代码吧，尽管那三行不同函数中所包含的数学公式看不懂。
用三段 140 字符以内的代码生成一张 1024×1024 的图片<转载>
weixin_34080903的博客
08-19 184
原文出处：http://www.matrix67.com/blog/archives/6039#more-6039 http://codegolf.stackexchange.com/questions/35569/tweetable-mathematical-art Kyle McCormick 在 StackExchange 上发起了一个叫做 Tweetable Mathematical...
基础作品制作思路.mp4
最新发布
07-09
。
完整的修正和优化代码.docx
07-09
完整的修正和优化代码.docx
用Kotlin语言实现Android版本的Tweetable Mathematical Art，用ImageView显示不同算法生成的图像。这里面难点在哪？
08-04
实现 Android 版本的 Tweetable Mathematical Art 可能的难点有以下几点： 1. 使用 Kotlin 语言编写 Android 应用程序的代码。如果您不熟悉 Kotlin，那么可能需要花费一些时间来学习基础知识。 2. 使用 ImageView ...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值