~~~~~ P1283 平板涂色 ~~~~~ 总题单链接
思路
~~~~~ 发现数据范围很小,考虑状压DP,用二进制下的 0 / 1 0/1 0/1表示每个矩形有没有涂色。
~~~~~ 设 d p [ i ] dp[i] dp[i] 为已经涂色的方块状态为 i i i 的情况下需要的最小次数。
~~~~~ 一种经典的转移就是枚举 i i i 的二进制下恰好少一个 1 的子状态,因为 i i i 的所有子状态都会转移成二进制下恰好少一个 1 1 1 的状态再转移。
~~~~~ 发现不太好转移,因为会根据(用来转移的子状态的最后一次涂色的颜色)和(当前涂色的矩形的颜色) 是否一样来分类,一样直接转移即可,而不一样要加 1 1 1。怎么办,再加上一维使得 d p [ i ] [ j ] dp[i][j] dp[i][j] 表示已经涂色的方块状态为 i i i ,且最后一次涂色的颜色为 j j j 的最小次数即可。
~~~~~
i
f
(
k
=
c
o
l
[
j
]
)
d
p
[
i
]
[
c
o
l
[
j
]
]
=
m
i
n
(
d
p
[
i
]
[
c
o
l
[
j
]
]
,
d
p
[
i
−
(
1
≪
j
)
]
[
k
]
)
if(k=col[j])dp[i][col[j]]=min(dp[i][col[j]],dp[i-(1\ll j)][k])
if(k=col[j])dp[i][col[j]]=min(dp[i][col[j]],dp[i−(1≪j)][k])
~~~~~
i
f
(
k
≠
c
o
l
[
j
]
)
d
p
[
i
]
[
c
o
l
[
j
]
]
=
m
i
n
(
d
p
[
i
]
[
c
o
l
[
j
]
]
,
d
p
[
i
−
(
1
≪
j
)
]
[
k
]
+
1
)
if(k\neq col[j])dp[i][col[j]]=min(dp[i][col[j]],dp[i-(1\ll j)][k]+1)
if(k=col[j])dp[i][col[j]]=min(dp[i][col[j]],dp[i−(1≪j)][k]+1)
代码
#include<bits/stdc++.h>
#define ll long long
#define INF 0x3f3f3f3f3f3f3f
using namespace std;
ll n,dp[100005][25];
vector<ll>vec[105],eg[25];
struct PIC{ll x1,y1,x2,y2,c;}pic[25];
inline ll lowbit(ll k){return k&-k;}
inline bool check(ll k){
for(ll j=0;j<n;j++){
if(k&(1ll<<j)){
for(ll v:eg[j])
if(!(k&(1ll<<v)))return 0;
}
}
return 1;
}
signed main(){
ios::sync_with_stdio(false);
cin>>n;
for(ll i=0;i<n;i++){
cin>>pic[i].x1>>pic[i].y1>>pic[i].x2>>pic[i].y2>>pic[i].c;
vec[pic[i].x2].push_back(i);
}
for(ll i=0;i<n;i++){
ll L=pic[i].y1,R=pic[i].y2;
for(ll it:vec[pic[i].x1]){
ll x=pic[it].y1,y=pic[it].y2;
if(!(y<=L||x>=R)){
eg[i].push_back(it);
}
}
}
memset(dp,0x3f,sizeof(dp));
for(ll i=1;i<(1ll<<n);i++){
if(!check(i))continue;
for(ll j=0;j<n;j++){
if(!((1ll<<j)&i))continue;
if(lowbit(i^lowbit(i))==0)dp[i][pic[j].c]=1;
for(ll k=1;k<=20;k++){
if(pic[j].c==k)dp[i][pic[j].c]=min(dp[i][pic[j].c],dp[i-(1ll<<j)][k]);
else dp[i][pic[j].c]=min(dp[i][pic[j].c],dp[i-(1ll<<j)][k]+1);
}
}
}
ll ans=INF;
for(ll i=0;i<n;i++)
ans=min(ans,dp[(1ll<<n)-1][i]);
cout<<ans;
return 0;
}