火车运煤问题（驴运萝卜问题）

#include <stdio.h>
#include <stdlib.h>
int main(){
    FILE * out;
    int total = 3000;

    int a[] = {1, 1, 1};  // 记录每段 长度， 总长度为 1000 km
    int t[] = {0, 0, 0};   // 记录每段路  返回经过的次数
    int x = 0, xb = 0;
    int xy = 0, xyb = 0;  // x + y
    int remain, remainb;
    int i, j;
    int max = 0;
    // 创建一个文件，用来记录结果
    if((out=fopen("res.txt", "a++"))==NULL){
        fprintf(stdout, "Cannot open \"res.txt\" file\n");
        exit(1);
    }
    // invariant: a[0] + a[1] + a[2] = 1000
    // invariant: 剩下的煤量超过剩下的距离
    for(i = 1; i <500; i++){
        a[0] = i;    // 第一段路程  0001
        t[0] = 2;
        xb = a[0] * t[0];
        xyb = a[0] * ( 2 * t[0] + 1);
        remainb = total - xyb;  // 记录 x, xy, remain 以便于在内层循环起始处恢复
        for(j = 1; j < 500; j++){
            x = xb, xy = xyb, remain = remainb;  // 每次都重新对 x, xy, remain赋值
            a[1] = j;           // 第二段路程     0002
            t[1] = remain % 1000 == 0? (remain/1000 - 1):(remain/1000);
            x += a[1] * t[1];
            xy += a[1] * ( 2 * t[1] + 1);
            if(x>=1000) continue;   // x + y < 3000, 所以2 * x < 2000, 即 x < 1000
            if(xy>=3000) continue;  // x + y < 3000
            remain = total - xy;

            a[2] = 1000 - i - j;   // 第三段路程    0003
            if(a[2] > remain) continue;
            t[2] = remain % 1000 == 0? (remain/1000 - 1):(remain/1000);
            x += a[2] * t[2];
            xy += a[2] * ( 2 * t[2] + 1);
            if(x>=1000) continue;  // x + y < 3000, 所以2 * x < 2000, 即 x < 1000
            if(xy>=3000) continue;  // x + y < 3000

            remain = total - xy;
            // 将 结果打印到 文件中，便于查看
            fprintf(out, "i = %d, j = %d, x = %d, xy = %d, r = %d\n", i,  j, x, xy, r);  
            if(remain > max) max = remain;
        }  // end of innner for loop j
    }  // end of outer for (i
    printf("max = %d\n", max);
    if(fclose(out)!=0){
        fprintf(stdout, "close failed");
    }
    getchar();
    return 0;