「Usaco2008 Jan」人工湖O(∩_∩)O 纯属的模拟+栈

题目描述

夏日那让人喘不过气的酷热将奶牛们的烦躁情绪推到了最高点．最终，约翰决定建一个人工湖供奶牛消暑之用． 
为了使湖看起来更加真实，约翰决定将湖的横截面建成N(1≤N≤105)个连续的平台高低错落的组合状，所有的平台从 左到右按1到N依次编号．当然咯，在湖中注入水后，这些平台都将被淹没．  
 平台i在设计图上用它的宽度wi(1≤Wi≤1000)和高度（你可以理解为该平台顶离约翰挖的地基的高度）Hi.（1≤Hi≤1000000）来描述的． 
所有平台的高度都是独一无二的．湖的边缘可以视为无限高的平台． 
下面给出了一张约翰的设计图：  
 
按约翰的设想，在坑挖好后，他会以1单位每分钟的速度往最低的那个平台上注水．水在离开水管后立即下落，直到撞到平台顶或是更早些时候注入的水然后，与所有常温下的水一样，它会迅速地流动、扩散．简单起见，你可以认为这些都是在瞬间完成的． 
下图展示了一个样例，样例数据即是本题的所给输入样例。 

约翰想知道，对于每一个平台，它的顶部是从哪个时刻开始，与水面的距离至少为1单位长度． 
注意：数据不保证答案全部在32位整型变量的范围内 

输入

第1行:1个整数，N ; 
第2~N+1行:第i+1行为2个用空格隔开的整数Wi、Hi，描述了第i个平台. 

输出

输出N行。  
第1~N行: 第i行为1个整数，表示平台i的顶到水面的距离从何时开始大于1 单位长度。 

样例输入 Copy

3
4 2
2 7
6 4

样例输出 Copy

4
50
26

解题过程：

• 这题的主要算法就是栈，不得不说用栈搞模拟还挺神奇的。。。
• 这题就是模拟从最低的平台注水，然后将最低的填满以后从最低的平台向两边扩展。
• 每次找最近的最低的平台h，然后将水填到h高度。
• 栈里存的是向外扩展的时候，有时会遇到高度递减的情况，这时并不能填水，但要把这些高度都递减（即扩展时的顺序）记录进栈。
• 然后遇到一个比水面高的平台h时，模拟倒水，水会挨个淹没最低的平台，即需要从栈顶一个一个出栈计算淹没时间，直至栈顶平台高度>h，此时h入栈。重复执行就可算出答案。

然后就是代码，感觉模拟这个东西就是考码力和思维强度的。。。。

 // luogu-judger-enable-o2
 #include<bits/stdc++.h>
 using namespace std;
 const int maxn=1e6+;
 const int inf=<<;
 long long  top;
 struct Pool
 {
     long long w;
     long long h;
     long long id;
     Pool(long long ww=,long long hh=, long long _id=):w(ww),h(hh),id(_id){};
 }pool[maxn],struck[maxn];
 long long n,m;
 long long ans[maxn];
 long long now;
 inline long long  read(){
     register long long  x=,f=;char ch=getchar();
     while(!isdigit(ch)){if(ch=='-')f=-;ch=getchar();}
     while(isdigit(ch)){x=(x<<)+(x<<)+ch-'';ch=getchar();}
     return (f==)?x:-x;
 }
 int main()
 {
     long long temp=inf,id;
     n=read();
     for(int i=;i<=n;i++)
     {
         pool[i].w=read();
         pool[i].h=read();
         pool[i].id=i;
         if(temp>pool[i].h)
         {
             temp=pool[i].h;
             id=i;
         }
     }
     struck[++top]=pool[id];
     struck[].h=inf;
     long long l=id;
     long long r=id;
     long long p;
     pool[].h=pool[n+].h=inf;
     for(int i=;i<=n;i++)
     {
         long long add=;
         if(pool[l-].h<pool[r+].h)
         {
             p=--l;
         }
         else
         {
             p=++r;
         }
         while(pool[p].h>struck[top].h&&top)
         {
             struck[top].w+=add;
             ans[struck[top].id]=now+struck[top].w;
             now+=struck[top].w*(min(pool[p].h,struck[top-].h)-struck[top].h);
             add=struck[top].w;
             top--;
         }
         pool[p].w+=add;
         struck[++top]=pool[p];
     }
     for(int i=;i<=n;i++)
     {
         printf("%lld\n",ans[i]);
     }
     return ;
 }

总结一下：

在想到某些思路时要敢于打，模拟可以用各种各样的方式去优化，然后注意一下一些细节就可以了。