FZU 1912 Divisibility by Thirty-six

题目链接：https://vjudge.net/problem/FZU-1912

题目大意：给定一个长度小于1000的数字串，选取任意个数字组成一个新的数字串，输出能被36整除的最大整数，最大整数不含前导0。不存在这样的整数则输出impossible。

题目思路：

①

数字很大，所以考虑能被36整除的特性。36的约数有2,3,4,6,9,12,18，约数除数的特性的并集就是36当除数的特性。

2：个位必定是偶数

3：各位之和一定是3的倍数

证明： abcd=1000a+100b+10c+d
=999a+99b+9c+a+b+c+d
因为999可被3整除，99能被3整除，9能被3整除
所以a+b+c+d能被3整除的话abcd就能被3整除

4：后两位的数一定是4的倍数 

证明：100能被4整除，所以num能被4整除只需num%100能被4整除。

6：2∪3 

9：各位之和一定是9的倍数 ，证明同3

12：3∪4 

18：2∪9

整理一下得到36当除数的性质，即后两位能被4整除并且各位之和是9的倍数，其他约数的性质都包括在这两个里面，所以一定正确。

 

②

考虑暴力求解，遍历所有能被4整除，且存在于数字串的两个数字作为后两位。得到数字串中各位数的总和，DFS暴力剔除，直到剩下的余数为0，记为一种答案。考虑如何记录最大的答案。首先，位数多的答案一定大。所以剔除的数字越少得到的答案越大。若剔除的位数相同，必然是字典序大的数字大。maxcnt[i]表示最大的答案数字为i的个数，每次遍历完所有情况判断更新即可。

如何反驳位数多可能会有前导0，但位数少可能没有前导0答案反而更大？ 

前导0的存在必定是因为maxcnt[1~9]=0，若是位数多的时候会有前导0，则位数少的时候一定也只会是只有0。

 

③

考虑什么情况下会出现前导0。

这里先占个坑，这里写的比较繁琐，有空了在来整理。

 

④

考虑到样例很多，所以剪枝极为重要。

（重要）cut 1：若是后两位是00，如果不存在某个数字串能被36整除，则其他数字充当后两位也一定不存在。如果存在的话，那么其他数字充当后两位的产生的能被36整除的数字串一定没有00来的大。

第一句话的证明：0和9是唯一两个一定不用剪的数字，因为他不影响和值是9的倍数。因此00放在后两位不动去遍历所有情况都得不到被36整除的串，那么其他数字放到后面一定不行。

第二句话的证明：00放在后两位去剔除其他位的可能性最多，因为把形如12,16的数放在后面，就会不能剔除1和2，1和6。所以说可能性最多一定会得到剔除次数最少的结果。同时00也是最小的两个数字，就算把其他数字放在后两位得到剔除次数和00一样多的结果，那字典序也一定没有00放在后面来的大。

cut 2：当余数是sum的时候，直接去找有没有sum，有的话就剔除这一个，不用考虑其他。 

cut 3：题目要求得到的数字最大，因此后两位是40，04只用取40,80和08只用取80。

cut 4：题目要求得到的数字最大，因此位数多的数字一定大，因此当当前剔除的次数超过最少剔除次数，就return。

 

代码：

#include <stdio.h>
#include <iostream>
#include <cstring>
typedef long long ll;
using namespace std;
int maxtwo;
int maxnum[10],cnt[10];
int num[22]={0,96,92,88,84,80,76,72,68,64,60,56,52,44,40,36,32,28,24,20,16,12};
int cut,pd,two;
void dfs(int sum,int depth)
{
    if(depth>cut)   return ;
    if(sum==0)
    {
        pd=1;
        if(cut==depth)
        {
            int flag=0;
            for(int i=9;i>=0;i--)
            {
                if(maxnum[i]<cnt[i])
                {
                    flag=1;
                    break;
                }
                else if(maxnum[i]>cnt[i])
                {
                    break;
                }
            }
            if(flag)
            {
                for(int i=0;i<=9;i++)    maxnum[i]=cnt[i];
                maxtwo=two;
            }
        }
        else
        {
            cut=depth;
            for(int i=0;i<=9;i++)  maxnum[i]=cnt[i];
            maxtwo=two;
        }
        return ;
    }
    if(cnt[sum])
    {
        cnt[sum]--;
        dfs(0,depth+1);
        cnt[sum]++;
    }
    else
    {
        for(int i=8;i>0;i--)
        {
        if(cnt[i])
            {
                cnt[i]--;
                dfs((sum-i+9)%9,depth+1);
                cnt[i]++;
            }
        }
    }
    return ;
}
int main()
{
    int t,ct=1,sum,a,b;
    cin>>t;
    string s;
    while(t--)
    {
        a=b=0;
        sum=0,pd=0,maxtwo=-1;
        cut=1005;
        memset(cnt,0,sizeof(cnt));
        memset(maxnum,-1,sizeof(maxnum));
        printf("Case %d: ",ct++);
        cin>>s;
        for(int i=0;i<s.size();i++)
        {
            sum=sum+(s[i]-'0');
            cnt[s[i]-'0']++;
            if(s[i]=='0')   a=1;
        }
        sum=sum%9;
        if(cnt[0]>1)  two=0,cnt[0]-=2,dfs(sum,0),cnt[0]+=2;
        if(pd)
        {
            for(int i=9;i>=1;i--)
            {
                if(maxnum[i])   b=1;
            }
            if(b==0)    printf("0");
            else
            {
                for(int i=9;i>=0;i--)
                {
                    for(int j=maxnum[i];j>0;j--)  printf("%d",i);
                }
                printf("00");
            }
        }
        else
        {
            for(int i=1;i<22;i++)
            {
                two=num[i];
                if(cnt[num[i]%10]&&cnt[num[i]/10])
                {
                    cnt[num[i]%10]--;
                    cnt[num[i]/10]--;
                    if(cnt[num[i]%10]>=0)   dfs(sum,0);  //if判断形如44 88是否存在两个
                    cnt[num[i]%10]++;
                    cnt[num[i]/10]++;
                }
            }
            if(pd)
            {
                    for(int i=9;i>=0;i--)
                    {
                        for(int j=maxnum[i];j>0;j--)  printf("%d",i);
                    }
                    printf("%d",maxtwo);

            }
            if(!pd&&a)  printf("0");
            if(!pd&&!a) printf("impossible");
        }
        printf("\n");
    }
    return 0;
}