[hdu1085]Holding Bin-Laden Captive!

卷积计数。

将所有面值的硬币写成多项式形式，用系数表示硬币个数，用指数表示硬币面值，则卷积一遍后系数相加，相同系数的式子累积，则再扫一遍后第一个系数为0的位置的系数便是答案。

#include <iostream>
#include <algorithm>
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <string>
#include <vector>
#include <cmath>
#include <queue>
#include <complex>
#include <set>
using namespace std;

const int N=16000;

typedef complex<double> c;

c c1[N],c2[N],c3[N],F[N];

int rev(int x,int n){
    int ret=0;
    for(int i=0;(1<<i)<n;i++)ret=(ret<<1)|((x&(1<<i))>0);
    return ret;
}

void fft(c *a,int n,int f){
    for(int i=0;i<n;i++)F[rev(i,n)]=a[i];
    for(int i=2;i<=n;i*=2){
        c wn=c(cos(2*acos(-1)*f/i),sin(2*acos(-1)*f/i));
        for(int j=0;j<n;j+=i){
            c w=1;
            for(int k=j;k<j+i/2;k++){
                c u=F[k],t=w*F[k+i/2];
                F[k]=u+t;
                F[k+i/2]=u-t;
                w*=wn;
            }
        }
    }
    for(int i=0;i<n&&f==-1;i++)F[i]/=n;
    for(int i=0;i<n;i++)a[i]=F[i];
}

int n1,n2,n3,n;

int main(){
    while(~scanf("%d%d%d",&n1,&n2,&n3)&&(n1||n2||n3)){
        n=1;
        while(n<=n1*1+n2*2+n3*5+10)n<<=1;
        for(int i=0;i<n;i++)c1[i]=c2[i]=c3[i]=0;
        for(int i=0;i<=n1;i++)c1[i*1]=1;
        for(int i=0;i<=n2;i++)c2[i*2]=1;
        for(int i=0;i<=n3;i++)c3[i*5]=1;

        fft(c1,n,1),fft(c2,n,1),fft(c3,n,1);
        for(int i=0;i<n;i++)c1[i]=c1[i]*c2[i]*c3[i];
        fft(c1,n,-1);
        for(int i=0;i<n;i++){
            if(int(c1[i].real()+0.5)==0){
                printf("%d\n",i);
                break;
            }
        }
    }
}