“景驰科技杯”2018年华南理工大学程序设计竞赛 A. 欧洲爆破(思维+期望+状压DP)...
生活随笔
收集整理的這篇文章主要介紹了
“景驰科技杯”2018年华南理工大学程序设计竞赛 A. 欧洲爆破(思维+期望+状压DP)...
小編覺得挺不錯(cuò)的,現(xiàn)在分享給大家,幫大家做個(gè)參考.
題目鏈接:https://www.nowcoder.com/acm/contest/94/A
題意:在一個(gè)二維平面上有 n 個(gè)炸彈,每個(gè)炸彈有一個(gè)坐標(biāo)和爆炸半徑,引爆它之后在其半徑范圍內(nèi)的炸彈也會爆炸,每個(gè)炸彈最多爆炸一次,每次隨機(jī)選一個(gè)未引爆的炸彈來引爆,問引爆所有炸彈的期望操作次數(shù)。
題解:先 dfs 把引爆每個(gè)炸彈之后會觸發(fā)的所有炸彈的狀態(tài)保存起來,當(dāng)作引爆該炸彈的下一個(gè)狀態(tài)。然后可以巧妙的考慮從后往前進(jìn)行狀壓 DP,顯然全部炸彈都引爆之后的操作次數(shù)為 0,此狀態(tài)則為初始狀態(tài),接一下考慮每個(gè)狀態(tài)引爆所有可能引爆的狀態(tài),把其轉(zhuǎn)移到的下一個(gè)狀態(tài)的期望操作次數(shù)加起來,最后除以該狀態(tài)下可以引爆的炸彈數(shù)則為該狀態(tài)的期望操作次數(shù)。
?
1 #include <bits/stdc++.h> 2 using namespace std; 3 #define ll long long 4 #define mst(a,b) memset((a),(b),sizeof(a)) 5 #define pi acos(-1) 6 #define pii pair<int,int> 7 const int INF = 0x3f3f3f3f; 8 const double eps = 1e-3; 9 const int MAXN = 1e5 + 10; 10 const int MAXM = 2e6 + 10; 11 const ll mod = 1e9 + 9; 12 13 int n; 14 ll x[25],y[25],r[25]; 15 ll inv[25],st[25]; 16 ll dp[1<<21]; 17 bool vis[25]; 18 19 ll pow_mod(ll a,ll m) { 20 ll res = 1; 21 while(m) { 22 if(m & 1) res = res * a % mod; 23 a = a * a % mod; 24 m /= 2; 25 } 26 return res; 27 } 28 29 ll dis(int i,int j) { 30 return (x[i]-x[j])*(x[i]-x[j]) + (y[i]-y[j])*(y[i]-y[j]); 31 } 32 33 ll num; 34 35 void dfs(int u) { 36 num |= 1<<u; 37 vis[u] = true; 38 for(int i=0; i<n; i++) { 39 if(vis[i] || dis(u,i)>r[u]*r[u]) continue; 40 dfs(i); 41 } 42 } 43 44 int main() 45 { 46 #ifdef local 47 freopen("data.txt","r",stdin); 48 // freopen("data.txt","w",stdout); 49 #endif 50 inv[0] = 1; 51 for(int i=1; i<=20; i++) 52 inv[i] = pow_mod(i,mod-2); 53 while(~scanf("%d",&n)) { 54 for(int i=0; i<n; i++) 55 scanf("%lld%lld%lld",&x[i],&y[i],&r[i]); 56 for(int i=0; i<n; i++) { 57 mst(vis,false); 58 num = 0; 59 dfs(i); 60 st[i] = num; 61 } 62 dp[(1<<n)-1] = 0; 63 for(ll i=(1<<n)-2; i>=0; i--) { 64 dp[i] = 0; 65 num = 0; 66 for(ll j=0; j<n; j++) { 67 if((1<<j)&i) { 68 num++; 69 continue; 70 } 71 dp[i] = (dp[i] + dp[i|st[j]]) % mod; 72 } 73 dp[i] = (dp[i] + n) % mod; 74 dp[i] = (dp[i] * inv[n-num]) % mod; 75 } 76 printf("%lld\n",dp[0]); 77 } 78 return 0; 79 }?
轉(zhuǎn)載于:https://www.cnblogs.com/scaulok/p/9770569.html
總結(jié)
以上是生活随笔為你收集整理的“景驰科技杯”2018年华南理工大学程序设计竞赛 A. 欧洲爆破(思维+期望+状压DP)...的全部內(nèi)容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: 血淋淋的事实告诉你:你为什么不应该在JS
- 下一篇: mysql explain的使用