jzoj4016-圈地为王【状压,bfs,几何】
生活随笔
收集整理的這篇文章主要介紹了
jzoj4016-圈地为王【状压,bfs,几何】
小編覺得挺不錯(cuò)的,現(xiàn)在分享給大家,幫大家做個(gè)參考.
正題
題目鏈接:https://jzoj.net/senior/#contest/show/3011/1
題目大意
n?mn*mn?m的格子,格子之間有道路,對于每個(gè)iii就走過最短的回路使得
解題思路
判斷一個(gè)點(diǎn)是否在一個(gè)多邊形內(nèi),我們可以往任何一個(gè)方向畫一條射線,如果與多邊形的交點(diǎn)為奇數(shù)那么就在,否則就不在。
那么我們考慮狀態(tài)壓縮fi,j,sf_{i,j,s}fi,j,s?表示走到(i,j)(i,j)(i,j)這個(gè)位置,往上畫線的交點(diǎn)為奇數(shù)的有用點(diǎn)格子集合為sss。
然后bfsbfsbfs轉(zhuǎn)移即可。
codecodecode
#include<cstdio> #include<cstring> #include<algorithm> #include<queue> #define lowbit(x) (x&-x) using namespace std; const int dx[4]={1,0,-1,0},dy[4]={0,1,0,-1}; struct point{int x,y,s; }val[50],no[50]; queue<point> q; int n,m,cnt1,cnt2,ans[50],f[55][55][1500]; char c[55]; int count_one(int x) {int ans=0;while(x){ans++;x-=lowbit(x);}return ans; } int main() {freopen("data.txt","r",stdin);cnt1=cnt2=-1;while(scanf("%s",c+1)!=EOF){n++;m=strlen(c+1);for(int j=1;j<=m;j++){if(c[j]=='I')val[++cnt1]=(point){n,j};if(c[j]=='X')no[++cnt2]=(point){n,j};}}memset(f,0x3f,sizeof(f));f[0][0][0]=0;q.push((point){0,0,0});while(!q.empty()){int x=q.front().x,y=q.front().y,s=q.front().s;q.pop();for(int k=0;k<4;k++){int zx=x+dx[k],zy=y+dy[k];if(zx<0||zx>n||zy<0||zy>m) continue;int tmp=s;if(k==1){for(int j=0;j<=cnt1;j++)if(val[j].y==y+1&&val[j].x>=x+1) tmp^=(1<<j); for(int j=0;j<=cnt2;j++)if(no[j].y==y+1&&no[j].x>=x+1) tmp^=(1<<j+cnt1+1); }if(k==3){for(int j=0;j<=cnt1;j++)if(val[j].y==y&&val[j].x>=x+1) tmp^=(1<<j); for(int j=0;j<=cnt2;j++)if(no[j].y==y&&no[j].x>=x+1) tmp^=(1<<j+cnt1+1); }if(f[x][y][s]+1<f[zx][zy][tmp]){q.push((point){zx,zy,tmp});f[zx][zy][tmp]=f[x][y][s]+1;}}}memset(ans,0x3f,sizeof(ans));cnt1++;for(int i=0;i<(1<<cnt1);i++){int d=count_one(i);ans[d]=min(ans[d],f[0][0][i]);}for(int i=1;i<=cnt1;i++)printf("%d\n",ans[i]); } 創(chuàng)作挑戰(zhàn)賽新人創(chuàng)作獎(jiǎng)勵(lì)來咯,堅(jiān)持創(chuàng)作打卡瓜分現(xiàn)金大獎(jiǎng)總結(jié)
以上是生活随笔為你收集整理的jzoj4016-圈地为王【状压,bfs,几何】的全部內(nèi)容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: jzoj4015-数列【循环节,数论】
- 下一篇: jzoj4017-逃跑【0/1分数规划,