Ural 1519. Formula 1 优美的插头DP

 

分类： ACM_Ural ACM_插头DP 2013-08-03 20:45  333人阅读  评论(0)  收藏  举报

今天早上学了插头DP的思想和最基础的应用，中午就开始敲了，岐哥说第一次写不要看别人代码，利用自己的理解一点点得写出来，这样才锻炼代码能力！于是下午慢慢地构思轮廓，一点点地敲出主体代码，其实是很磨蹭的，由于要考虑好多东西，而且昨晚2点睡的有点困，最后终于磨蹭出来了，第一次的代码搓没关系，自己写的才重要。然后果然不出我所料，调试到了晚上才A了（一个郁闷的错误）。。。A的感觉真的是爽呀，虽然搞了差不多一天。当然自己写了自己想的代码后也要把代码优化，不然队友看不懂自己代码就囧了。。。

插头DP，建议大家想学的好好看看陈丹琦的国家集训队论文，这是个优美的DP。

http://www.docin.com/p-46797997.html

 #include <stdio.h>
 #include <string.h>

 #define LL __int64

 const int mod = 10007;

 //   哈希表
 struct HASH{
     int head[mod+10], E, next[80000];
     LL val[80000], cnt[80000];

     void init() {
         memset(head, -1, sizeof(head));
         E = 0;
     }

     int findhash(LL x) {
         return (x%mod + mod)%mod;
     }

     void add(LL x, LL sum) {
         int u = findhash(x);
         for(int i = head[u];i != -1;i = next[i]) if(val[i] == x) {
             cnt[i] += sum;
             return ;
         }
         val[E] = x;
         cnt[E] = sum;
         next[E] = head[u];
         head[u] = E++;
     }

 }biao1, biao2;

 int c[22], n, m, d[22];
 //  编码
 void get(LL x) {
     for(int i = m+1;i >= 1; i--) {
         c[i] = x&7;
         x /= 8;
     }
 }
 //  解码
 LL getval() {
     LL ret = 0;
     for(int i = 1;i <= m+1; i++) {
         ret |= d[i];
         ret *= 8;
     }
     ret /= 8;
     return ret ;
 }
 //   转化成最小表示法
 void change() {
     int vis[22];
     memset(vis, 0, sizeof(vis));
     int num = 1;
     for(int i = 1;i <= m+1;i ++) {
         if(!d[i])   continue;
         if(!vis[d[i]]) {
             vis[d[i]] = num;
             d[i] = num++;
         }
         else {
             d[i] = vis[d[i]];
         }
     }
 }

 void fuzhi() {
     for(int i = 1;i <= m+1;i ++)    d[i] = c[i];
 }

 char s[22][22];

 int main() {
     int i, j, k, l;
     while(scanf("%d%d", &n, &m) != -1) {
         for(i = 1;i <= n; i++)
             scanf("%s", s[i]+1);
         int tot = 0;
         for(i = 1;i <= n; i++)
             for(j = 1;j <= m; j++)
             if(s[i][j] == '.')  tot++;
         if(tot%2==1 || tot < 4) {
             puts("0");
             continue;
         }
         int tox  = -1, toy = -1;
         for(i = 1;i <= n; i++)
         for(j = 1;j <= m; j++) if(s[i][j] == '.') {
             tox = i;
             toy = j;
         }
         biao1.init();
         biao1.add(0, 1);
         LL ans = 0;
         for(i = 1;i <= n; i++){
             for(j = 0;j <= m; j++){
                 biao2.init();
                 for(l = 0;l < biao1.E; l++)  {
                     get(biao1.val[l]);
                     if(j == m) {
                         for(int ii = 2;ii <= m+1; ii++) d[ii] = c[ii-1];
                         d[1] = 0;
                         change();
                         LL now = getval();
                         biao2.add(now, biao1.cnt[l]);
                         continue;
                     }
                     if(c[j+1] && !c[j+2]) {  //  有左插头无上插头
                         if(s[i][j+1] != '.')  continue;
                         if(j+2 <= m) {
                             fuzhi();
                             d[j+1] = 0;d[j+2] = c[j+1];
                             change();
                             LL now = getval();
                             biao2.add(now, biao1.cnt[l]);
                         }
                         if(i < n) {
                             fuzhi();
                             change();
                             LL now = getval();
                             biao2.add(now, biao1.cnt[l]);
                         }
                     }
                     else if(!c[j+1] && c[j+2]) {  //  有上插头无左插头
                         if(s[i][j+1] != '.')  continue;
                         if(i < n) {
                             fuzhi();
                             d[j+1] = c[j+2]; d[j+2] = 0;
                             change();
                             LL now = getval();
                             biao2.add(now, biao1.cnt[l]);
                         }
                         if(j+2 <= m) {
                             fuzhi();
                             change();
                             LL now = getval();
                             biao2.add(now, biao1.cnt[l]);
                         }
                     }
                     else if(!c[j+1] && !c[j+2]) { //   左和上都无插头
                         if(s[i][j+1] != '.') {
                             fuzhi();
                             change();
                             LL now = getval();
                             biao2.add(now, biao1.cnt[l]);
                             continue;
                         }
                         if(j+2 <= m && i < n) {
                             fuzhi();
                             d[j+1] = d[j+2] = 13;
                             change();
                             LL now = getval();
                             biao2.add(now, biao1.cnt[l]);
                         }
                     }
                     else { //　　左和上都有插头　，　　要判断左和上插头是否连通
                         if(c[j+2] == c[j+1]) {
                             int tot = 0;
                             for(int ii = 1;ii <= m+1; ii++) if(c[ii])
                                 tot++;
                             if(tot == 2 && i == tox && j+1 == toy) ans += biao1.cnt[l];
                         }
                         else {
                             if(s[i][j+1] != '.')    continue;
                             fuzhi();
                             for(int ii = 1;ii <= m+1; ii++) if(ii != j+1 && ii != j+2 && d[ii] == d[j+1]) {
                                 d[ii] = d[j+2];
                                 break;
                             }
                             d[j+1] = d[j+2] = 0;
                             change();
                             LL now = getval();
                             biao2.add(now, biao1.cnt[l]);
                         }
                     }
                 }
                 biao1 = biao2;
             }
         }
         printf("%I64d\n", ans);
     }
     return 0;
 }