ccpc长春场j题

J题：
对于这个题目我首先考虑的是dp[L]，表示x轴长度为L时的方案总数，然后状态转移时考虑L+1，但是无论如何都会产生一定的后效性。但是，只要求出以L+1点为右边界的圆的所有情况，这样就没有后效性了（不以L+1点为右边界的圆都可以在L中体现出来，所以L+1为右边界的圆是L+1状态的特有圆）而以L+1为右边界的圆只可能影响前10个点的状态，所以最多影响到L-9.所以对dp数组新增一维表示前10个点的状态，假设为j，j是一个10位的整数（1<<10)，j的第i位为0或1表示第i个点是否被包含在圆内（边界是0，题目允许相切）。所以我们枚举所有L+1点右边界的圆（半径为1-5的圆，所以枚举1<<5次，状态为pos，同样是一个10位的整数（1<<10),第i位为1表示左边界为i），再判断这些圆的边界是否和前状态j相交，即判断pos状态的左边界是否和j状态的1相同位置，如果相同即为相交。如果没有相交，则可以从j状态转换到pos状态。状态转移的就是将pos数组上的最大圆全覆盖到j状态形成一个新的状态。如下

代码中，mask数组用来存储最大圆信息。fix数组用来存储初始圆的信息。（如果转移过程中包括初始圆则跳过）

#include<iostream>
#include<stdio.h>
#include<vector>
#include<map>
#include<queue>
#include<stack>
#include<stdlib.h>
#include<string>
#include<cstring>
#include<algorithm>
using namespace std;
#define MAX 200005
#define ull unsigned long long
typedef long long ll;
const ll mo = 1e9 + 7;
const double eps = 1e-8;
const int maxn = 1e6 + 3;
const int MK = 1023;
 
int n, m;
int pos[40],mask[40];
 
vector<int> fix[1003];
 
ll dp[1003][1 << 10];
 
 
int main()
{
    cin >> n >> m;
 
    int c, r;
    for (int i = 0; i < m; i++)
    {
        cin >> c >> r;
        c++;
        fix[c + r].push_back(2 * r);
    }
    memset(mask, -1, sizeof(mask));
    for (int i = 1; i <= 31; i++)
    {
        int tmp=0;
        for (int j = 4; j >= 0;j--)
        {
            tmp = (tmp * 2 + ((i >> j) & 1)) * 2;
            if (mask[i] == -1 && ((i >> j) & 1))
            {
                mask[i] = (1 << (j * 2 + 1)) - 1;
            }
        }
        pos[i] = tmp;
    }
    dp[0][MK] = 1;
    for (int i = 0; i <= n; i++)
    {
        for (int j = 0; j <= MK; j++)
        {
            if (fix[i+1].size() == 0)
            {
                dp[i + 1][(j << 1) & MK] += dp[i][j];
                dp[i + 1][(j << 1) & MK] %= mo;
            }
            for (int k = 1; k <= 31; k++) {
                bool flag = 0;
                for (int to : fix[i + 1]) {
                    if (!((pos[k] >> (to - 1)) & 1)) 
                        flag = 1;
                    if (flag) 
                        break;
                }
                if (flag || (j & pos[k]) != 0) continue;
                dp[i + 1][(((j | mask[k]) << 1)) & MK] += dp[i][j];
                dp[i + 1][(((j | mask[k]) << 1)) & MK] %= mo;
            }
 
        }
 
    }
    ll ans = 0;
    for (int i = 0; i <= MK; i++)
    {
        ans += dp[n + 1][i];
        ans %= mo;
    }
 
    cout << ans << endl;