回 帖 发 新 帖 刷新版面

主题:[讨论]NOIP2006 第一题(能量项链)!

NO.1能量项链
【问题描述】
在Mars星球上,每个Mars人都随身佩带着一串能量项链。在项链上有N颗能量珠。能量珠是一颗有头标记与尾标记的珠子,这些标记对应着某个正整数。并且,对于相邻的两颗珠子,前一颗珠子的尾标记一定等于后一颗珠子的头标记。因为只有这样,通过吸盘(吸盘是Mars人吸收能量的一种器官)的作用,这两颗珠子才能聚合成一颗珠子,同时释放出可以被吸盘吸收的能量。如果前一颗能量珠的头标记为m,尾标记为r,后一颗能量珠的头标记为r,尾标记为n,则聚合后释放的能量为 (Mars单位),新产生的珠子的头标记为m,尾标记为n。
需要时,Mars人就用吸盘夹住相邻的两颗珠子,通过聚合得到能量,直到项链上只剩下一颗珠子为止。显然,不同的聚合顺序得到的总能量是不同的,请你设计一个聚合顺序,使一串项链释放出的总能量最大。
例如:设N=4,4颗珠子的头标记与尾标记依次为(2,3) (3,5) (5,10) (10,2)。我们用记号⊕表示两颗珠子的聚合操作,(j⊕k)表示第j,k两颗珠子聚合后所释放的能量。则第4、1两颗珠子聚合后释放的能量为:
(4⊕1)=10*2*3=60。
这一串项链可以得到最优值的一个聚合顺序所释放的总能量为
((4⊕1)⊕2)⊕3)=10*2*3+10*3*5+10*5*10=710。

【输入文件】
输入文件的第一行是一个正整数N(4≤N≤100),表示项链上珠子的个数。第二行是N个用空格隔开的正整数,所有的数均不超过1000。第i个数为第i颗珠子的头标记(1≤i≤N),当i<N时,第i颗珠子的尾标记应该等于第i+1颗珠子的头标记。第N颗珠子的尾标记应该等于第1颗珠子的头标记。
至于珠子的顺序,你可以这样确定:将项链放到桌面上,不要出现交叉,随意指定第一颗珠子,然后按顺时针方向确定其他珠子的顺序。
【输出文件】
输出文件只有一行,是一个正整数E(E≤2.1*109),为一个最优聚合顺序所释放的总能量。
【输入样例】
4
2 3 5 10
【输出样例】
710

回复列表 (共4个回复)

沙发

动态规划算法如下:
第一种方法:
#include<stdio.h>
#include<stdlib.h>
#include<memory.h>
#include<time.h>
#define MAX 200

int main()
{   
    clock_t begin,end;
    FILE *in,*out;
    long m[MAX][MAX],mar[MAX];
    long E,i,j,k,l,n,max;
    memset(m,0,sizeof(m));
    begin=clock();
    in=fopen("D:\\energy10.in","r");
    fscanf(in,"%ld",&n);
     
    for(i=0;i<n;i++)
    fscanf(in,"%ld",&mar[i]);
    fclose(in);
    
    //scanf("%ld",&mar[i]);
    for(k=1;k<n;k++)
     for(i=0;i<n;i++)
     {
       j=i+k;
       max=0;
       for(l=i;l<j;l++)
       {
         max=m[i][l%n]+m[(l+1)%n][j%n]+mar[i]*mar[(l+1)%n]*mar[(j+1)%n];
         if(max>m[i][j%n])
         m[i][j%n]=max;
       }    
     }
    E=m[0][n-1];
    for(i=1,j=0;i<n;i++,j++)
    if(m[i][j]>E)
    E=m[i][j];
    out=fopen("D:\\energy10.result","w");
    fprintf(out,"%ld",E);
    fclose(out);
    end=clock();
    printf("run %f seconds!",(double)(end-begin)/CLOCKS_PER_SEC);
    return 0;
}

    

板凳

第二种算法:
#include<stdio.h>
#include<memory.h>
#include<stdlib.h>
#include<time.h>
#define MAX 200

int main()
{
  FILE *in,*out;
  clock_t begin,end;
  long m[2*MAX][2*MAX],mar[2*MAX];
  long i,j,k,h,n,max,E=0;
   
  begin=clock();
  in=fopen("D:\\energy10.in","r");
  fscanf(in,"%ld",&n);
  for(i=1;i<=n;i++)
  fscanf(in,"%ld",&mar[i]);
  fclose(in);
  for(i=n+1,j=1;i<=2*n;i++,j++)
  mar[i]=mar[j];
  for(k=2;k<=n;k++)
    for(i=1;i<=2*n-k;i++)
    {
      j=i+k-1;
      m[i][j]=m[i+1][j]+mar[i]*mar[i+1]*mar[j+1];
      max=0;
      for(h=i+1;h<j;h++)
       {
        max=m[i][h]+m[h+1][j]+mar[i]*mar[h+1]*mar[j+1];
        if(max>m[i][j])
        m[i][j]=max;
       }
      }
      
   out=fopen("D:\\energy10.out","w");
   for(i=1,j=n;i<=n;i++,j++)
   if(m[i][j]>E)
   E=m[i][j];
   fprintf(out,"%ld",E);
   fclose(out);
   end=clock();
   printf("\nRUN %f  seconds !",(double)(end-begin)/CLOCKS_PER_SEC);
   return 0;
}

3 楼

最麻烦的是破环  可以用X:=(I+K-1) MOD N+1来解决 其余简单

4 楼

不用破也可以

我来回复

您尚未登录,请登录后再回复。点此登录或注册