[BZOJ3688]折线统计(dp+树状数组)

这篇博客详细介绍了如何使用动态规划(DP)和树状数组来解决BZOJ3688题目中的折线统计问题。博主首先阐述了题目的背景,接着深入解析了解题思路,重点讲解了状态转移方程f(i, j, 0/1),表示划分i段,以第j段结束,且最后一段上升或下降的方案数。最后,博客提供了实现代码,帮助读者理解并应用这种解题方法。" 127241226,10786536,vue3项目配置eslint+prettier教程,"['前端开发', 'Vue', 'JavaScript', 'VSCode', '代码风格']

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

题目:

我是超链接

题解:

f(i,j,0/1)表示划分了i段,以第j个结尾,最后一段是上升0\下降1的方案数

f(i,j,0)=k=i+1j1y(k)<y(j)f(i,k,0)+k=i+1j1y(k)<y(j)f(i1,k,1) f ( i , j , 0 ) = ∑ k = i + 1 j − 1 ∑ y ( k ) < y ( j ) f ( i , k , 0 ) + ∑ k = i + 1 j − 1 ∑ y ( k ) < y ( j ) f ( i − 1 , k , 1 )

f(i,j,1)=k=i+1j1y(k)>y(j)f(i1,k,0)+k=i+1j1y(k)>y(j)f(i,k,1) f ( i , j , 1 ) = ∑ k = i + 1 j − 1 ∑ y ( k ) > y ( j ) f ( i − 1 , k , 0 ) + ∑ k = i + 1 j − 1 ∑ y ( k ) > y ( j ) f ( i , k , 1 )

两个限制条件,第一个边做边维护,第二个用bit维护前缀和,01分开维护,滚动数组
然后就是O(nlogn)了

代码:

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
const int mod=100007;
const int N=50005;
const int inf=100000;
struct hh{int x,y;}p[N];
int cmp(hh a,hh b){return a.x<b.x;}
int f[12][N][2],c[inf+5][2][2];
void add(int loc,int v,int id,int jd)
{
    for (int i=loc;i<=inf;i+=i&(-i))
      c[i][id][jd]=(c[i][id][jd]+v)%mod;
}
int qurry(int loc,int id,int jd)
{
    int ans=0;
    for (int i=loc;i>=1;i-=i&(-i))
      ans=(ans+c[i][id][jd])%mod;
    return ans;
}
int main()
{
    freopen("line.in","r",stdin);
    freopen("line.out","w",stdout);
    int n,k;scanf("%d%d",&n,&k);
    for (int i=1;i<=n;i++) scanf("%d%d",&p[i].x,&p[i].y);
    sort(p+1,p+n+1,cmp);

    for (int i=1;i<=n;i++) f[0][i][0]=f[0][i][1]=1;

    for (int i=1;i<=k;i++)
    {
        memset(c,0,sizeof(c));
        add(p[i].y,f[i-1][i][0],(i-1)&1,0);
        add(p[i].y,f[i-1][i][1],(i-1)&1,1);
        for (int j=i+1;j<=n;j++)
        {
            f[i][j][0]=(f[i][j][0]+qurry(p[j].y-1,i&1,0)+qurry(p[j].y-1,(i-1)&1,1))%mod;
            f[i][j][1]=(f[i][j][1]+qurry(inf,(i-1)&1,0)-qurry(p[j].y-1,(i-1)&1,0)+qurry(inf,i&1,1)-qurry(p[j].y,i&1,1)+mod)%mod;
            //注意这里一减可能是负数,要+mod
            add(p[j].y,f[i][j][0],i&1,0);
            add(p[j].y,f[i-1][j][0],(i-1)&1,0);
            add(p[j].y,f[i][j][1],i&1,1);
            add(p[j].y,f[i-1][j][1],(i-1)&1,1);
        }
    }
    int ans=0;
    for (int i=1;i<=n;i++) ans=(ans+f[k][i][0]+f[k][i][1])%mod;
    printf("%d",ans);
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值