蓝桥杯-动态规划Sereja and Squares

/*
    蓝桥杯上一道中文题，题目意思就不解释了
    最开始看到题目的时候，虽然已经标明是动态规划了，我还是想暴力深搜一波（明知道会超时）
    写到一半感觉好麻烦，就没搞完。
    既然已经说明是动态规划了，就用动态搞搞呗。接下来说说正确思路。
    题目说擦去了全部的大写字母，和部分小写字母，而且要求大小写配对
    既然是配对就肯定要求字母数是偶数，也就是奇数直接输出0就行了
    而且大小写字母数肯定是一样的为n/2（谁叫是配对呢，嫐、嬲是非法的/doge）
    因为小写字母在前的原因，小写字母可以随便填（x除外），大写根据小写的填就行了
    那么就可以建立一个数组f[i][j]，表示前i个位置，已经填好了j个大写字母
    那么就要考虑f[i][j]是怎么构成的
    当第i个地方是需要填小写字母的话，就由f[i-1][j]过来
    当第i个地方是需要填大写字母的话，就由f[i-1][j-1]过来
    综合一下就是f[i][j]=f[i-1][j]+f[i-1][j-1]
    最后结果就是f[n][n/2]，前n个数填写了n/2个大写字母
    然后我们就可以矩阵快速幂敲代码了
    （没有啦，开个玩笑，n给到了100000根本没法开那么大的二维数组，
    我印象中数组一次性能开的空间在主函数内只能开到1000000，外是100000000）

    没法开那么大的数组的话，就要想别的方法，
    我们发现在上面的表达式里面，右边全部是i-1，左边全部是i
    后面的数把前面的数覆盖了的话对计算没什么影响，
    这时候我们可以考虑数组重复使用，也就是滚动数组
    用一层循环去表示i，用一维数组f[j]来表示填好了j个大写字母
    这样大写字母的填写方式就计算完了

    小写字母的就比较简单了，因为小写字母一定是n/2个，
    输入不是问号就是小写字母（我们就用p表示输入的小写字母）
    那么小写字母就需要填n/2-p个，再加上25个字母随便填
    也就是(n/2-p)^25
    最后和大写字母的一乘取模就就是结果了，但是这样还是只有90分
    ((n/2-p)^25*f[n][n/2])%mod

    这是因为大小写字母相乘可能直接爆long long了所以要用快速乘
    计算小写字母的时候也会爆，为了简单用了快速乘加循环来解决快速幂

    后面看了别的大佬的解析，不用long long用unsigned int可以不用取模
    这样简单是简单了，但是我这个菜鸡根本想不到，还是快速乘好用/doge
*/
#include <iostream>
#include <string.h>
#include <stdio.h>
#include <math.h>
#define ll long long
#define mod 4294967296
using namespace std;

ll f[200005]= {1},ans=1;
char c[100005];

ll quick(ll a,ll n)
{
    ll sum=0;
    while(n>=1)
    {
        if(n&1)
            sum=(sum+a)%mod;
        a=(a+a)%mod;
        n=n>>1;
    }
    return sum;
}
int main()
{
    int n;
    scanf("%d",&n);
    getchar();
    for(int i=1; i<=n; i++)
        scanf("%c",&c[i]);
//    for(int i=0;i<n;i++)
//        cout<<c[i]<<' ';
    if(n%2)
        printf("0\n");
    else
    {
        int m=n/2;
        int p=0;
        for(int i=1; i<=n; i++)
        {
            if(c[i]=='?')
            {
                for(int j=i/2; j>=i-m&&j; j--)
                {
                    f[j]=(f[j]+f[j-1])%mod;
                    //cout<<f[j]<<endl;
                }
            }
            else
                p++;
        }
        for(int i=0; i<n/2-p; i++)
        {
//            ans=ans*25;
//            ans=ans%mod;
            ans=quick(ans,25);
        }
        ans=quick(ans,f[n/2]);
//        ans=ans*f[n/2];
//        ans=ans%mod;
        cout<<ans<<endl;
    }
    return 0;
}
/*
4
a???

4
abc?

6
abc???
*/