BZOJ4240: 有趣的家庭菜园
Aug 18, 2015 08:32:43 PM
题解:
首先预处理出$l_i$表示第$i$个数左侧比第$i$个数大的个数,$r_i$右侧同理.
首先脑补出我们应该将序列排成先单调不降,后单调不增的形式.
我们从大到小将数插入序列,考虑插入某个数,比这个数大的已经全部插入完毕了,由于排成那种形式,我们只能将这个数放在当前序列的头或尾.
若放在开头,则带来$l_i$的逆序对;否则带来$r_i$的逆序对.
显然每个数之间都是独立的,我们直接贪心就行了.
时间复杂度$O(nlogn)$.
代码:
#include<cstdio>
#include<cstring>
#include<cctype>
#include<iostream>
#include<algorithm>
using namespace std;
typedef long long ll;
#define N 300010
int n,_n,a[N],b[N];
inline int getc(){
static const int L=1<<15;
static char buf[L],*S=buf,*T=buf;
if(S==T){
T=(S=buf)+fread(buf,1,L,stdin);
if(S==T)
return EOF;
}
return*S++;
}
inline int getint(){
int c;
while(!isdigit(c=getc()));
int x=c-'0';
while(isdigit(c=getc()))
x=(x<<1)+(x<<3)+c-'0';
return x;
}
int A[N],l[N],r[N];
inline int ask(int x){
int re=0;
for(;x;x-=x&-x)
re+=A[x];
return re;
}
inline void upd(int x){
for(;x<=_n;x+=x&-x)
++A[x];
}
int main(){
n=getint();
int i,j;
for(i=1;i<=n;++i)
a[i]=b[i]=getint();
sort(b+1,b+n+1);
_n=unique(b+1,b+n+1)-b-1;
for(i=1;i<=n;++i)
a[i]=lower_bound(b+1,b+_n+1,a[i])-b;
for(i=1;i<=n;++i){
l[i]=i-1-ask(a[i]);
upd(a[i]);
}
memset(A,0,sizeof A);
for(i=n;i>=1;--i){
r[i]=n-i-ask(a[i]);
upd(a[i]);
}
long long ans=0;
for(i=1;i<=n;++i)
ans+=min(l[i],r[i]);
cout<<ans<<endl;
return 0;
}
BZOJ3174:[Tjoi2013]拯救小矮人 贪心+dp
Jan 23, 2015 04:17:52 PM
思路:
假设在能逃出相同多数目的小矮人的情况下,我们更愿意让逃生能力更强的小矮人留到最后.
逃生能力就是指\(a_i+b_i\),也即在下面垫的总高度确定时更容易达到更高的高度.
这样我们就有了一个贪心选择顺序,我们设\(f[i]\)表示逃生了\(i\)个小矮人时下面垫的总高度的最大值,并依据这个判定能否再逃出一人即可.
#include<cstdio>
#include<cstring>
#include<cctype>
#include<iostream>
#include<algorithm>
using namespace std;
#define N 2010
struct Case{int a,b;bool operator<(const Case&B)const{return a+b<B.a+B.b;}}S[N];
int f[N];
int main(){
int n,m;scanf("%d",&n);register int i,j;for(i=1;i<=n;++i)scanf("%d%d",&S[i].a,&S[i].b);
sort(S+1,S+n+1);
scanf("%d",&m);
memset(f,-1,sizeof f);for(f[0]=0,i=1;i<=n;++i)f[0]+=S[i].a;
int ans=0;
for(i=1;i<=n;++i){
for(j=ans;j>=0;--j){
if(f[j]+S[i].b>=m)f[j+1]=max(f[j+1],f[j]-S[i].a);
}if(f[ans+1]>=0)++ans;
}
printf("%d\n",ans);
return 0;
}
BZOJ2799:[Poi2012]Salaries 贪心
Jan 23, 2015 03:54:24 PM
思路:
我们维护以每一个节点为根的子树内部有多少个没有确定权值的点,另外,我们利用一个数组维护每个权值已经给了哪个节点.
我们从小到大遍历每个权值,若这个权值当前并没有确定给哪个节点,则将其压入栈中;否则我们在对应节点的子树中进行一系列确定操作:
我们维护栈中有多少个权值,以及有多少个自由权值.(这个是什么一会再说)
若未确定的节点数正好等于当前的自由权值数与栈中的权值总数,显而易见这些权值都应该分配给这棵子树.但是只有只有一个儿子的情况才能确定这个权值.所以我们不断在链上去找即可.这个过程可能确定了一些权值,随后我们将栈清空,并令自由权值的数目为0.因为他们只能被拘束在这颗子树内,而不能去别的地方.
若不然,我们将自由权值及栈中的权值分配给这棵子树所需要的数目,剩下的全部变成自由权值.
我们发现栈中的权值和自由权值很像,但是为什么不都搞成自由权值,而是非要维护一个栈呢?
这个问题显而易见:自由权值都是可以互相交换的(等价的),而栈中的元素显然不是!
于是\(O(n)\)水.
#include<cstdio>#include<cstring>#include<cctype>#include<iostream>#include<algorithm>#include<vector>using namespace std;namespace Fio{ inline int getc(){ static const int L=1<<15;static char buf[L],*S=buf,*T=buf; if(S==T){T=(S=buf)+fread(buf,1,L,stdin);if(S==T)return EOF;} return*S++; } template<typename T>inline void Get(T&x){ int c;while(!isdigit(c=getc()));x=c-'0';while(isdigit(c=getc()))x=(x<<1)+(x<<3)+c-'0'; }}#define N 1000010vector<int>v[N];int pa[N],w[N];int q[N],fr,ta,ins[N],size[N],stk[N],top,tot,sum;int main(){ int n;Fio::Get(n); register int i,j; for(i=1;i<=n;++i){ Fio::Get(pa[i]),Fio::Get(w[i]);if(i==pa[i])w[i]=n; if(w[i])q[ta++]=i,ins[w[i]]=i;else v[pa[i]].push_back(i); } int u; while(fr<ta){ u=q[fr++]; for(j=0;j<v[u].size();++j)q[ta++]=v[u][j]; } for(i=n-1;i>=0;--i){ u=q[i]; size[u]=1; for(j=0;j<v[u].size();++j)size[u]+=size[v[u][j]]; } int tot=0; for(i=1;i<=n;++i){ if(!(u=ins[i])){stk[++top]=i;continue;} int sum=0,ac=0,pretop=top; for(j=0;j<v[u].size();++j)sum+=size[v[u][j]]; if(sum==tot+top){ for(;v[u].size()==1&⊤u=v[u][0]) ins[w[v[u][0]]=stk[top--]]=v[u][0],++ac,--sum; if(sum)tot=top=0; } else{ if(sum)tot=tot+top-sum,top=0; } } for(i=1;i<=n;++i)printf("%d\n",w[i]); return 0;}