題目鏈接：https://www.nowcoder.com/acm/contest/94/A

題意：在一個二維平面上有 n 個炸彈，每個炸彈有一個坐標和爆炸半徑，引爆它之后在其半徑范圍內的炸彈也會爆炸，每個炸彈最多爆炸一次，每次隨機選一個未引爆的炸彈來引爆，問引爆所有炸彈的期望操作次數。

題解：先 dfs 把引爆每個炸彈之后會觸發的所有炸彈的狀態保存起來，當作引爆該炸彈的下一個狀態。然后可以巧妙的考慮從后往前進行狀壓 DP，顯然全部炸彈都引爆之后的操作次數為 0，此狀態則為初始狀態，接一下考慮每個狀態引爆所有可能引爆的狀態，把其轉移到的下一個狀態的期望操作次數加起來，最后除以該狀態下可以引爆的炸彈數則為該狀態的期望操作次數。

?

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 }

?

轉載于:https://www.cnblogs.com/scaulok/p/9770569.html

總結

以上是生活随笔為你收集整理的“景驰科技杯”2018年华南理工大学程序设计竞赛 A. 欧洲爆破（思维+期望+状压DP）...的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。