读取in.bmp文件，判断其中哪些点是连通的

【问题描述】读取in.bmp文件，判断其中哪些点是连通的，并求各个连通分量的面积。
BMP是常见的图像存储格式。
如果用来存黑白图像（颜色深度=1），则其信息比较容易读取。

与之相关的数据：
（以下偏移均是从文件头开始）
偏移：10字节, 长度4字节： 图像数据真正开始的位置。
偏移：18字节, 长度4字节： 位图的宽度，单位是像素。
偏移：22字节, 长度4字节： 位图的高度，单位是像素。

从图像数据开始处，每个像素用1个二进制位表示。从图片的底行开始，一行一行向上存储。
Windows规定图像文件中一个扫描行所占的字节数必须是4字节的倍数，不足的位均以 0 填充。例如，图片宽度为45像素，实际上每行会占用8个字节。
可以通过Windows自带的画图工具生成和编辑二进制图像。 需要在“属性”中选择“黑白”，指定为二值图像。 可能需要通过 查看 | 缩放 | 自定义... 把图像变大比例一些，更易于操作。
图像的左下角为图像数据的开始位置。白色对应1，黑色对应0
我们可以定义：两个点距离如果小于2个像素，则认为这两个点连通。 也就是说：以一个点为中心的九宫格中，围绕它的8个点与它都是连通的。
程序的目标是：根据给定的黑白位图，分析出所有独立的连通分量，输出每个连通分量的面积。

所谓面积，就是一个连通部分所含黑色像素点的个数。

输入数据固定存在in.bmp中，先输出in.bmp图片中点的信息（白点用空格表示，黑点用*表示）；

再输出in.bmp图片中各个连通分量的面积（输出的连通分量面积，按由小到大的顺序）。

【输入样例】

【输出样例】

 

#include <iostream>
#include <bitset>
#include <vector>
#include <algorithm>
using namespace std;
const int N = 1024;
char vis[N][N] = {0};
char g[N][N] = {0};
int di[] = {-1, 1, 0, 0, -1, -1, 1, 1};
int dj[] = {0, 0, -1, 1, -1, 1, -1, 1};
int width, height, cnt; // cnt用来表示当前连通分量的大小
void dfs(int i, int j) {
	if (i < 0 || i >= height || j < 0 || j >= width || vis[i][j]) { // 超出图边界的情况或已经被访问过
		return;
	}
	vis[i][j] = 1; // 标记这个点已经访问过了
	if (g[i][j] != '*') { // 如果不是黑点，直接退出
		return;
	}
	cnt++;
	for (int dir = 0; dir < 8; dir++) { // 对周围8个点做dfs
		dfs(i + di[dir], j + dj[dir]);
	}
}
int getDataFrom(FILE* fp, int offset) {
	int data;
	fseek(fp, offset, 0);
	fread(&data, sizeof(int), 1, fp);
	return data;
}
void readBMP(FILE* fp) {
	int start = getDataFrom(fp, 10);
	fseek(fp, start, 0);
	int fill = (32 - width & 31) & 31;
	printf("%d+%d=%d\n", width, fill, width + fill);
	for (int i = height - 1; i >= 0; i--) {
		for (int j = 0; j < width + fill; ) {
			char c;
			fread(&c, 1, 1, fp); // 读1个字节
			bitset<8> b(c);
			for (int bit = b.size() - 1; bit >= 0; bit--, j++) {
				if (b[bit]) { // 是白色
					g[i][j] = ' ';
				} else { // 是黑色
					g[i][j] = '*';
				}
			}
		}
	}
}
void printBMP() {
	for (int i = 0; i < height; i++) {
		for (int j = 0; j < width; j++) {
			printf("%c", g[i][j]);
		}
		printf("\n");
	}
}
void calConnectedPart() {
	vector<int>ans;
	// int a[128] = {0}, k = 0;
	for (int i = 0; i < height; i++) { // 从每个顶点开始做一次dfs，一次dfs只能遍历他所在的连通分量
		for (int j = 0; j < width; j++) {
			cnt = 0;
			dfs(i, j);
			if (cnt != 0) {
				ans.push_back(cnt);
				// a[k++] = cnt;
			}
		}
	}
	sort(ans.begin(), ans.end()); // c++实现的快排
	printf("The area of each connected part is:%d", ans[0]);
	for (int i = 1; i < ans.size(); i++) {
		printf(",%d", ans[i]);
	}
}
int main() {
	FILE* fp = fopen("in.bmp", "rb"); // 打开文件
	width = getDataFrom(fp, 18); // 读取长宽数据
	height = getDataFrom(fp, 22);
	readBMP(fp);
	fclose(fp);
	printBMP();
	calConnectedPart();
	return 0;
}