欢迎您访问程序员文章站本站旨在为大家提供分享程序员计算机编程知识!
您现在的位置是: 首页  >  IT编程

POJ1862 Stripies 贪心 B

程序员文章站 2022-05-30 21:58:24
POJ 1862 Stripies https://vjudge.net/problem/POJ-1862 题目: Our chemical biologists have invented a new very useful form of life called stripies (in fac ......

poj 1862 stripies

https://vjudge.net/problem/poj-1862

题目:

    our chemical biologists have invented a new very useful form of life called stripies (in fact, they were first called in russian - polosatiki, but the scientists had to invent an english name to apply for an international patent). the stripies are transparent amorphous amebiform creatures that live in flat colonies in a jelly-like nutrient medium. most of the time the stripies are moving. when two of them collide a new stripie appears instead of them. long observations made by our scientists enabled them to establish that the weight of the new stripie isn't equal to the sum of weights of two disappeared stripies that collided; nevertheless, they soon learned that when two stripies of weights m1 and m2 collide the weight of resulting stripie equals to 2*sqrt(m1*m2). our chemical biologists are very anxious to know to what limits can decrease the total weight of a given colony of stripies. 
    you are to write a program that will help them to answer this question. you may assume that 3 or more stipies never collide together. 

input

    the first line of the input contains one integer n (1 <= n <= 100) - the number of stripies in a colony. each of next n lines contains one integer ranging from 1 to 10000 - the weight of the corresponding stripie.

output

    the output must contain one line with the minimal possible total weight of colony with the accuracy of three decimal digits after the point.

sample input

3
72
30
50

sample output

120.000

分析:

贪心题目

题目本身不难,总共就几种策略,wa几发蒙也能蒙过了

 

问题在于这个题为什么这么做是正确的

我在做的时候就顺便用txt进行了证明,毕竟是练习,也为了保证一遍过

POJ1862 Stripies 贪心 B

首先写了这个,然后意识到不是相加,是直接把v1和v2变成另外一种,所以这是错误的

 

  POJ1862 Stripies 贪心 B

然后列了几种情况,当两个数的时候,产生的结果是不固定的,可能大,可能小,也可能与其中一个相等,这里就可以想到用增加变量然后用字母来代替数值进行分析

POJ1862 Stripies 贪心 B

v1v2v3

然后进行第一步,如果12合并,如果23合并,如果13合并

 POJ1862 Stripies 贪心 B

一行一个情况

在这其中保证

 POJ1862 Stripies 贪心 B

那么最后的三个结果就是

 POJ1862 Stripies 贪心 B

看起来比较复杂emmmm

那就简单化简吧,假设第一个和第二个相等

 POJ1862 Stripies 贪心 B

化简

 POJ1862 Stripies 贪心 B

然后取平方

 POJ1862 Stripies 贪心 B

诶是不是发现什么了

都会有

 POJ1862 Stripies 贪心 B

这个东西,

唯一的区别是v1,v2,v3

那么显然是

 POJ1862 Stripies 贪心 B

倒推:

 POJ1862 Stripies 贪心 B

得证

放全部分析过程:

POJ1862 Stripies 贪心 B

 

v1+v2 -> 2*sqrt(v1*v2)

sqrt(v1^2+v2^2+2*v1*v2) -> sqrt(4*v1*v2)

1 100

20

72 50

120

2 50

20

5 20

20

v1 v2 v3

 

2*sqrt(v1*v2) v3

2*sqrt(v1*v3) v2

2*sqrt(v2*v3) v1

v1<v2<v3

2*sqrt(2*sqrt(v1*v2)*v3)

2*sqrt(2*sqrt(v1*v3)*v2)

2*sqrt(2*sqrt(v2*v3)*v1)

2*sqrt(2*sqrt(v1*v2)*v3)=2*sqrt(2*sqrt(v1*v3)*v2)

sqrt(v1*v3)*v2=sqrt(v1*v2)*v3

v1*v2*v3*v3  v1*v3*v2*v2

v1*v2*v3

v1*v2*v3*v1 min

也就是sqrt(v2*v3)*v1

再倒推就是2*sqrt(2*sqrt(v2*v3)*v1)

再倒推就是先处理两个大的

得证

ac代码:

POJ1862 Stripies 贪心 B
 1 #include <stdio.h>
 2 #include <math.h>
 3 #include <string.h>
 4 #include <algorithm>
 5 #include <iostream>
 6 #include <string>
 7 #include <time.h>
 8 #include <queue>
 9 #include <string.h>
10 #define sf scanf
11 #define pf printf
12 #define lf double
13 #define ll long long
14 #define p123 printf("123\n");
15 #define pn printf("\n");
16 #define pk printf(" ");
17 #define p(n) printf("%d",n);
18 #define pln(n) printf("%d\n",n);
19 #define s(n) scanf("%d",&n);
20 #define ss(n) scanf("%s",n);
21 #define ps(n) printf("%s",n);
22 #define sld(n) scanf("%lld",&n);
23 #define pld(n) printf("%lld",n);
24 #define slf(n) scanf("%lf",&n);
25 #define plf(n) printf("%lf",n);
26 #define sc(n) scanf("%c",&n);
27 #define pc(n) printf("%c",n);
28 #define gc getchar();
29 #define re(n,a) memset(n,a,sizeof(n));
30 #define len(a) strlen(a)
31 #define ll long long
32 #define eps 1e-6
33 using namespace std;
34 double a[100500];
35 bool cmp(double a, double b){
36     return a > b;
37 }
38 int main() {
39     int n = 0;
40     s(n);
41     for(int i = 0; i < n; i ++){
42         slf(a[i]);
43     }
44     sort(a,a+n,cmp);
45     double temp = a[0];
46     for(int i = 1; i< n ;i ++){
47         temp = 2.0*sqrt(a[i]*temp);
48     }
49     pf("%.3lf\n",temp);
50     return 0;
51 }
view code