woj4768~woj4770
名字长了点,,不过的确是三个线段树。
见这个博客T3。
签到题。分两个集合。煮熟和没煮熟。同时还要分id。我们用优先队列区分是否煮熟。用树状数组存储每种有多少个煮熟的还在锅里。用vector装每个食物的加入序列(自然有序),然后用isout记录每种食物煮熟了几个,这样可以知道vector当前到哪一位了。
然后模拟。中间有个二分。
#include<bits/stdc++.h>
using namespace std;
#define in read()
#define int long long
int mx;
int in{
int cnt=0,f=1;char ch=0;
while(!isdigit(ch)){
ch=getchar();if(ch=='-')f=-1;
}
while(isdigit(ch)){
cnt=cnt*10+ch-48;
ch=getchar();
}return cnt*f;
}
struct node{int id,tim;};
bool operator <(node a,node b){return a.tim>b.tim;}
int T,n,qq;
priority_queue<node> q;
vector<int> v[1000003];
int isout[1000003],s[1000003];
int t[1000003];
int lowbit(int x){return x&(-x);}
void modify(int x,int key){while(x<=mx){t[x]+=key;x+=lowbit(x);}}
int query(int x){
int sum=0;
while(x){sum+=t[x];x-=lowbit(x);}
return sum;
}
void workyazid(){
if(!t[mx]){cout<<"Yazid is angry."<<'\n';return;}
int l=0,r=mx,now=1;
while(l<r){
int mid=(l+r)>>1;
if(t[mid]<now)l=mid+1,now-=t[mid];
else r=mid;
}
modify(l,-1);cout<<l<<'\n';
}
void workyzzkal(int id,int t){
if(query(id)-query(id-1)){modify(id,-1);cout<<"Succeeded!"<<'\n';return;}
if(isout[id]==v[id].size()){cout<<"YJQQQAQ is angry."<<'\n';return;}
cout<<v[id][isout[id]]-t<<'\n';
}
signed main(){
//freopen("b.in","r",stdin);freopen("b.out","w",stdout);
T=in;
while(T--){mx=1;
memset(t,0,sizeof(t));memset(isout,0,sizeof(isout));
n=in;for(int i=1;i<=n;i++)v[i].clear(),isout[i]=0,t[i]=0,s[i]=in;
while(!q.empty())q.pop();qq=in;
while(mx<=n)mx<<=1;
while(qq--){
int t=in;int op=in;
while(!q.empty()){if(q.top().tim>t)break;isout[q.top().id]++;modify(q.top().id,1);q.pop();}
if(op==0){int id=in;v[id].push_back(s[id]+t);q.push((node){id,t+s[id]});}
if(op==1){workyazid();}
if(op==2){int id=in;workyzzkal(id,t);}
if(op==3){int l=in;int r=in;int gu=query(r)-query(l-1);cout<<gu<<'\n';}
}
}
return 0;
}
这道题有zxyoi的n方log的好算法。不过还有更骚的O(n)的结论方法。
结论:对于一个三角剖分,将有边相邻的三角形连边,可以形成一棵树。这棵树上的任何一个简单路径所代表的三角形,都可以被一条连接两个端点的直线穿过且仅穿过这些三角形。
所以这道题转化成了,点权为1,求一个大概是“人”形的路线,就是说允许分叉一次,的最大值。。
树形动规即可,,最大次大值。
#include<bits/stdc++.h>
using namespace std;
#define in read()
int in{
int cnt=0,f=1;char ch=0;
while(!isdigit(ch)){
ch=getchar();if(ch=='-')f=-1;
}
while(isdigit(ch)){
cnt=cnt*10+ch-48;
ch=getchar();
}return cnt*f;
}
int cnt[1003][1003];
int s[1003][1003][2];
struct node{
int a,b,c;
}a[1003];
int n;
int gu[4];
int first[1003],nxt[2003],to[2003],tot;
void add(int a,int b){
//cout<<a<<" "<<b<<endl;
nxt[++tot]=first[a];first[a]=tot;to[tot]=b;
}int d[1003],ans;
void dfs(int u,int fa){
for(int i=first[u];i;i=nxt[i]){
int v=to[i];if(v==fa)continue;dfs(v,u);d[u]=max(d[u],d[v]);
}d[u]++;
}
void dfs2(int u,int fa,int out){
int mx[2]={0,0};
for(int i=first[u];i;i=nxt[i]){
int v=to[i];if(v==fa)continue;
if(d[v]>mx[0])mx[1]=mx[0],mx[0]=d[v];
else if(d[v]>mx[1])mx[1]=d[v];
}ans=max(ans,mx[0]+mx[1]+1+out);
for(int i=first[u];i;i=nxt[i]){
int v=to[i];if(v==fa)continue;
dfs2(v,u,1+max(out,(d[v]==mx[0])?mx[1]:mx[0]));
}
}
signed main(){
n=in;
for(int i=1;i<=n-2;i++){
gu[1]=in+1;gu[2]=in+1;gu[3]=in+1;sort(gu+1,gu+4);a[i]={gu[1],gu[2],gu[3]};
//cout<<"# "<<gu[1]<<" "<<gu[2]<<" "<<gu[3]<<endl;
s[gu[1]][gu[2]][cnt[gu[1]][gu[2]]++]=i;s[gu[1]][gu[3]][cnt[gu[1]][gu[3]]++]=i;s[gu[2]][gu[3]][cnt[gu[2]][gu[3]]++]=i;
}
for(int i=1;i<n;i++){
for(int j=i+1;j<=n;j++){
//cout<<i<<" "<<j<<" "<<s[i][j][0]<<" "<<s[i][j][1]<<endl;
if(!s[i][j][0]||!s[i][j][1])continue;
add(s[i][j][0],s[i][j][1]);add(s[i][j][1],s[i][j][0]);
}
}
dfs(1,0);dfs2(1,0,0);cout<<ans;
return 0;
}