hdu 4038 模拟+贪心 Stone

题意：

        给一个数列，可以对这数列进行两种操作：A.使其中的一个数加1；B.增加1到这个数列。要求这个数列的乘积最大。

思路：贪心+模拟。

我觉得不够清晰。再理一遍。

     1. 负数个数为奇数的时候我们把最大的负数变成0；

     2.0是我们在这题当中很讨厌的东西（其他题也有卡0的，比如0！= 1，1对1的逆元是0），所以我们先将所有的精力集中对付0。

     3.当0没有了的时候，我们其次讨厌1。为什么呢，因为1->2是2倍的收益，相对其他的数以及其他的添加1的操作来说，她用1的收益最大。我们要得到最大收益，所以，把1全部变成2。

     4.没有1,0的时候，我们把2全部变成3。是3/2的收益。

     5.没有2,1,0的时候，我们是不是要把所有的3变成4呢？不是的。为什么这里就不是，上面就是呢？是这样的。当操作m=2时，我们用这两个二将两个1变成2收益是4 > 2，将两个2变成3收益是1.5*1.5 > 2，将两个3变成4的收益是1.333333*1.333333 = 1.777777 < 2，还不如将这两个一加在一起变成2加到原来的序列里面；当m=3时，我们用这3个二将3个1变成2收益是8 > 3，将3个2变成3收益是1.5*1.5 *1.5 > 3，将两个3变成4的收益是1.333333*1.333333 * 1.333333 < 3；以此类推。

     6.由以上推得，做完以上操作之后如果m = 1，则加到最小的那个非负数上；m>1时，如果m%3==0,则加m / 3个3就行；(m - 1) % 3 == 0时，加(m-1)/3 - 1个3，再加个4；(m - 2) % 3 == 0，加(m - 2)/3个3，再加个2。

到此，终于完了。。。。

做题过程：

       将思路理清楚之后就容易了，调了一会，1A。

#include <iostream>
#include <cstdio>
#include <algorithm>
#define mod 1000000007
using namespace std;
int t,n,a[100010],odd_num,m,indx,Minpos,Min;
__int64 ans;
void init(){
    odd_num = 0;    ans = 1;
    Min = mod;  indx = 1;//初始化！！！
}
long long pow_mod(int a,int n,int p)
{
    long long ret=1;
    long long A=a;
    while(n)
    {
        if (n & 1)
            ret=(ret*A)%p;
        A=(A*A)%p;
        n>>=1;
    }
    return ret;
}
int main()
{
    scanf("%d",&t);
    for(int ca = 1; ca <= t; ca ++){
        scanf("%d%d",&n,&m); init();
        for(int i = 1; i <= n;i ++){
            scanf("%d",a + i);
            if(a[i] < 0) odd_num++;
        }
        sort(a + 1, a + 1 + n);
        if(odd_num & 1){//负数个数为奇数的时候，去改最大的负数
            int j = 1;
            while(a[j] < 0) j ++;   indx = -- j;//
//            cout << "j" << j << endl;
            if(a[j] + m > 0) m += a[j], a[j] = 0;
            else a[j] += m,m = 0;
        }
        if(m > 0){
//            for(int j = 1; j <= n - 1; j ++)
//                printf("%d ",a[j]);  printf("%d\n",a[n]);
//            cout << indx << endl;
            for(int j = indx; j <= n && m > 0; j ++)
                if(a[j] == 0) a[j] ++, m --;
            for(int j = indx; j <= n && m > 0; j ++)
                if(a[j] == 1) a[j] ++, m --;
            for(int j = indx; j <= n && m > 0; j ++)
                if(a[j] == 2) a[j] ++, m --;

//            for(int j = 1; j <= n - 1; j ++)
//                printf("%d ",a[j]);  printf("%d\n",a[n]);

            for(int j = 1; j <= n; j ++)
                ans = (ans * a[j]) % mod;
            if(m <= 0){
                printf("Case %d: %I64d\n",ca,ans);
            }else{
                if(m == 1){
                    for(int j = 1; j <= n; j ++)
                        if(Min > a[j]) Min = a[j], Minpos = j;
                    ans /= a[Minpos];
                    ans = (ans * (a[Minpos] + 1)) % mod;
                    printf("Case %d: %I64d\n",ca,ans);
                }else if(m % 3 == 0){
                    ans = (ans * pow_mod(3,m/3,mod)) % mod;
                    printf("Case %d: %I64d\n",ca,ans);
                }else if((m - 1) % 3 == 0){
                    ans = (ans * pow_mod(3,(m - 1)/3 - 1,mod)) % mod;
                    ans = (ans * 4) % mod;
                    printf("Case %d: %I64d\n",ca,ans);
                }else if((m - 2) % 3 == 0){
                    ans = (ans * pow_mod(3,(m - 2)/3,mod)) % mod;
                    ans = (ans * 2) % mod;
                    printf("Case %d: %I64d\n",ca,ans);
                }
            }
        }
        else{
            for(int j = 1; j <= n; j ++)
                ans = (ans * a[j]) % mod;
            printf("Case %d: %I64d\n",ca,ans);
        }
    }
    return 0;
}