?»??»?#include<stdio.h>
#include<math.h>
double a0,h1,a1,a2,c1,c2;
double b1,b2;
double x1,x0,x2,x3,x4;
double w1,w2,w3,w4;
double a,b,c,d,e,f,g,h,i,j,k,l,m,n,o,p,x,E;
double f1,f2;
double fx,k1=0.618;

int main(int argc, char* argv[])
{
 printf("fx=p*(a*x1+b*x2)^c+d*(e*x3-f*x4)^g+h*(i*x2-j*x3)^k+l*(m*x1+n*x4)^o:\n");
 printf("input a-o\n");
 printf("a=");
 scanf("%lg",&a);
 printf("b=");
 scanf("%lg",&b);
 printf("c=");
 scanf("%lg",&c);
 printf("d=");
 scanf("%lg",&d);
 printf("e=");
 scanf("%lg",&e);
 printf("f=");
 scanf("%lg",&f);
 printf("g=");
 scanf("%lg",&g);
 printf("h=");
 scanf("%lg",&h);
 printf("i=");
 scanf("%lg",&i);
 printf("j=");
 scanf("%lg",&j);
 printf("k=");
 scanf("%lg",&k);
 printf("l=");
 scanf("%lg",&l);
 printf("m=");
 scanf("%lg",&m);
 printf("n=");
 scanf("%lg",&n);
 printf("o=");
 scanf("%lg",&o);
 printf("p=");
 scanf("%lg",&p);
 printf("°_©l­È=");
 scanf("%lg",&a0);
 printf("¦¬´M¨B¥ï=");
 scanf("%lg",&h1);
 printf("fx=%g*(%g*x1+%g*x2)^%g+%g*(%g*x3-%g*x4)^%g+%g*(%g*x2-%g*x3)^%g+%g*(%g*x1+%g*x4)^%g\n",a,b,c,d,e,f,g,h,i,j,k,l,m,n,o,p);
 x1=a0;
 x2=a0;
 x3=a0;
 x4=a0;
 w1=x1+10*x2;
 w2=x3-x4;
 w3=x2-2*x3;
 w4=x1-x4;
 f1=p*(pow(w1,c)+d*pow(w2,g)+h*pow(w3,k)+l*pow(w4,o));
 x1=a0+h1;
 x2=a0+h1;
 x3=a0+h1;
 x4=a0+h1;
 f2=p*(pow(w1,c)+d*pow(w2,g)+h*pow(w3,k)+l*pow(w4,o));
  if (f2<f1)
  {
    aa:
       h1=2*h1;
       a2=a2+h1;
       f1=f2;
       f2=f2;
          if(f2<f1)
          {
            a=a1;
            b=a2;
          }
          else
          {
            a1=a2+h1;
            goto aa;
          }
  }
  else
  {
    h1=-h1/4;
    bb:
    a1=a1+h1;
    f2=f1;
    f1=f1;

        if (f2<=f1)
        {
          a=a1;
          b=a2;
        }
        else
        {
         a2=a1+h1;
         h1=2*h1;
         goto bb;
        }
printf("ºë«×­ÈE=");
scanf("%lg",&E);
b1=b-E*(b-a);
x1=b1;
x2=b1;
x3=b1;
x4=b1;
w1=x1+10*x2;
w2=x3-x4;
w3=x2-2*x3;
w4=x1-x4;
f1=p*(pow(w1,c)+d*pow(w2,g)+h*pow(w3,k)+l*pow(w4,o));
cc:
   b2=a+k1*(b-a);
   x1=b2;
   x2=b2;
   x3=b2;
   x4=b2;
   w1=x1+10*x2;
   w2=x3-x4;
   w3=x2-2*x3;
   w4=x1-x4;
   f2=p*(pow(w1,c)+d*pow(w2,g)+h*pow(w3,k)+l*pow(w4,o));
dd:
 if(fabs((f2-f1)/f1)<E)
 {
  if(fabs((b2-b1)/b1)<E)
   {
     if(f1<f2)
      {
       printf("f1");
      }
       else
       {
       printf("f2");
       }
   }
  else
    {
     a=b1;
     b=b2;
     goto cc;
    }
  else
  {
  if (f1>f2)
   {
    a=b1;
    b1=b2;
    f1=f2;
    b2=a+k1*(b-a);
    goto cc;
   }
   else
    {
     b=b2;
     b1=b2;
     f1=f2;
     b1=b-k1*(b-a);
     x1=b1;
     x2=b1;
     x3=b1;
     x4=b1;
     w1=x1+10*x2;
     w2=x3-x4;
     w3=x2-2*x3;
     w4=x1-x4;
     f1=p*(pow(w1,c)+d*pow(w2,g)+h*pow(w3,k)+l*pow(w4,o));
     goto dd;
    }
  }
 }
        //getch();

 }
 return 0;
}