bzoj千题计划164：bzoj5123: 线段树的匹配

http://www.lydsy.com/JudgeOnline/upload/201712/prob12.pdf

dp[len][0/1] 表示节点表示区间长度为len，节点选/不选的 最大匹配

sum[len][0/1] 表示对应dp[len][0/1]的方案数

这里选节点即选节点与其父节点的边

设区间长度为len，左子区间长度为L，右子区间长度为R

这个节点选，那么左右子节点都不能选

dp[len][1]=1+dp[L][0]+dp[R][0]

sum[len][1]=sum[L][0]*sum[R][0]

这个节点不选，有3种情况：

左右子节点都不选：

dp[len][0]=dp[L][0]+dp[R][0]

sum[len][0]=sum[L][0]+sum[R][0]

选左子节点：

dp[len][0]=dp[L][1]+dp[R][0]

sum[len][0]=sum[L][1]+sum[R][0]

选右子节点：

dp[len][0]=dp[L][0]+dp[R][1]

sum[len][0]=sum[L][0]+sum[R][1]

如果dp[len][0] 在三种情况中有相同的，sum[len][0]要累加

len虽然是1e18，但只会用log种，所以用map

#include<map>
#include<cstdio>
#include<cstring>

using namespace std;

typedef long long LL;

const int mod=;

//int dp[100001][2];
//int sum[10001][2];

map<LL,LL>dp[];
map<LL,LL>sum[];

void dfs(LL len)
{
    if(len==)
    {
        dp[][]=;
        sum[][]=;
        dp[][]=;
        sum[][]=;
        return;
    }
    LL R=len>>;
    LL L=len-R;
    if(dp[].find(L)==dp[].end()) dfs(L);
    if(dp[].find(R)==dp[].end()) dfs(R);
    dp[][len]=+dp[][L]+dp[][R];
    sum[][len]=sum[][L]*sum[][R]%mod;
    dp[][len]=dp[][len]-;
    sum[][len]=sum[][len];
    if(dp[][L]+dp[][R]>dp[][len])
    {
        dp[][len]=dp[][L]+dp[][R];
        sum[][len]=sum[][L]*sum[][R]%mod;
    }
    else if(dp[][L]+dp[][R]==dp[][len])
    {
        sum[][len]=(sum[][len]+sum[][L]*sum[][R])%mod;
    }
    if(dp[][L]+dp[][R]>dp[][len])
    {
        dp[][len]=dp[][L]+dp[][R];
        sum[][len]=sum[][L]*sum[][R]%mod;
    }
    else if(dp[][L]+dp[][R]==dp[][len])
    {
        sum[][len]=(sum[][len]+sum[][L]*sum[][R])%mod;
    }
}

int main()
{
    LL n;
    scanf("%I64d",&n);
    dfs(n);
    printf("%I64d %I64d",dp[][n],sum[][n]);
}