3193: [JLOI2013]地形生成


Time Limit: 10 Sec  Memory Limit: 128 MB
Submit: 459  Solved: 223
[Submit][Status][Discuss]

Description


 
 
最近IK正在做关于地形建模的工作。其中一个工作阶段就是把一些山排列成一行。每座山都有各不相同的标号和高度。为了遵从一些设计上的要求,每座山都设置了一个关键数字,要求对于每座山,比它高且排列在它前面的其它山的数目必须少于它的关键数字。
 显然满足要求的排列会有很多个。对于每一个可能的排列,IK生成一个对应的标号序列和等高线序列。标号序列就是按顺序写下每座山的标号。等高线序列就是按顺序写下它们的高度。例如有两座山,这两座山的一个合法排列的第一座山的标号和高度为1和3,而第二座山的标号和高度分别为2和4,那么这个排列的标号序列就是1 2,而等高线序列就是3 4.
 现在问题就是,给出所有山的信息,IK希望知道一共有多少种不同的符合条件的标号序列和等高线序列。
 

Input


输入第一行给出山的个数N。接下来N行每行有两个整数,按照标号从1到N的顺序分别给出一座山的高度和关键数。
 

Output


 
输出两个用空格分隔开的数,第一个数是不同的标号序列的个数,第二个数是不同的等高线序列的个数。这两个答案都应该对2011取模,即输出两个答案除以2011取余数的结果
 

Sample Input



Sample Output


 

HINT


对于所有的数据,有1<=N<=1000,所有的数字都是不大于109的正整数。

分析:


先考虑第一个问题(先假设没有山高度相同):

我们把所有山按高度从大到小排序,设这个集合为S,设一个空集合为V。

我们把山按顺序放入集合V。

考虑第i座山放进去的时候能放的位置,因为现在前面有i - 1座山比它高,再加上本来的第i个位置,它一共有i个位置可以放。这个山的关键值如果为k,它只能在前min(k,i)个位置里挑。

所以第i个山的放的位置的组合为ci = min(k,i)种。ans = c1 * c2 …… * cn;

但是我们现在有山的高度相同,且关键值不同。两座高度相同的山能放的最靠后的位置,一定是关键值大的那个越靠后。所以我们先按高度从大到小,然后高度相同的关键值从小到大排序出S集合。

设[x,y]这段区间的山高度相同。所以 ci = min(x,k) + i - num。因为一座山在满足关键值的情况下,可以也放在高度相同的山前面。

所以最后 得出 ans1 = c1 * c2 …… * cn;

然后再看第二个问题,会产生重复情况的只有高度相同的山放一起的情况。

我们依然排序出V集合。按顺序插入山。

定义dp状态f[i][j],表示高度相同的山里面前i座,放在前j个位置的方案数。

因为我们是插入的山,比如说在一座高度为3的山后面再插入一座高度为3的山,此时是多了一种方案,是不会重复的。

状态也是很好转移设区间[x,y]山高度相同

f[i][j] = sum(f[i - 1][k] );k <= min(data[i].d,x) - 1;

因为是按顺序来刷表的可以省掉一维。

又因为f[i][j - 1] = sum(f[i][j - 1]);f[i][j] = sum(f[i - 1][k] );k <= min(data[i - 1].d,x) - 1;所以f[i][j] = f[i - 1][j] + f[i][j - 1]

最终方程dp[i] = dp[i - 1] + dp[i]

最后把每一区间[x,y]相同的山的ci = dp[0] + dp[1] +……dp[min(data[y].d,x) - 1]乘起来。

ans2 = c1 * c2 …… * cn;

就可以了。

AC代码:


# include <iostream>
# include <cstdio>
# include <cstring>
# include <algorithm>
using namespace std;
const int mod = ;
const int N = ;
int dp[N],h[N],n,ans1 = ,ans2 = ;
struct Mountain{
int h;
int d;
bool operator <(const Mountain & other)const{
if(h == other.h)return d < other.d;
return h > other.h;
}
}data[N];
void Init(){
memset(dp,,sizeof dp);
dp[] = ;
}
void read(){
scanf("%d",&n);
for(int i = ;i <= n;i++){
scanf("%d %d",&data[i].h,&data[i].d);
}
sort(data + ,data + n + );
int num = ;
for(int i = ;i <= n;i++){
if(data[i].h != data[num].h)num = i;
(ans1 *= min(num,data[i].d) + i - num) %= mod;
}
printf("%d",ans1);
}
void work(){
int pos;
for(int i = ;i <= n;i++){
pos = i;
while(data[pos].h == data[i].h && pos <= n)pos++;
pos--;Init();
for(int j = i;j <= pos;j++){
for(int k = ;k <= min(data[j].d,i) - ;k++){
(dp[k] = dp[k - ] + dp[k]) %= mod;
}
}
int sum = ;
for(int j = ;j <= min(data[pos].d,i) - ;j++)(sum += dp[j]) %= mod;
(ans2 *= sum) %= mod;
i = pos;
}
printf(" %d\n",ans2);
}
int main()
{
read();
work(); }

最新文章

  1. fragment 重叠问题
  2. Convert
  3. 【转载】javascript与C#的语法区别
  4. 《所用到的AJAX技术基础》
  5. matlab 扩大虚拟内存
  6. Swift中可选类型(Optional)的用法 以及? 和 ! 的区别 (转载博客,知识分享)
  7. SGU 135.Drawing Lines
  8. swift 重载 泛式 inout的使用
  9. 使用Post方法模拟登陆爬取网页
  10. C#删除区域实现透明
  11. css3的动画特效--元素旋转
  12. 3. Java面向对象之泛型-指定多个泛型
  13. Laravel-任务调度
  14. Description Resource Path Location Type Java compiler level does not match the version of(编译问题)
  15. python机器可读数据-csv
  16. 【Spark 深入学习-08】说说Spark分区原理及优化方法
  17. Loj10154 选课
  18. 第一个map reduce程序
  19. java的类继承(与c++对比)
  20. 转Google Protocol Buffer 的使用和原理

热门文章

  1. 什么是Entitlement
  2. SQL数据库移植到ARM板步骤
  3. js运行机制(线程)
  4. 网新恩普(T 面试)
  5. DLL动态链接库的创建
  6. PRJ0003 : Error spawning &#39;midl.exe&#39;
  7. linux 搜索文本
  8. js hover 下拉框
  9. Windows 命令收集
  10. 枚举为何不能设置成public?