码迷,mamicode.com
首页 > 其他好文 > 详细

【BZOJ】【1021】【SHOI2008】Dept循环的债务

时间:2015-04-12 22:43:51      阅读:277      评论:0      收藏:0      [点我收藏+]

标签:

DP


  去膜拜题解了>_>玛雅原来是动规……

  让我先理解一下为什么要用动规:这个题根据钱数推方案其实是无从下手的……(线性规划?……事实证明我想多了)

  啦~我们先来看个超级简化版的问题:怎么判无法还清?正着判很麻烦对不对= =(其实是我没想……)

  那么我们倒着来考虑:有哪些状态是我们通过交换钱币能够到达的,这个可以递推对不>_>

  现在我们就知道哪些状态我们是可以到达的了……再多想一下……递推……如果我们依次考虑每种面额的交换策略,顺便也就知道了我们到达这个状态的最小交换次数对吧?

  原因:因为每种面值只有6种策略(给出钱的是一个人,3种,给出钱的是两个人,3种),且不同面值之间的交换互不影响(这不是废话吗!我手里有多少张10块的难道还会影响我手里100块的张数?可是我一开始绕进这个弯里没走出来TAT)

  (P.S.看到这里你就可以去写题了,因为我动规方程比较奇葩……自己想的时候又把自己坑了

  我是让$f[i][j][k]$表示当前a手里有 j 块钱,b手里有 k 块钱,现在该考虑第 i 种面额的钞票(但是还没有对第 i 种钞票进行分配!!我忘了啊啊啊!!!看了题解以后跟人家的方法搞混了?>_>)

  那么我是分3种情况讨论的(将上面的6种情况组合了一下)……1. a给b、c或是b、c给a ; 2.b给a、c或是a、c给b;3.c给a、b或是a、b给c。

  明显可行的状态比数组大小小得多,所以就用了队列来保存/扩展可行状态(类似记忆化搜索的策略吧)

  (不过我这种姿势明显是比较捉急的……不过也能过>_>想看 快/短 一点的做法请戳@regina8023 @rausen

  自我感觉我写这题的动规的姿势还是蛮特别的,而且还是第666个提交 2333~沾沾自喜一下(>_<)~

技术分享
  1 /**************************************************************
  2     Problem: 1021
  3     User: Tunix
  4     Language: C++
  5     Result: Accepted
  6     Time:784 ms
  7     Memory:52136 kb
  8 ****************************************************************/
  9  
 10 //BZOJ 1021
 11 #include<vector>
 12 #include<cstdio>
 13 #include<cstring>
 14 #include<cstdlib>
 15 #include<iostream>
 16 #include<algorithm>
 17 #define rep(i,n) for(int i=0;i<n;++i)
 18 #define F(i,j,n) for(int i=j;i<=n;++i)
 19 #define D(i,j,n) for(int i=j;i>=n;--i)
 20 #define pb push_back
 21 using namespace std;
 22 inline int getint(){
 23     int v=0,sign=1; char ch=getchar();
 24     while(ch<0||ch>9){ if (ch==-) sign=-1; ch=getchar();}
 25     while(ch>=0&&ch<=9){ v=v*10+ch-0; ch=getchar();}
 26     return v*sign;
 27 }
 28 const int M=1e6+10,INF=~0u>>2;
 29 typedef long long LL;
 30 /******************tamplate*********************/
 31 const int money[8]={0,1,5,10,20,50,100};
 32 int a[8],b[8],c[8],sa,sb,sc,ea,eb,ec;
 33 int n,m,f[8][1001][1001],tot;
 34 bool inq[8][1001][1001];
 35 struct node{
 36     int x,y,z;
 37 }Q[M];
 38 void dp(){
 39     F(i,1,7) F(j,0,tot) F(k,0,tot) f[i][j][k]=INF;
 40     int l=0,r=-1;
 41     Q[++r]=(node){1,sa,sb}; f[1][sa][sb]=0; inq[1][sa][sb]=1;
 42     while(l<=r){
 43         node now=Q[l++];
 44         int x=now.x,y=now.y,z=now.z;
 45         inq[x][y][z]=0;
 46         if (x==7) continue;
 47         if (f[x+1][y][z]>f[x][y][z]){
 48             f[x+1][y][z]=f[x][y][z];
 49             if (!inq[x+1][y][z]){
 50                 Q[++r]=(node){x+1,y,z};
 51                 inq[x+1][y][z]=1;
 52             }
 53         }
 54         F(i,-b[x],a[x]) F(j,-c[x],a[x]-i){
 55             int ty=y-(i+j)*money[x],tz=z+i*money[x];
 56             if (ty<0) break;
 57             if (f[x+1][ty][tz]>f[x][y][z]+abs(i)+abs(j)){
 58                 f[x+1][ty][tz]=f[x][y][z]+abs(i)+abs(j);
 59                 if (!inq[x+1][ty][tz]){
 60                     Q[++r]=(node){x+1,ty,tz};
 61                     inq[x+1][ty][tz]=1;
 62                 }
 63             }
 64         }
 65         F(i,-a[x],b[x]) F(j,-c[x],b[x]-i){
 66             int ty=y+i*money[x],tz=z-(i+j)*money[x];
 67             if (tz<0) break;
 68             if (f[x+1][ty][tz]>f[x][y][z]+abs(i)+abs(j)){
 69                 f[x+1][ty][tz]=f[x][y][z]+abs(i)+abs(j);
 70                 if (!inq[x+1][ty][tz]){
 71                     Q[++r]=(node){x+1,ty,tz};
 72                     inq[x+1][ty][tz]=1;
 73                 }
 74             }
 75         }
 76         F(i,-a[x],c[x]) F(j,-b[x],c[x]-i){
 77             int ty=y+i*money[x],tz=z+j*money[x];
 78             if (ty+tz>tot) break;
 79             if (f[x+1][ty][tz]>f[x][y][z]+abs(i)+abs(j)){
 80                 f[x+1][ty][tz]=f[x][y][z]+abs(i)+abs(j);
 81                 if (!inq[x+1][ty][tz]){
 82                     Q[++r]=(node){x+1,ty,tz};
 83                     inq[x+1][ty][tz]=1;
 84                 }
 85             }
 86         }
 87     }
 88 //  F(i,1,7) F(j,0,tot) F(k,0,tot)
 89 //      if (f[i][j][k]!=INF)printf("f[%d][%d][%d]=%d\n",i,j,k,f[i][j][k]);
 90     if (f[7][ea][eb]==INF||ea<0||eb<0||ec<0||ea+eb+ec>tot) puts("impossible");
 91     else printf("%d\n",f[7][ea][eb]);
 92 }
 93  
 94  
 95 int main(){
 96 #ifndef ONLINE_JUDGE
 97     freopen("1021.in","r",stdin);
 98     freopen("1021.out","w",stdout);
 99 #endif
100     int x1=getint(),x2=getint(),x3=getint();
101     F(i,1,6) {a[7-i]=getint();sa+=a[7-i]*money[7-i];}
102     F(i,1,6) {b[7-i]=getint();sb+=b[7-i]*money[7-i];}
103     F(i,1,6) {c[7-i]=getint();sc+=c[7-i]*money[7-i];}
104     tot=sa+sb+sc;
105     ea=sa-x1+x3; eb=sb-x2+x1; ec=sc-x3+x2;
106     dp();
107     return 0;
108 }
View Code

 

 

1021: [SHOI2008]Debt 循环的债务

Time Limit: 1 Sec  Memory Limit: 162 MB
Submit: 666  Solved: 336
[Submit][Status][Discuss]

Description

Alice、 Bob和Cynthia总是为他们之间混乱的债务而烦恼,终于有一天,他们决定坐下来一起解决这个问题。不过,鉴别钞票的真伪是一件很麻烦的事情,于是他 们决定要在清还债务的时候尽可能少的交换现金。比如说,Alice欠Bob 10元,而Cynthia和他俩互不相欠。现在假设Alice只有一张50元,Bob有3张10元和10张1元,Cynthia有3张20元。一种比较直 接的做法是:Alice将50元交给Bob,而Bob将他身上的钱找给Alice,这样一共就会有14张钞票被交换。但这不是最好的做法,最好的做法 是:Alice把50块给Cynthia,Cynthia再把两张20给Alice,另一张20给Bob,而Bob把一张10块给C,此时只有5张钞票被 交换过。没过多久他们就发现这是一个很棘手的问题,于是他们找到了精通数学的你为他们解决这个难题。

Input

输 入的第一行包括三个整数:x1、x2、x3(-1,000≤x1,x2,x3≤1,000),其中 x1代表Alice欠Bob的钱(如果x1是负数,说明Bob欠了Alice的钱) x2代表Bob欠Cynthia的钱(如果x2是负数,说明Cynthia欠了Bob的钱) x3代表Cynthia欠Alice的钱(如果x3是负数,说明Alice欠了Cynthia的钱)接下来有三行,每行包括6个自然数: a100,a50,a20,a10,a5,a1 b100,b50,b20,b10,b5,b1 c100,c50,c20,c10,c5,c1 a100表示Alice拥有的100元钞票张数,b50表示Bob拥有的50元钞票张数,以此类推。另外,我们保证有 a10+a5+a1≤30,b10+b5+b1≤30,c10+c5+c1≤30,而且三人总共拥有的钞票面值总额不会超过1,000。

Output

如果债务可以还清,则输出需要交换钞票的最少张数;如果不能还清,则输出“impossible”(注意单词全部小写,输出到文件时不要加引号)。

Sample Input

输入一
10 0 0
0 1 0 0 0 0
0 0 0 3 0 10
0 0 3 0 0 0
输入二
-10 -10 -10
0 0 0 0 0 0
0 0 0 0 0 0
0 0 0 0 0 0

Sample Output

输出一
5
输出二
0

HINT

对于100%的数据,x1、x2、x3 ≤ |1,000|。

Source

[Submit][Status][Discuss]

【BZOJ】【1021】【SHOI2008】Dept循环的债务

标签:

原文地址:http://www.cnblogs.com/Tunix/p/4420832.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!