Google code jam 2008 R1A - Milkshakes

柔情只为你懂 2021-12-13 03:49 220阅读 0赞

这是Google code jam 2008 R1A的第二道题，关于如何最大限度顾客需求的问题。这可以归类为满足性问题，它是一个NP-完全问题。然而这道题有自己的特征，可以根据题目给定的约束条件，采用比较简单的算法就能得出结果。最后给出源码。  
  
**Problem**

You own a milkshake shop. There are N different flavors that you can prepare, and each flavor can be prepared "malted" or "unmalted". So, you can make 2N different types of milkshakes.

Each of your customers has a set of milkshake types that they like, and they will be satisfied if you have at least one of those types prepared. At most one of the types a customer likes will be a "malted" flavor.

You want to make N batches of milkshakes, so that:

*  There is exactly one batch for each flavor of milkshake, and it is either malted or unmalted.
 *  For each customer, you make at least one milkshake type that they like.
 *  The minimum possible number of batches are malted.

Find whether it is possible to satisfy all your customers given these constraints, and if it is, what milkshake types you should make.

If it is possible to satisfy all your customers, there will be only one answer which minimizes the number of malted batches.

**Input**

*  One line containing an integer **C**, the number of test cases in the input file.

For each test case, there will be:

*  One line containing the integer **N**, the number of milkshake flavors.
 *  One line containing the integer **M**, the number of customers.
 *  **M** lines, one for each customer, each containing:
    
     *  An integer **T** >= 1, the number of milkshake types the customer likes, followed by
     *  **T** pairs of integers "**X Y**", one for each type the customer likes, where **X** is the milkshake flavor between **1** and **N** inclusive, and **Y** is either 0 to indicate unmalted, or 1 to indicated malted. Note that:
        
         *  No pair will occur more than once for a single customer.
         *  Each customer will have at least one flavor that they like (T >= 1).
         *  Each customer will like at most one malted flavor. (At most one pair for each customer has Y = 1).
    
    All of these numbers are separated by single spaces.

**Output**

*  **C** lines, one for each test case in the order they occur in the input file, each containing the string "Case \#**X**: " where **X** is the number of the test case, starting from 1, followed by:
    
     *  The string "**IMPOSSIBLE**", if the customers' preferences cannot be satisfied; **OR**
     *  **N** space-separated integers, one for each flavor from **1** to **N**, which are 0 if the corresponding flavor should be prepared unmalted, and 1 if it should be malted.

**Limits**

Small dataset

C = 100 
1 <= N <= 10 
1 <= M <= 100

Large dataset

C = 5 
1 <= N <= 2000 
1 <= M <= 2000

The sum of all the **T** values for the customers in a test case will not exceed 3000.

**Sample**

<table> 
 <tbody> 
 <tr> 
 <td> Input &nbsp; </td> 
 <td> Output &nbsp; </td> 
 </tr> 
 <tr> 
 <td><code>2 5 3 1 1 1 2 1 0 2 0 1 5 0 1 2 1 1 0 1 1 1 </code></td> 
 <td><code>Case #1: 1 0 0 0 0 Case #2: IMPOSSIBLE </code></td> 
 </tr> 
 </tbody> 
</table>

In the first case, you must make flavor \#1 malted, to satisfy the first customer. Every other flavor can be unmalted. The second customer is satisfied by getting flavor \#2 unmalted, and the third customer is satisfied by getting flavor \#5 unmalted.

In the second case, there is only one flavor. One of your customers wants it malted and one wants it unmalted. You cannot satisfy them both.

**Analysis**  
On the surface, this problem appears to require solving the classic problem "Satisfiability," the canonical example of an NP-complete problem. The customers represent clauses, the milkshake flavors represent variables, and malted and unmalted flavors represent whether the variable is negated.

We are not evil enough to have chosen a problem that hard! The restriction that makes this problem easier is that the customers can only like at most one malted flavor (or equivalently, the clauses can only have at most one negated variable.)

Using the following steps, we can quickly find whether a solution exists, and if so, what the solution is.

1.  Start with every flavor unmalted and consider the customers one by one.
2.  If there is an unsatisfied customer who only likes unmalted flavors, and all those flavors have been made malted, then no solution is possible.
3.  If there is an unsatisfied customer who has one favorite malted flavor, then we must make that flavor malted. We do this, then go back to step 2.
4.  If there are no unsatisfied customers, then we already have a valid solution and can leave the remaining flavors unmalted.

Notice that whenever we made a flavor malted, we were forced to do so. Therefore, the solution we got must have the minimum possible number of malted flavors.

With clever data structures, the above algorithm can be implemented to run in linear time.

More information:

[The Satisfiability problem][] \- [Horn clauses][]

**Source Code**

![None.gif][] \#include  < iostream > 
![None.gif][] 
![None.gif][] using   namespace  std; 
![None.gif][] 
![None.gif][] \#define  Rep(i,n) for (int i(0),\_n(n); i<\_n; ++i) 
![None.gif][] 
![ExpandedBlockStart.gif][]![ContractedBlock.gif][] struct  Flavor ![dot.gif][] \{ 
![InBlock.gif][]    int X; 
![InBlock.gif][]    char Y; 
![ExpandedBlockEnd.gif][]\} ; 
![None.gif][] 
![ExpandedBlockStart.gif][]![ContractedBlock.gif][] struct  Customer ![dot.gif][] \{ 
![InBlock.gif][]    int T; 
![InBlock.gif][]    Flavor\* F; 
![ExpandedSubBlockStart.gif][]![ContractedSubBlock.gif][]    Customer() ![dot.gif][]\{ 
![InBlock.gif][]        F = NULL; 
![ExpandedSubBlockEnd.gif][]    \} 
![ExpandedSubBlockStart.gif][]![ContractedSubBlock.gif][]    ~Customer() ![dot.gif][]\{ 
![ExpandedSubBlockStart.gif][]![ContractedSubBlock.gif][]        if(NULL!=F) ![dot.gif][]\{ 
![InBlock.gif][]            delete\[\] F; 
![InBlock.gif][]            F = NULL; 
![ExpandedSubBlockEnd.gif][]        \} 
![ExpandedSubBlockEnd.gif][]    \} 
![ExpandedSubBlockStart.gif][]![ContractedSubBlock.gif][]    void Init(int t) ![dot.gif][]\{ 
![InBlock.gif][]        T = t; 
![InBlock.gif][]        F = new Flavor\[T\]; 
![ExpandedSubBlockEnd.gif][]    \} 
![ExpandedSubBlockStart.gif][]![ContractedSubBlock.gif][]    void SetFlavor(int i, int X, int Y) ![dot.gif][]\{ 
![InBlock.gif][]        F\[i\].X = X; 
![InBlock.gif][]        F\[i\].Y = Y; 
![ExpandedSubBlockEnd.gif][]    \} 
![ExpandedSubBlockStart.gif][]![ContractedSubBlock.gif][]    int GetFlavorX(int i) ![dot.gif][]\{ 
![InBlock.gif][]        return F\[i\].X; 
![ExpandedSubBlockEnd.gif][]    \} 
![ExpandedSubBlockStart.gif][]![ContractedSubBlock.gif][]    int GetFlavorY(int i) ![dot.gif][]\{ 
![InBlock.gif][]        return F\[i\].Y; 
![ExpandedSubBlockEnd.gif][]    \} 
![ExpandedSubBlockStart.gif][]![ContractedSubBlock.gif][]    bool IsSatisfied() ![dot.gif][]\{ 
![InBlock.gif][]        return T==0; 
![ExpandedSubBlockEnd.gif][]    \} 
![ExpandedSubBlockStart.gif][]![ContractedSubBlock.gif][]    void Satisfy() ![dot.gif][]\{ 
![InBlock.gif][]        T=0; 
![ExpandedSubBlockStart.gif][]![ContractedSubBlock.gif][]        if(NULL!=F) ![dot.gif][]\{ 
![InBlock.gif][]            delete\[\] F; 
![InBlock.gif][]            F = NULL; 
![ExpandedSubBlockEnd.gif][]        \} 
![ExpandedSubBlockEnd.gif][]    \} 
![ExpandedSubBlockStart.gif][]![ContractedSubBlock.gif][]    bool IsSatisfing(int i, int \*f) ![dot.gif][]\{ 
![InBlock.gif][]        return f\[F\[i\].X\]==F\[i\].Y; 
![ExpandedSubBlockEnd.gif][]    \} 
![ExpandedSubBlockStart.gif][]![ContractedSubBlock.gif][]    void SetMalted(int i, int \*f) ![dot.gif][]\{ 
![InBlock.gif][]        f\[F\[i\].X\] = 1; 
![ExpandedSubBlockEnd.gif][]    \} 
![ExpandedBlockEnd.gif][]\} ; 
![None.gif][] 
![None.gif][] int  main() 
![ExpandedBlockStart.gif][]![ContractedBlock.gif][] ![dot.gif][] \{ 
![InBlock.gif][]    int C; 
![InBlock.gif][]    FILE \*fp = fopen("A.out", "w"); 
![InBlock.gif][]    scanf("%d", &C); 
![ExpandedSubBlockStart.gif][]![ContractedSubBlock.gif][]    Rep(c, C) ![dot.gif][]\{ 
![InBlock.gif][]        int N; 
![InBlock.gif][]        scanf("%d", &N); 
![InBlock.gif][]        int\* f = new int\[N\+1\]; 
![ExpandedSubBlockStart.gif][]![ContractedSubBlock.gif][]        Rep(i ,N\+1) ![dot.gif][]\{ 
![InBlock.gif][]            f\[i\]=0; 
![ExpandedSubBlockEnd.gif][]        \} 
![InBlock.gif][]        int M; 
![InBlock.gif][]        scanf("%d", &M); 
![InBlock.gif][]        Customer\* customer = new Customer\[M\]; 
![ExpandedSubBlockStart.gif][]![ContractedSubBlock.gif][]        Rep(m, M) ![dot.gif][]\{ 
![InBlock.gif][]            int T; 
![InBlock.gif][]            scanf("%d", &T); 
![InBlock.gif][]            customer\[m\].Init(T); 
![ExpandedSubBlockStart.gif][]![ContractedSubBlock.gif][]            Rep(t, T) ![dot.gif][]\{ 
![InBlock.gif][]                int X, Y; 
![InBlock.gif][]                scanf("%d%d", &X,&Y); 
![InBlock.gif][]                customer\[m\].SetFlavor(t,X,Y); 
![ExpandedSubBlockEnd.gif][]            \} 
![ExpandedSubBlockEnd.gif][]        \} 
![InBlock.gif][]        bool findSolution = true; 
![InBlock.gif][]        int m = 0; 
![ExpandedSubBlockStart.gif][]![ContractedSubBlock.gif][]        while(m<M) ![dot.gif][]\{ 
![ExpandedSubBlockStart.gif][]![ContractedSubBlock.gif][]            if(customer\[m\].IsSatisfied()) ![dot.gif][]\{ 
![InBlock.gif][]                m\++; 
![InBlock.gif][]                continue; 
![ExpandedSubBlockEnd.gif][]            \} 
![InBlock.gif][]            bool malted = false; 
![InBlock.gif][]            int idx; 
![InBlock.gif][]            bool satisfied = false; 
![ExpandedSubBlockStart.gif][]![ContractedSubBlock.gif][]            Rep(t, customer\[m\].T) ![dot.gif][]\{ 
![ExpandedSubBlockStart.gif][]![ContractedSubBlock.gif][]                if(customer\[m\].GetFlavorY(t)==1) ![dot.gif][]\{ 
![InBlock.gif][]                    malted = true; 
![InBlock.gif][]                    idx = t; 
![ExpandedSubBlockEnd.gif][]                \} 
![ExpandedSubBlockStart.gif][]![ContractedSubBlock.gif][]                if(customer\[m\].IsSatisfing(t,f)) ![dot.gif][]\{ 
![InBlock.gif][]                    satisfied = true; 
![ExpandedSubBlockEnd.gif][]                \} 
![ExpandedSubBlockEnd.gif][]            \} 
![ExpandedSubBlockStart.gif][]![ContractedSubBlock.gif][]            if(!satisfied) ![dot.gif][]\{ 
![ExpandedSubBlockStart.gif][]![ContractedSubBlock.gif][]                if(malted) ![dot.gif][]\{ 
![InBlock.gif][]                    customer\[m\].SetMalted(idx,f); 
![InBlock.gif][]                    customer\[m\].Satisfy(); 
![InBlock.gif][]                    m=0; 
![ExpandedSubBlockStart.gif][]![ContractedSubBlock.gif][]                \} else ![dot.gif][]\{ 
![InBlock.gif][]                    findSolution = false; 
![InBlock.gif][]                    break; 
![ExpandedSubBlockEnd.gif][]                \} 
![ExpandedSubBlockStart.gif][]![ContractedSubBlock.gif][]            \} else ![dot.gif][]\{ 
![InBlock.gif][]                m\++; 
![ExpandedSubBlockEnd.gif][]            \} 
![ExpandedSubBlockEnd.gif][]        \} 
![InBlock.gif][]        fprintf(fp,"Case \#%d: ", c\+1); 
![ExpandedSubBlockStart.gif][]![ContractedSubBlock.gif][]        if(findSolution) ![dot.gif][]\{ 
![ExpandedSubBlockStart.gif][]![ContractedSubBlock.gif][]            Rep(i ,N) ![dot.gif][]\{ 
![InBlock.gif][]                fprintf(fp,"%d ", f\[i\+1\]); 
![ExpandedSubBlockEnd.gif][]            \} 
![InBlock.gif][]            fprintf(fp,"\\n"); 
![ExpandedSubBlockStart.gif][]![ContractedSubBlock.gif][]        \} else ![dot.gif][]\{ 
![InBlock.gif][]            fprintf(fp,"IMPOSSIBLE\\n"); 
![ExpandedSubBlockEnd.gif][]        \} 
![InBlock.gif][]        delete\[\] customer; 
![InBlock.gif][]        delete\[\] f; 
![InBlock.gif][] 
![ExpandedSubBlockEnd.gif][]    \} 
![InBlock.gif][]    fclose(fp); 
![ExpandedBlockEnd.gif][]\} 
![None.gif][]

转载于:https://www.cnblogs.com/clive/archive/2009/08/13/Google\_code\_jam\_2008\_R1A\_Milkshakes.html

[The Satisfiability problem]: http://www.google.com/search?q=Satisfiability+problem
[Horn clauses]: http://www.google.com/search?q=Horn+clauses
[None.gif]: /images/20211210/6b128536f4a642dd8a55e1b27978e441.png
[ExpandedBlockStart.gif]: http://www.cppblog.com/Images/OutliningIndicators/ExpandedBlockStart.gif
[ContractedBlock.gif]: http://www.cppblog.com/Images/OutliningIndicators/ContractedBlock.gif
[dot.gif]: /images/20211210/0102f04c49424c3e97c42a6ebe3d3ba5.png
[InBlock.gif]: /images/20211210/7ed216420ffd46a7bfe866bdfb1b8738.png
[ExpandedBlockEnd.gif]: /images/20211210/6161ac02f821455087d6b847a7046ea0.png
[ExpandedSubBlockStart.gif]: http://www.cppblog.com/Images/OutliningIndicators/ExpandedSubBlockStart.gif
[ContractedSubBlock.gif]: http://www.cppblog.com/Images/OutliningIndicators/ContractedSubBlock.gif
[ExpandedSubBlockEnd.gif]: /images/20211210/2db8785f68294c8aa07c514ed8199ea1.png