用c语言写一个程序，要求数据管理：程序中用来表示三维物体的数据必须存放在磁盘文件中（示例文件见Bunny.off），三维物体由多个三角形组成，文件中包含了三角形各个顶点的坐标位置以及构成三角形的点的信息。可以考虑定义一个三角形类，对类中元素进行操作。

好的，下面是一个简单的示例程序，用于读取三角形数据文件，并输出每个三角形的面积和体积。注意，这个程序假定三角形的顶点坐标已经按照逆时针方向排列好了。

#include <stdio.h>
#include <stdlib.h>
#include <math.h>

typedef struct {
  float x, y, z;
} Vector;

typedef struct {
  int p1, p2, p3;
} Triangle;

float dot(Vector v1, Vector v2) {
  return v1.x * v2.x + v1.y * v2.y + v1.z * v2.z;
}

Vector cross(Vector v1, Vector v2) {
  Vector result;
  result.x = v1.y * v2.z - v1.z * v2.y;
  result.y = v1.z * v2.x - v1.x * v2.z;
  result.z = v1.x * v2.y - v1.y * v2.x;
  return result;
}

float length(Vector v) {
  return sqrtf(v.x * v.x + v.y * v.y + v.z * v.z);
}

float area(Vector v1, Vector v2, Vector v3) {
  Vector cross_product = cross((Vector) {v2.x - v1.x, v2.y - v1.y, v2.z - v1.z},
                                (Vector) {v3.x - v1.x, v3.y - v1.y, v3.z - v1.z});
  return 0.5f * length(cross_product);
}

float volume(Vector v1, Vector v2, Vector v3, Vector v4) {
  Vector cross_product = cross((Vector) {v2.x - v1.x, v2.y - v1.y, v2.z - v1.z},
                                (Vector) {v3.x - v1.x, v3.y - v1.y, v3.z - v1.z});
  return dot(cross_product, (Vector) {v4.x - v1.x, v4.y - v1.y, v4.z - v1.z}) / 6.0f;
}

int main(int argc, char **argv) {
  if (argc < 2) {
    printf("Usage: %s <filename>\n", argv[0]);
    return 1;
  }
  
  FILE *file = fopen(argv[1], "r");
  if (file == NULL) {
    printf("Failed to open file %s\n", argv[1]);
    return 1;
  }
  
  int num_vertices, num_triangles, num_edges;
  fscanf(file, "OFF\n%d %d %d\n", &num_vertices, &num_triangles, &num_edges);
  
  Vector *vertices = (Vector*)malloc(num_vertices * sizeof(Vector));
  Triangle *triangles = (Triangle*)malloc(num_triangles * sizeof(Triangle));
  
  for (int i = 0; i < num_vertices; i++) {
    fscanf(file, "%f %f %f\n", &vertices[i].x, &vertices[i].y, &vertices[i].z);
  }
  
  for (int i = 0; i < num_triangles; i++) {
    int num_vertices;
    fscanf(file, "%d %d %d %d\n", &num_vertices, &triangles[i].p1, &triangles[i].p2, &triangles[i].p3);
    if (num_vertices != 3) {
      printf("Error: unsupported polygon type\n");
      return 1;
    }
    
    Vector v1 = vertices[triangles[i].p1];
    Vector v2 = vertices[triangles[i].p2];
    Vector v3 = vertices[triangles[i].p3];
    
    float area_triangle = area(v1, v2, v3);
    printf("Triangle %d: area = %f\n", i, area_triangle);
    
    if (i < num_triangles - 1) {
      for (int j = i + 1; j < num_triangles; j++) {
        int p1 = triangles[i].p1, p2 = triangles[i].p2, p3 = triangles[i].p3;
        int q1 = triangles[j].p1, q2 = triangles[j].p2, q3 = triangles[j].p3;
        
        if (p1 == q1 || p1 == q2 || p1 == q3 ||
            p2 == q1 || p2 == q2 || p2 == q3 ||
            p3 == q1 || p3 == q2 || p3 == q3) {
          Vector v4 = vertices[p1 == q1 || p1 == q2 || p1 == q3 ? p1 :
                                p2 == q1 || p2 == q2 || p2 == q3 ? p2 :
                                p3 == q1 || p3 == q2 || p3 == q3 ? p3 : -1];
          
          float volume_tetrahedron = volume(v1, v2, v3, v4);
          printf("Tetrahedron %d-%d: volume = %f\n", i, j, volume_tetrahedron);
        }
      }
    }
  }
  
  fclose(file);
  free(vertices);
  free(triangles);
  
  return 0;
}

这个程序假定输入文件的格式是二进制 OFF 格式，可以按照以下命令进行编译和运行：

gcc -o program program.c
./program Bunny.off

在这个程序中，我们定义了一个 `Vector` 结构体来表示三维向量，以及一个 `Triangle` 结构体来表示三角形。然后，我们读取输入文件，将顶点和三角形数据存储在内存中，并且对于每个三角形，计算其面积，并且检查它是否和其它三角形共享一个或多个顶点。如果是，我们就可以计算由这两个三角形和一个共享点构成的四面体的体积。