T 1 T_1 T1——seq(3018)
Description:
有一个长度为
n
n
n的序列
A
A
A,每次操作可以将相邻的两个数
A
i
,
A
i
+
1
A_i,A_{i+1}
Ai,Ai+1删去,在这个它们位置之间插入
m
a
x
{
A
i
,
A
i
+
1
}
max\{A_i,A_{i+1}\}
max{Ai,Ai+1},而这样操作的代价为
m
a
x
{
A
i
,
A
i
+
1
}
max\{A_i,A_{i+1}\}
max{Ai,Ai+1}。求将序列
A
A
A的长度变为
1
1
1的最小代价。
n
≤
1
0
6
,
A
i
≤
1
0
9
n\le 10^6,A_i\le 10^9
n≤106,Ai≤109
Solution:
- 奇妙的贪心…(然而考试时只是无脑ST表,结果MLE 5分了)
- 待证明…
Code:
#include<bits/stdc++.h>
using namespace std;
#define REP(i,f,t) for(int i=(f),i##_end_=(t);i<=i##_end_;++i)
#define ll long long
template<class T>inline void Rd(T &x){
x=0;char c;
while((c=getchar())<48);
do x=(x<<1)+(x<<3)+(c^48);
while((c=getchar())>47);
}
const int N=1e6+2;
int n;
int A[N];
ll ans;
int main(){
// freopen("seq.in","r",stdin);
// freopen("seq.out","w",stdout);
Rd(n);
REP(i,1,n){
Rd(A[i]);
if(i>1) ans+=max(A[i],A[i-1]);
}
printf("%lld\n",ans);
return 0;
}
T 2 T_2 T2——install(3019)
Description:
有
n
n
n个软件和总容量为
m
m
m的磁盘,每个软件有三个属性容量
w
i
w_i
wi,价值
v
i
v_i
vi,依赖
d
i
d_i
di,这里依赖是指对于
i
i
i来说,只有
d
i
d_i
di安装了,才会产生价值。求最大价值。
n
≤
100
,
m
≤
500
,
v
i
≤
1000
n\le 100,m\le500,v_i\le 1000
n≤100,m≤500,vi≤1000
Solution:
- 首先,对于都没有依赖的情况,就是一个裸的背包问题
- 而对于有依赖的情况,很明显是有环的情况,那么我们就可以tarjan缩点,使之为一个DAG,这样就与树上的情况是一样的。
- 那么就是树上的背包, d p [ i ] [ j ] dp[i][j] dp[i][j]表示以 i i i为子树的用了体积 j j j的最大价值。
Code:
#include<bits/stdc++.h>
using namespace std;
#define REP(i,f,t) for(int i=(f),i##_end_=(t);i<=i##_end_;++i)
#define SREP(i,f,t) for(int i=(f),i##_end_=(t);i<i##_end_;++i)
#define DREP(i,f,t) for(int i=(f),i##_end_=(t);i>=i##_end_;--i)
#define ll long long
template<class T>inline bool chkmax(T &x,T y){return x<y?x=y,1:0;}
template<class T>inline bool chkmin(T &x,T y){return x>y?x=y,1:0;}
const int N=502,M=1002;
int n,m;
int w[N<<1],v[N<<1],fa[N];
struct p20{
int t[22];
void solve(){
int ans=0;
SREP(s,0,1<<n){
memset(t,0,sizeof(t));
t[0]=1;
int sum=0,res=0;
REP(i,1,n) if(s&(1<<(i-1))) t[i]=1;
REP(i,1,n){
if(t[i]){
sum+=w[i];
if(t[fa[i]]) res+=v[i];
}
}
if(sum<=m) chkmax(ans,res);
}
printf("%d\n",ans);
}
}p1;
struct p50{
int dp[502];
void solve(){
REP(i,1,n) DREP(j,m,w[i]) chkmax(dp[j],dp[j-w[i]]+v[i]);
int ans=0;
REP(i,0,m) chkmax(ans,dp[i]);
printf("%d\n",ans);
}
}p2;
struct p100{
int head[N],qwq;
struct edge{
int to,nxt;
}E[N<<2];
void addedge(int x,int y){E[qwq]=(edge){y,head[x]};head[x]=qwq++;}
int dfn[N<<1],low[N<<1],tim;
int stk[N<<1],top;
int Id[N<<1],tot;
bool vis[N<<1];
int In[N<<1];
int W[N],V[N];
int dp[N<<1][M];
void tarjan(int x){
dfn[x]=low[x]=++tim;
stk[++top]=x;
vis[x]=1;
for(int i=head[x];~i;i=E[i].nxt){
int y=E[i].to;
if(!dfn[y]) tarjan(y),chkmin(low[x],low[y]);
else if(vis[y]) chkmin(low[x],dfn[y]);
}
if(dfn[x]==low[x]){
tot++;
do{
Id[stk[top]]=tot;
vis[stk[top]]=0;
W[tot]+=w[stk[top]];
V[tot]+=v[stk[top]];
}while(x!=stk[top--]);
}
}
void dfs(int x){
REP(i,W[x],m) dp[x][i]=V[x];
for(int i=head[x];~i;i=E[i].nxt){
int y=E[i].to;
dfs(y);
DREP(k,m-W[x],0) REP(j,0,k) chkmax(dp[x][k+W[x]],dp[y][j]+dp[x][k+W[x]-j]);
}
}
void solve(){
qwq=0;
memset(head,-1,sizeof(head));
REP(i,1,n) if(fa[i]) addedge(fa[i],i);
REP(i,1,n) if(!dfn[i]) tarjan(i);
qwq=0;
memset(head,-1,sizeof(head));
REP(i,1,n) if(Id[fa[i]]!=Id[i]) addedge(Id[fa[i]],Id[i]),In[Id[i]]++;
REP(i,1,tot) if(!In[i]) addedge(0,i);
dfs(0);
printf("%d\n",dp[0][m]);
}
}p3;
int main(){
// freopen("install.in","r",stdin);
// freopen("install.out","w",stdout);
scanf("%d%d",&n,&m);
REP(i,1,n) scanf("%d",&w[i]);
REP(i,1,n) scanf("%d",&v[i]);
bool flag=0;
REP(i,1,n) scanf("%d",&fa[i]),flag|=fa[i];
if(!flag) p2.solve();
else p3.solve();
return 0;
}
T 3 T_3 T3——book(1763)
Description:
有
n
n
n个书架,每个书架有
m
m
m个位置,对于每个位置都放着不同的书,现在有两种操作:(1)若在某个书架上,某书的左边或右边为空,则可以将该书左移或右移;(2)可以将一本书拿出放在空位置上。
现在给出书架的初态和末态,求需要最少的操作2次数。
若无法完成,即输出-1。
n
,
m
≤
1000
n,m\le 1000
n,m≤1000,书的编号
≤
n
⋅
m
\le n\cdot m
≤n⋅m
Solution:
- 首先,对于操作1的情况,如果它的初态和末态的相对位置是不一样的,那么一定要换,可能需要操作2
- 模拟小数据发现,只需要关心该行是否还需要操作2,而需要操作2一定要它的初态和末态的相对位置不同。
- 而真正需要操作2的,即初态和末态不在同一个行时,我们可以建图跑环,发现只要有存在不同的位置,那么一定要进行一次操作2,
- 而末态位置的书架是满的,那么一定还要一次操作将该书架空出一个位置来。
Code:
#include<bits/stdc++.h>
using namespace std;
#define REP(i,f,t) for(int i=(f),i##_end_=(t);i<=i##_end_;++i)
#define SREP(i,f,t) for(int i=(f),i##_end_=(t);i<i##_end_;++i)
#define DREP(i,f,t) for(int i=(f),i##_end_=(t);i>=i##_end_;--i)
#define ll long long
template<class T>inline bool chkmax(T &x,T y){return x<y?x=y,1:0;}
template<class T>inline bool chkmin(T &x,T y){return x>y?x=y,1:0;}
template<class T>inline void Rd(T &x){
x=0;char c;
while((c=getchar())<48);
do x=(x<<1)+(x<<3)+(c^48);
while((c=getchar())>47);
}
const int N=1002;
int n,m;
int A[N][N],B[N][N];
typedef pair<int,int>pii;
#define fi first
#define se second
struct pwater{
pii pos1[N*N],pos2[N*N];
vector<pii>row[N];
vector<int>tmp;
vector<int>::iterator it;
int cnt[N];
int Mxid;
int qwq,head[N*N];
struct edge{
int to,nxt;
}E[(N*N)<<1];
void addedge(int x,int y){E[qwq]=(edge){y,head[x]};head[x]=qwq++;}
bool vis[N];
int move(int x){//看相对位置
tmp.clear();
sort(row[x].begin(),row[x].end());
SREP(i,0,row[x].size()){
it=lower_bound(tmp.begin(),tmp.end(),row[x][i].se);
if(it==tmp.end()) tmp.push_back(row[x][i].se);
else *it=row[x][i].se;
}
int res=0;
if(tmp.size()==row[x].size()) return 0;
if((int)row[x].size()==m) res=1;
res+=row[x].size()-tmp.size();
return res;
}
int dfs(int x){
int res=(cnt[x]==m);
vis[x]=1;
for(int i=head[x];~i;i=E[i].nxt){
int y=E[i].to;
if(vis[y])continue;
res&=dfs(y);
}
return res;
}
void solve(){
bool flag=1;
bool zero=0;
REP(i,1,n) REP(j,1,m) {
if(A[i][j]!=B[i][j]) flag=0;
if(!A[i][j]) zero=1;
pos1[A[i][j]]=(pii){i,j};
pos2[B[i][j]]=(pii){i,j};
chkmax(Mxid,A[i][j]);
cnt[i]+=(A[i][j]>0);
}
if(flag){puts("0");return;}
if(!zero){puts("-1");return;}
memset(head,-1,sizeof(head));
int ans=0;
REP(i,1,Mxid){
if(pos1[i].fi==pos2[i].fi){//同一行 可能可以移动
row[pos1[i].fi].push_back((pii){pos1[i].se,pos2[i].se});
}
else {//至少要拿出
addedge(pos1[i].fi,pos2[i].fi);
addedge(pos2[i].fi,pos1[i].fi);
++ans;
}
}
REP(i,1,n) ans+=move(i);
REP(i,1,n) if(!vis[i] && ~head[i]) ans+=dfs(i);
printf("%d\n",ans);
}
}p2;
int main(){
// freopen("book.in","r",stdin);
// freopen("book.out","w",stdout);
Rd(n),Rd(m);
REP(i,1,n) REP(j,1,m) Rd(A[i][j]);
REP(i,1,n) REP(j,1,m) Rd(B[i][j]);
p2.solve();
return 0;
}