杭州网赛的1005,题意挺简单,每个球有x,y,z三个坐标,球i得三个坐标同事小于等于球j的坐标时,可以用i去撞击j,撞击后i消失,j继续去撞下一个球。问最多可以撞多少个球和有多少种不同的方法。比赛结束后也一直没什么好的思路,刚开始想着用拓扑排序去搞结果发现极限情况的边太多记录不下..后来看了 http://troysky.com/hdu-4742-pinball-game-3d/ 这篇题解,又看了看一些大神的代码,大致是明白怎么做了。
总体思路其实还是个Dp,dp[k]是一个pair<int,int>用来表示当前激活的球是k时,撞击的球数和撞击的方法数.具体做法用的是归并排序的思想,假定现在要搞区间[l,r],(l,r)表示当前要处理的球的编号是l到r (这里的球编号可以看做x坐标的哈希值,编号大的球x坐标一定不小于编号小的球)),区间[l,mid]已经处理完毕,那么对区间l,r按照y坐标排个序,这样y也变得有序了,之后还剩下z坐标,这里可以用线段树或者树状数组去维护一个最优值,对于排序以后的序列,从前到后依次遍历,若dt[i].id小于等于mid,那么dp[dt[a[i].id]]已经处理过了,直接按dt[i].z的值更新到树状数组里即可,对于dt[i].id大于mid的球,根据dt[i].z,需要在0--dt[i].z中中找到dp值最大的,直接用树状数组查询即可。之后递归去搞区间[mid+1,r]就ok了。因为用到了树状数组,所以读取数据的时候要对z坐标做一个哈希。之后再按x为第一关键字排个序,记录一下每个球的id。
代码如下:
#include <iostream>
#include <algorithm>
#include <cstdio>
#include <memory.h>
#include <cmath>
#include <algorithm>
#include <queue>
#include <stack>
#include <map>
using namespace std;
typedef long long ll;
const ll mod=(1<<30)-1;
const int maxn=100000+100;
inline int getint()
{
char c=getchar();
int t=0;
while (c<'0' || c>'9') c=getchar();
while (c >='0' && c<='9')
{
t=t*10+c-'0';
c=getchar();
}
return t;
}
int tt,n,m,p,q;
struct node
{
int x,y,z,id;
node(){}
node(int a,int b,int c)
{
x=a;
y=b;
z=c;
}
bool operator<(const node& p)const
{
// return x<p.x||(x==p.x && y<p.y)||(x&&p.x&&y==p.y && z<p.z);
if (x!=p.x) return x<p.x;
if (y!=p.y) return y<p.y;
return z<p.z;
}
}dt[maxn],tp[maxn];
int valz[maxn];
typedef pair<int,int> pii;
pii dp[maxn];
void merge(pii& p,pii q)
{
if (p.first<q.first) p=q;
else if (p.first==q.first)
{
p.second+=q.second;
p.second=p.second&mod;
}
}
struct BIT
{
pii c[maxn];
void init()
{
memset(c,0,sizeof c);
}
inline int lowbit(int x)
{
return x&-x;
}
void add(int pos,pii t)
{
for(;pos<=n; pos+=lowbit(pos))
{
merge(c[pos],t);
}
}
void clear(int pos)
{
for(;pos<=n;pos+=lowbit(pos))
{
c[pos]=pii(0,0);
}
}
pii query(int pos)
{
pii ret(0,0);
for(;pos; pos-=lowbit(pos))
{
merge(ret,c[pos]);
}
return ret;
}
}bit;
void slove(int l,int r)
{
if (l==r) return;
int mid=(l+r)>>1;
slove(l,mid);
int sp=0;
for (int i=l; i<=r;i++)
{
tp[sp]=dt[i];
tp[sp].x=0;
sp++;
}
sort(tp,tp+sp);
for (int i=0; i<sp; i++)
{
if(tp[i].id<=mid)
bit.add(tp[i].z,dp[tp[i].id]);
else
{
pii res=bit.query(tp[i].z);
res.first++;
merge(dp[tp[i].id],res);
}
}
for (int i=0; i<sp; i++)
{
if (tp[i].id<=mid) bit.clear(tp[i].z);
}
slove(mid+1,r);
}
int main()
{
// freopen("in.txt","r",stdin);
tt=getint();
while(tt--)
{
bit.init();
n=getint();
int tolz=0;
for (int i=0; i<n; i++)
{
dt[i].x=getint();
dt[i].y=getint();
dt[i].z=getint();
valz[tolz++]=dt[i].z;
}
sort(valz,valz+tolz);
tolz=unique(valz,valz+tolz)-valz;
for (int i=0; i<n;i++)
{
dt[i].z=lower_bound(valz,valz+tolz,dt[i].z)-valz+1;
}
sort(dt,dt+n);
for (int i=0;i<n;i++)
{
dt[i].id=i;
}
for (int i=0;i<=n;i++)
dp[i]=pii(1,1);
slove(0,n-1);
pii ans(0,0);
for (int i=0; i<n; i++)
merge(ans,dp[i]);
printf("%d %d\n",ans.first,ans.second);
}
return 0;
}