uva 10827 - Maximum sum on a torus

Problem H
Maximum sum on a torus
Input: Standard Input

Output: Standard Output

A grid that wraps both horizontally and vertically is called a torus. Given a torus where each cell contains an integer, determine the sub-rectangle with the largest sum. The sum of a sub-rectangle is the sum of all the elements in that rectangle. The grid below shows a torus where the maximum sub-rectangle has been shaded.

1

-1

0

0

-4

2

3

-2

-3

2

4

1

-1

5

0

3

-2

1

-3

2

-3

2

4

1

-4

Input

The first line in the input contains the number of test cases (at most 18). Each case starts with an integer N (1≤N≤75) specifying the size of the torus (always square). Then follows N lines describing the torus, each line containing N integers between -100 and 100, inclusive.

Output

For each test case, output a line containing a single integer: the maximum sum of a sub-rectangle within the torus.

Sample Input Output for Sample Input

2
5
1 -1 0 0 -4
2 3 -2 -3 2
4 1 -1 5 0
3 -2 1 -3 2
-3 2 4 1 -4
3
1 2 3
4 5 6
7 8 9
15

45


Problem setter: Jimmy Mårdell

Special Thanks: Derek Kisman, Md. Kamruzzaman

uva108加强版,矩形是循环的,补上三个,直觉让我们想到先枚举所有可能出现的矩形,然后套上108的代码,果断TLD,稍微优化了下,还是枚举所有行的组合,再枚举计算最大子序列的起始点,卡着时间过去了
#include<stdio.h>
#include<string.h>
#define inf -9999999
int b[200],n;
int count()
{
 int i,j,f=1,max=inf,sum;
 for (i=1;i<=n;i++)
 {
  if  (b[i]>max) max=b[i];
  if (b[i]>=0)  f=0;
 }
 if (f) return max;
 for (j=1;j<=n;j++)
 {
  sum=0;
  for (i=0;i<n;i++)
  {
   sum=sum+b[i+j];
   if (sum>max) max=sum;
   if (sum<0) sum=0;
  }
 }
 return max;
}
int main()
{
 int T,i,j,p,q,k,t,a[200][200],max,ans,m;
 scanf("%d\n",&T);
 while (T--)
 {
  scanf("%d",&n);
  for (i=1;i<=n;i++)
  for (j=1;j<=n;j++)
  {
   scanf("%d",&a[i][j]);
   a[i+n][j]=a[i][j];
   a[i][j+n]=a[i][j];
   a[i+n][j+n]=a[i][j];
  }
  max=inf;
  for (k=0;k<n;k++)
  {
   for (i=1;i<=2*n;i++)
   {
    if (i+k<=2*n)
    {
      memset(b,0,sizeof(b));
      for (j=i;j<=i+k;j++)
      for (t=1;t<=2*n;t++)
      b[t]+=a[j][t];
      ans=count();
      if (ans>max) max=ans;
    }
   }
  }
  printf("%d\n",max);
 }
 return 0;
}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
extern UFUNEXPORT int UF_MODL_ask_face_data( tag_t face ,/* <I> Face obj_id / int * type ,/ <O> Face type is NX surface type code 16 = cylinder 17 = cone 18 = sphere 19 = revolved (toroidal) 20 = extruded 22 = bounded plane 23 = fillet (blend) 43 = b-surface 65 = offset surface 66 = foreign surface / double point[] ,/ <O,len:3> Point information is returned according to the value of type as follows. Plane = Position in plane Cylinder= Position on axis Cone = Position on axis Sphere = Center position Torus = Center position Revolved = Position on axis / double dir[] ,/ <O,len:3> Direction information is returned according to the value of type as follows. Plane = Normal direction Cylinder= Axis direction Cone = Axis direction Torus = Axis direction Revolved = Axis direction / double box[] ,/ <O,len:6> Face boundary. The coordinates of the opposite corners of a rectangular box with sides parallel to X, Y, and Z axes (Absolute Coordinate System) are returned. The box contains the specified face and is usually close to the minimum possible size, but this is not guaranteed. box[0] = Xmin box[1] = Ymin box[2] = Zmin box[3] = Xmax box[4] = Ymax box[5] = Zmax / double * radius ,/ <O> Face major radius: For a cone, the radius is taken at the point[3] position on the axis. For a torus, the radius is taken at the major axis. / double * rad_data ,/ <O> Face minor radius: only a torus or cone has rad_data as a minor radius. For a cone, rad_data is the half angle in radians. For a torus, rad_data is taken at the minor axis. / int * norm_dir / <O> Face normal direction: +1 if the face normal is in the same direction as the surface normal (cross product of the U- and V-derivative vectors), -1 if reversed. */ ); 详解
06-01
`UF_MODL_ask_face_data` 函数是 UG/Open C API 中的一个函数,用于查询面(face)的数据。以下是该函数的详细解释: 函数原型: ```c extern UFUNEXPORT int UF_MODL_ask_face_data( tag_t face, /* <I> 面对象的 ID */ int *type, /* <O> 面的类型,NX 表面类型代码 16 = 圆柱体,17 = 圆锥体,18 = 球体,19 = 旋转体(环形),20 = 挤出体,22 = 平面,23 = 圆弧,43 = B 曲面,65 = 偏移曲面,66 = 外部曲面 */ double point[], /* <O,len:3> 根据类型返回点信息,具体如下:平面=平面上的位置,圆柱体=轴上的位置,圆锥体=轴上的位置,球体=中心位置,旋转体=轴上的位置 */ double dir[], /* <O,len:3> 根据类型返回方向信息,具体如下:平面=法线方向,圆柱体=轴向,圆锥体=轴向,旋转体=轴向 */ double box[], /* <O,len:6> 面边界,返回一个矩形盒子的对角线上的两个点的坐标(参考绝对坐标系),盒子包含了指定的面并且通常是最小可能的,但是不能保证。 */ double *radius, /* <O> 面的主要半径:对于圆锥体,半径是在轴上的点[3]处取的。对于环体,半径是在主轴上取的。 */ double *rad_data, /* <O> 面的次要半径:只有环体或圆锥体具有 rad_data 作为次要半径。对于圆锥体,rad_data 是半角的一半,对于环体,rad_data 是在副轴上取的。 */ int *norm_dir /* <O> 面的法线方向:+1 表示面法线与曲面法线(U- 和 V- 导数向量的叉积)相同,-1 表示相反。 */ ); ``` 参数说明: - `face`:要查询的面对象的 ID。 - `type`:返回该面的类型,NX 表面类型代码。 - `point[]`:根据类型返回点信息。具体如下: - 平面:平面上的位置。 - 圆柱体:轴上的位置。 - 圆锥体:轴上的位置。 - 球体:中心位置。 - 旋转体:轴上的位置。 - `dir[]`:根据类型返回方向信息。具体如下: - 平面:法线方向。 - 圆柱体:轴向。 - 圆锥体:轴向。 - 旋转体:轴向。 - `box[]`:面边界,返回一个矩形盒子的对角线上的两个点的坐标(参考绝对坐标系),盒子包含了指定的面并且通常是最小可能的,但是不能保证。 - `radius`:返回面的主要半径。对于圆锥体,半径是在轴上的点[3]处取的。对于环体,半径是在主轴上取的。 - `rad_data`:返回面的次要半径。只有环体或圆锥体具有 rad_data 作为次要半径。对于圆锥体,rad_data 是半角的一半,对于环体,rad_data 是在副轴上取的。 - `norm_dir`:返回面的法线方向。+1 表示面法线与曲面法线(U- 和 V- 导数向量的叉积)相同,-1 表示相反。 函数返回值: - `0` 表示成功。 - 非 `0` 值表示失败。 使用示例: ```c #include <uf_modl.h> #include <stdio.h> int main() { tag_t face; int type; double point[3], dir[3], box[6], radius, rad_data; int norm_dir; // 获取一个面的 ID // ... int ret = UF_MODL_ask_face_data(face, &type, point, dir, box, &radius, &rad_data, &norm_dir); if (ret != 0) { printf("查询面数据失败!\n"); return 1; } printf("面类型:%d\n", type); printf("点坐标:(%lf, %lf, %lf)\n", point[0], point[1], point[2]); printf("方向向量:(%lf, %lf, %lf)\n", dir[0], dir[1], dir[2]); printf("盒子坐标:(%lf, %lf, %lf), (%lf, %lf, %lf)\n", box[0], box[1], box[2], box[3], box[4], box[5]); printf("主要半径:%lf\n", radius); printf("次要半径:%lf\n", rad_data); printf("法线方向:%d\n", norm_dir); return 0; } ``` 该示例中,我们首先获取了一个面对象的 ID `face`,然后调用 `UF_MODL_ask_face_data` 函数查询该面的数据。函数返回成功后,我们打印出面的类型、点坐标、方向向量、盒子坐标、主要半径、次要半径和法线方向。

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值