确实如y总所说
状态压缩dp真不难
其可以一眼看出的范围
和必须暴力枚举的思路也不是很难想
你问我为啥今天才三题
不知道 感觉假期中摆烂还是很正常吧!
明天想出去玩了 咕咕咕
我这样真的能拿au吗 估计挺难的吧 哈哈
AcWing 292. 炮兵阵地
说hard 也不是很 hard
感觉和昨天一样的套路
本来以为四层循环会爆掉
好吧 并没有 而且我这里用的 g[i-2]都没有越界 很神秘
还学到了滚动数组
很厉害呢 先照常写 写过了再加个& 就可以了呢!
#include <bits/stdc++.h>
using namespace std;
const int N = 110;
const int M = 1<<10;
const int mod = 1e8;
#define int long long
#define endl '\n'
#define Endl '\n'
#define _ 0
#define inf 0x3f3f3f3f3f3f3f3f
#define fast ios::sync_with_stdio(false);cin.tie(nullptr);
int max(int x,int y){return x>y?x:y;}
int min(int x,int y){return x<y?x:y;}
int g[N],f[2][M][M],n,m,cnt[M];
vector<int>state;
bool check(int x){
for(int i=0;i<m;i++){
if((x>>i&1)&&((x>>(i+1)&1)||(x>>(i+2)&1)))return false;
}
return true;
}
int count(int x){
int res=0;
for(int i=0;i<m;i++)if(x>>i&1)res++;
return res;
}
signed main(){
fast
cin>>n>>m;
for(int i=1;i<=n;i++){
for(int j=1;j<=m;j++){
char x;cin>>x;
if(x=='H')g[i]+=1<<(j-1);
}
}
for(int i=0;i<1<<m;i++){
if(check(i)){
state.push_back(i);
cnt[i]=count(i);
}
}
for(int i=1;i<=n;i++){
for(auto a:state){
for(auto b:state){
for(auto c:state){//j i k i-1
if(i>=2&&(a&b)==0&&(a&c)==0&&(b&c)==0&&(g[i]&a)==0&&(g[i-1]&b)==0&&(g[i-2]&c)==0)
f[i&1][a][b]=max(f[i&1][a][b],f[i-1&1][b][c]+cnt[a]);
}
}
}
}
int ans=0;
for(int i=0;i<1<<m;i++)
for(int j=0;j<1<<m;j++)
ans=max(ans,f[n&1][i][j]);
cout<<ans<<endl;
return ~~(0^_^0);
}
524. 愤怒的小鸟
卧槽 怎么还有卡精度的题啊
说实话这道题 so hard
思路:
我们知道有n个点 那么我们可以搞出n2条抛物线(注意有些抛物线可以只过一个点
然后就类似与重复覆盖问题
一般我们都是搜索
不过今天用了状压优化 相当于
因为我们这个与顺序无关 所以我们可以根据每个状态不同 考虑刷表 枚举一个没有被打死的猪即可
思路挺简单的 可是就是细节有点多
path存的是过了i j 两点的状态 卡精度 a要<0 枚举两个点时 不能在一条直线上
#include <bits/stdc++.h>
using namespace std;
const int N = 20;
const int M = 1<<18;
const int mod = 1e8;
#define int long long
#define endl '\n'
#define Endl '\n'
#define _ 0
#define x first
#define y second
#define inf 0x3f3f3f3f3f3f3f3f
#define fast ios::sync_with_stdio(false);cin.tie(nullptr);
int max(int x,int y){return x>y?x:y;}
int min(int x,int y){return x<y?x:y;}
int f[M],path[N][N],n,m;//path[]存的是过i j点线的状态
pair<double,double>q[N];
int cmp(double x, double y)
{
if (fabs(x - y) < 1e-6) return 0;
if (x < y) return -1;
return 1;
}
signed main(){
fast
int t;cin>>t;
while(t--){
cin>>n>>m;
for(int i=0;i<n;i++)cin>>q[i].x>>q[i].y;
memset(path,0,sizeof path);
for(int i=0;i<n;i++){
path[i][i]=1<<i;
for(int j=0;j<n;j++){
double x1=q[i].x,y1=q[i].y,x2=q[j].x,y2=q[j].y;
double a=((y1/x1)-(y2/x2))/(x1-x2);
if(!cmp(x1,x2))continue;
if (cmp(a,0) >= 0) continue;
double b=(y1/x1)-a*x1;
int state=0;
for(int k=0;k<n;k++){
double x=q[k].x,y=q[k].y;
if(!cmp(a*x*x+b*x,y))state+=1<<k;
}
path[i][j]=state;
}
}
memset(f,0x3f3f,sizeof f);
f[0]=0;
for(int i=0;i+1<1<<n;i++){
int x=0;
for(int j=0;j<n;j++){
if((i>>j&1)==0){
x=j;
break;
}
}
for(int j=0;j<n;j++){
f[i|path[x][j]]=min(f[i|path[x][j]],f[i]+1);
}
}
cout<<f[(1<<n)-1]<<Endl;
}
return ~~(0^_^0);
}
P3959 [NOIP2017 提高组] 宝藏
hard
我们都知道暴力怎么打
每个点作为起点爆搜一次
所以状压dp也就自然而然出来了
我们因为有这个条件
新开发一条道路的代价是 \mathrm{L} \times \mathrm{K}L×K。其中 LL 代表这条道路的长度,KK 代表从赞助商帮你打通的宝藏屋到这条道路起点的宝藏屋所经过的宝藏屋的数量(包括赞助商帮你打通的宝藏屋和这条道路起点的宝藏屋)
所以要维护一维dis 自然还要维护一维点 最后再维护一维状态
状态转移自然:
f[j][1<<(j-1)|s][depth+1]=sum+dis[i]*g[i][j];
#include <bits/stdc++.h>
using namespace std;
const int N = 15;
const int M = 1<<12;
const int mod = 1e8;
#define int long long
#define endl '\n'
#define Endl '\n'
#define _ 0
#define inf 0x3f3f3f3f3f3f3f3f
#define fast ios::sync_with_stdio(false);cin.tie(nullptr);
int max(int x,int y){return x>y?x:y;}
int min(int x,int y){return x<y?x:y;}
int n,m,ans=inf,g[N][N],dis[N],f[N][M][N];
void dfs(int s,int sum,int depth){
if(sum>=ans)return;
if(s==(1<<n)-1){ans=sum;return;}
for(int i=1;i<=n;i++){
if(!(1<<(i-1)&s))continue;
for(int j=1;j<=n;j++){
if(!(1<<(j-1)&s)&&g[i][j]<2e9){
if(f[j][1<<(j-1)|s][depth+1]<=sum+dis[i]*g[i][j])continue;
f[j][1<<(j-1)|s][depth+1]=sum+dis[i]*g[i][j];
dis[j]=dis[i]+1;
dfs(1<<(j-1)|s,f[j][1<<(j-1)|s][depth+1],depth+1);
}
}
}
}
signed main(){
fast
cin>>n>>m;
memset(g,0x3f3f,sizeof g);
for(int i=1;i<=m;i++){
int x,y,w;cin>>x>>y>>w;
g[x][y]=g[y][x]=min(w,g[x][y]);
}
for(int i=1;i<=n;i++){
memset(dis,0,sizeof dis);
memset(f,0x3f3f,sizeof f);
dis[i]=1;dfs(1<<(i-1),0,0);
}
cout<<ans<<endl;
return ~~(0^_^0);
}
话说这还是一道ac率不高的紫题欸