【区间DP】矩阵乘法

矩阵乘法

问题描述
    一个A x B的矩阵乘以一个B x C的矩阵将得到一个A x C的矩阵，时间复杂度为A x B x C。矩阵乘法满足结合律（但不满足交换律）。顺序给出n个矩阵的大小，请问计算出它们的乘积的最少需要花费多少时间。

输入数据
    第一行输入一个正整数n，表示有n个矩阵。
    接下来m行每行两个正整数Xi,Yi，其中第i行的两个数表示第i个矩阵的规模为Xi x Yi。所有的Xi、Yi<=100。输入数据保证这些矩阵可以相乘。

输出数据
    输出最少需要花费的时间。

样例输入
3
10 100
100 5
5 50

样例输出
7500

样例说明
    顺序计算总耗时7500；先算后两个总耗时75000。

时间限制
    各测试点1秒

内存限制
    你的程序将被分配32MB的运行空间

数据范围
    n<=100。

============

这道题其实是一道区间动规

所需要的是

1.先枚举区间长度

2.再枚举起点

3.最后枚举中点

-----------------

所遇到的错误（PASCAL语言）

运行错误216：存取非法主要是因为在定义范围外.

============

var x,y:array[1..100]of longint; f:array[1..100,1..100]of longint; n:longint;procedure init;begin assign(input,'matrix.in'); assign(output,'matrix.out'); reset(input); rewrite(output);end;procedure terminate;begin close(input); close(output); halt;end;function min(a,b:longint):longint;begin if a>b then exit(b); exit(a);end;procedure main;var i,j,k:longint;begin readln(n); fillchar(f,sizeof(f),$7); for i:=1 to n do begin readln(x[i],y[i]); f[i,i]:=0; end; for k:=2 to n do //枚举区间长度 for i:=1 to n-k+1 do //枚举起点. for j:=i to i+k-2 do //枚举中点. begin f[i,i+k-1]:=min(f[i,i+k-1],f[i,j]+f[j+1,i+k-1]+x[i]*x[j+1]*y[i+k-1]); end; writeln(f[1,n]);end;begin init; main; terminate;end.