【XSY1285】【BZOJ3814】【清华集训2014】简单回路 状压DP

最新推荐文章于 2021-09-28 10:30:51 发布
最新推荐文章于 2021-09-28 10:30:51 发布
阅读量647 收藏 1
点赞数
分类专栏： DP DP--状压DP
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/ez_yww/article/details/78155949
版权

题目描述

　　给你一个 n×m 的网格图和 k 个障碍，有q个询问，每次问你有多少个不同的不经过任何一个障碍点且经过 (x,y) 与 (x+1,y) 之间的简单回路

　　 n≤1000,m≤6,q≤10000

题解

　　简单插头DP

　　先用DP求出前面 i 行的的方案数和后面i行的方案数，询问时暴力合并。

　　其实状态数很少的（不到 200 ）

　　时间复杂度： O(???)

代码

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cstdlib>
#include<ctime>
#include<utility>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int,int> pii;
ll p=1000000007;
int c[100010];
int d[1010];
int ld[1010][10];
int rd[1010][10];
int cnt;
int n,m,k,q;
int la[10];
int ra[10];
int set(int s,int x,int v)
{
    s&=~(3<<(2*(x-1)));
    s|=v<<(2*(x-1));
    return s;
}
int get(int s,int x)
{
    return (s>>(2*(x-1)))&3;
}
int st[10];
void dfs(int x,int s)
{
    if(x>m+1)
    {
        int i;
        int now=0;
        for(i=1;i<=m+1;i++)
        {
            if(get(s,i)==1)
                st[++now]=i;
            else if(get(s,i)==2)
            {
                la[i]=st[now];
                ra[st[now]]=i;
                now--;
            }
            if(now<0)
                return;
        }
        if(now>0)
            return;
        for(i=1;i<=m+1;i++)
            fprintf(stderr,"%d",get(s,i));
        fprintf(stderr,"\n");
        d[++cnt]=s;
        c[s]=cnt;
        memcpy(ld[cnt],la,sizeof la);
        memcpy(rd[cnt],ra,sizeof ra);
        return;
    }
    int i;
    for(i=0;i<=2;i++)
        dfs(x+1,set(s,x,i));
}
void add(ll &x,ll y)
{
    x=(x+y)%p;
}
struct dp
{
    int a[1010][10];
    ll f[1010][10][200];
    void solve()
    {
        f[1][1][1]=1;
        int i,j,k;
        for(i=1;i<=n;i++)
        {
            for(j=1;j<=m;j++)
                for(k=1;k<=cnt;k++)
                {
                    if(!f[i][j][k])
                        continue;
                    int s=d[k];
                    int l=get(s,j);
                    int r=get(s,j+1);
                    if(a[i][j])
                    {
                        if(!l&&!r)
                            add(f[i][j+1][k],f[i][j][k]);
                        continue;
                    }
                    if(!l)
                        if(!r)
                        {
                            add(f[i][j+1][k],f[i][j][k]);
                            int v=set(s,j,1);
                            v=set(v,j+1,2);
                            add(f[i][j+1][c[v]],f[i][j][k]);
                        }
                        else if(r==1)
                        {
                            add(f[i][j+1][k],f[i][j][k]);
                            int v=set(s,j,1);
                            v=set(v,j+1,0);
                            add(f[i][j+1][c[v]],f[i][j][k]);
                        }
                        else
                        {
                            add(f[i][j+1][k],f[i][j][k]);
                            int v=set(s,j,2);
                            v=set(v,j+1,0);
                            add(f[i][j+1][c[v]],f[i][j][k]);
                        }
                    else if(l==1)
                    {
                        if(!r)
                        {
                            add(f[i][j+1][k],f[i][j][k]);
                            int v=set(s,j,0);
                            v=set(v,j+1,1);
                            add(f[i][j+1][c[v]],f[i][j][k]);
                        }
                        else if(r==1)
                        {
                            int v=set(s,j,0);
                            v=set(v,j+1,0);
                            v=set(v,rd[k][j+1],1);
                            add(f[i][j+1][c[v]],f[i][j][k]);
                        }
                        else
                        {
//                          fprintf(stderr,"orzzjt\n");
                        }
                    }
                    else
                    {
                        if(!r)
                        {
                            add(f[i][j+1][k],f[i][j][k]);
                            int v=set(s,j,0);
                            v=set(v,j+1,2);
                            add(f[i][j+1][c[v]],f[i][j][k]);
                        }
                        else if(r==1)
                        {
                            int v=set(s,j,0);
                            v=set(v,j+1,0);
                            add(f[i][j+1][c[v]],f[i][j][k]);
                        }
                        else
                        {
                            int v=set(s,j,0);
                            v=set(v,j+1,0);
                            v=set(v,ld[k][j],2);
                            add(f[i][j+1][c[v]],f[i][j][k]);
                        }
                    }
                }
            for(j=1;j<=cnt;j++)
                if(!get(d[j],m+1))
                    add(f[i+1][1][c[d[j]<<2]],f[i][m+1][j]);
        }
    }
};
dp a,b;
int tot;
int e1[10];
int e2[10];
int op[1010];
int f[1010];
struct list
{
    int h[210];
    int v[100010];
    int t[100010];
    int n;
    list()
    {
        n=0;
        memset(h,0,sizeof h);;
    }
    void add(int x,int y)
    {
        n++;
        v[n]=y;
        t[n]=h[x];
        h[x]=n;
    }
};
list e;
int find(int x)
{
    return f[x]==x?x:f[x]=find(f[x]);
}
int check(int x,int y)
{
    if(x==53&&y==53)
        int xxx=1;
    int i;
    for(i=1;i<=m+1;i++)
    {
        bool x1=get(d[x],i);
        bool y1=get(d[y],i);
        if(x1^y1)
            return 0;
    }
    for(i=1;i<=m+1;i++)
        f[i]=i;
    int fa;
    for(i=1;i<=m+1;i++)
        if(get(d[x],i))
        {
            fa=find(i);
            if(get(d[x],i)==1)
                f[find(i)]=find(rd[x][i]);
            if(get(d[y],i)==1)
                f[find(i)]=find(rd[y][i]);
        }
    for(i=1;i<=m+1;i++)
        if(get(d[x],i)&&find(i)!=fa)
            return 0;
    return 1;
}
int main()
{
    freopen("d2t1.in","r",stdin);
    freopen("d2t1.out","w",stdout);
    scanf("%d%d%d",&n,&m,&k);
    int j,l;
    int i,x,y;
    for(i=1;i<=k;i++)
    {
        scanf("%d%d",&x,&y);
        a.a[x][y]=b.a[n-x+1][y]=1;
    }
    dfs(1,0);
    for(i=1;i<=cnt;i++)
        for(j=1;j<=cnt;j++)
            if(check(i,j))
                e.add(i,j);
    a.solve();
    b.solve();
    scanf("%d",&q);
    for(i=1;i<=q;i++)
    {
        scanf("%d%d",&x,&y);
        ll ans=0;
        for(j=1;j<=cnt;j++)
        {
            int s=d[j];
            if(get(s,m+1))
                continue;
            if(!get(s,y))
                continue;
            int bo=1;
            tot=0;
            for(l=1;l<=m;l++)
                if(get(s,l))
                {
                    if((a.a[x][l]||a.a[x+1][l]))
                    {
                        bo=0;
                        break;
                    }
                    tot++;
                    e1[tot]=l;
                    e2[tot]=get(s,l);
                }
            if(!bo)
                continue;
            for(l=e.h[j];l;l=e.t[l])
                ans=(ans+a.f[x+1][1][c[s<<2]]*b.f[n-x+1][1][c[d[e.v[l]]<<2]])%p;
        }
        printf("%lld\n",ans);
    }
    return 0;
}
【BZOJ3814】【UOJ39】【清华集训2014】简单回路
cz_xuyixuan的博客
02-28 569
【题目链接】BZOJUOJ【思路要点】直接进行\(Q\)次基于连通性的动态规划可以获得60分。注意到一定需要经过的边只有竖直的，而没有水平的。我们可以在完成一行的动态规划时顺便记录在第1行到第\(i\)行的区域中，第\(i\)行的插头为\(Mask\)的方案数以及第\(i\)行到第\(N\)行的区域中，第\(i\)行的插头为\(Mask\)的方案数。在询问时将经过某条边的配对的插头对应的方案数相乘...
【XSY1536】【BZOJ3522】【BZOJ4543】【POI2014】Hotel 树形DP 长链剖分 启发式合并
ez_yww的博客
08-12 659
题目描述​　　给你一棵树，求有多少个组点满足x≠y,x≠z,y≠z,distx,y=distx,z=disty,zx\neq y,x\neq z,y\neq z,dist_{x,y}=dist_{x,z}=dist_{y,z}
【清华集训2014】 uoj #40 卡常数 - KD-Tree
Mys_C_K的博客
10-02 509
总之就是各种乱七八糟的KDT卡常即可。。。 #include&amp;amp;lt;iostream&amp;amp;gt; #include&amp;amp;lt;cstring&amp;amp;gt; #include&amp;amp;lt;cstdio&amp;amp;gt; #include&amp;amp;lt;algorithm&amp;amp;gt; #include&amp;amp;lt;cmath&amp;am
清华集训2014 day2 task1 简单回路
weixin_34221276的博客
01-10 140
题目 如题。 算法 就是刚学习的插头DP。 从前往后和从后往前分别进行一次DP。 要点 合法的括号序列只有103个 如何合并两次dp的信息 一开始犯傻了，以为当且仅当两个轮廓线的状态相同才是合法的方案。其实很容易举出反例。 如果直接枚举的话，每次询问的时间复杂度是\(O(103^2 m)\)。 为了加快速度，可以把所有合法的方案先列举出来（就是预处理），只有\(103^2\)个。每次询问的复杂度优...
uoj#164. 【清华集训2015】V
weixin_30617561的博客
04-26 134
传送门：http://uoj.ac/problem/164 思路：科学的题面： 请你写一个数据结构支持以下功能： 1：区间[l,r]加x 2：区间[l,r]减x并和0取max 3：区间覆盖 4：单点询问 5：单点历史最大值询问 线段树维护分段函数 标记就是一个二元组(a,b)表示标记生效后x=max(x+a,b) 1操作就是打(x,0)的标记 2就是(-x,0...
【uoj#37/bzoj3812】[清华集训2014]主旋律 状压dp+容斥原理
weixin_30901729的博客
03-30 135
题目描述 求一张有向图的强连通生成子图的数目对 $10^9+7$ 取模的结果。 题解 状压dp+容斥原理 设 $f[i]$ 表示点集 $i$ 强连通生成子图的数目，容易想到使用总方案数 $2^{sum[i]}$ 减去不为强连通图的方案数得到强连通图的方案数，其中 $sum[i]$ 表示点集 $i$ 中边的数目。 考虑什么样的图不是强连通图：缩点后入度为0的强连通分量对应的点...
【XSY2536】【BZOJ2655】calc DP 数学 拉格朗日插值
ez_yww的博客
08-16 1451
题目大意​　　一个序列a1,…,ana_1,\ldots,a_n是合法的，当且仅当： ​　　长度为给定的nn。 ​　　a1,…,ana_1,\ldots,a_n都是[1,m][1,m]中的整数。 ​　　a1,…,ana_1,\ldots,a_n互不相等。 ​　　一个序列的值定义为它里面所有数的乘积，即a1×a2×⋯×ana_1\times a_2\times\cdots\times a_n。
OFDM.rar_XSY_ofdm qpsk awgn_ofdm+qpsk_两径_两径信道
09-21
《OFDM技术在QPSK调制下的性能分析——基于AWGN及多径衰落信道》 OFDM（Orthogonal Frequency Division Multiplexing，正交频分复用）是一种高效的数字调制技术，广泛应用于现代无线通信系统，如4G LTE、5G NR等。...
[CSA35G][XSY3318]Counting Quests （DP）
Edward The Bunny 的博客
09-28 241
XSY3318 CSA35G 对于一个询问区间的集合SSS，求出每一个数被哪些区间覆盖了，记为SiS_iSi​。 要能保证猜出选中数，当且仅当每个数的SiS_iSi​互不相同。 考虑求出不满足要求的集合SSS的个数。 首先可以观察得到SiS_iSi​的一个性质：若a<b<c<da<b<c<da<b<c<d且Sa=Sc,Sb=SdS_a=S_c,S_b=S_dSa​=Sc​,Sb​=Sd​，则必有Sa=Sb=Sc=SdS_a=S_b=S_c=S_dSa​=
[XSY] PQ树（区间DP）
Edward The Bunny 的博客
03-09 477
PQ树 对于非叶子节点，它只能调换它的儿子的顺序，却不能调换整颗子树的顺序 合法的PQ树都能表示出1,2,3,…,n 这个排列 有以上两点可以得出：PQ树的每棵子树对应的一定都是一段连续的数字区间 考虑区间DP，设f(l,r)f(l,r)f(l,r)表示数字l~r构成一棵PQ树的方案数 数字l~r 能够构成一棵子树 当且仅当 数字l~r 在题目给出的所有排列中都在一段连续区间中，我们定义这样的 [l,r] 是“好的” 若[l,r]不是"好的"，f(l,r)=0f(l,r)=0f(l,r)=0 若[l,r]
UOJ#37.【清华集训2014】主旋律 状压DP
NeighThorn的博客
10-09 1383
状压DP+容斥原理
【LOJ2323】「清华集训 2017」小 Y 和地铁
cz_xuyixuan的博客
09-29 545
【题目链接】 点击打开链接 【思路要点】 很不错的脑洞题。 附上官方题解。 时间复杂度 O(T∗2N2)O(T*2^{\frac{N}{2}})O(T∗22N​) 。 【代码】 #include&lt;bits/stdc++.h&gt; using namespace std; const int MAXN = 105; const int INF = 1e9; typede...
【清华集训2014】【BZOJ3817】Sum
CreationAugust is 14 years old forever
02-06 2293
Description给定正整数N，R。求 Input第一行一个数 T，表示有 T 组测试数据。 接下来 T 行，每行两个正整数 n,r。 Output输出 T 行，每行一个整数表示答案。 Sample Input33 53 63 7 Sample Output31-1 HINT对于 100% 的数据，满足 n≤10^9,r≤10^4,T≤10^4。Source2015年国家集训队测
CF11D：A Simple Task（状压dp & 图）
junior19的博客
07-31 687
D. A Simple Task time limit per test 2 seconds memory limit per test 256 megabytes input standard input output standard output Given a simple graph, output the number of
poj 2288 Islands and Bridges解题报告-状态压缩dp
weixin_30908103的博客
04-17 105
题目链接：http://poj.org/problem?id=2288 题目描述：哈密尔顿路问题。n个点，每一个点有权值，设哈密尔顿路为C1C2...Cn，Ci的权值为Vi，一条哈密尔顿路的值分为三部分计算： 1.每一个点的权值之和 2.对于图中的每一条CiCi+1，加上Vi*Vi+1 3.对于路径中的连续三个点：CiCi+1Ci+2，若在图中，三点构成三角形，则要加上Vi*Vi+1*V...
[转载]有向图无向图是否有环判断
weixin_30782331的博客
04-20 663
问题：给出一个算法，用它来确定一个给定的无向图G=(V,E)中是否包含一个回路。所给出的算法的运行时间为O(V)，这一时间独立于|E|解答：我们都知道对于一个无向图而言，如果它能表示成一棵树，那么它一定没有回路，并且有|E|=|V|-1，如果在这个树上添加一条边，那么就构成了回路，如果在这个树中去掉一个边就成了森林（注意：如果只是限定|E|<|V|-1它不一定是森林，它当中可能存在无向连通子...
loj#2326. 「清华集训 2017」简单数据结构【动态规划+set】
cdsszjj的博客
05-31 890
传送门 解题思路： 感觉复杂度是证不来的……不过跑得挺快…… 注意序列长度不超过20，所以也比较暴力，就是记录每个值的位置，以每个点开头的最长长度和每个长度对应的最小结尾位置，用set维护以每个点为结尾的最长序列开头位置，用桶维护答案……向左加就枚举倍数更新，向右加就枚举约数更新…… #include&amp;lt;bits/stdc++.h&amp;gt; using namespace std; in...
CodeForces 21D Traveling Graph 状压dp+欧拉回路
九野的博客
07-01 1685
题目链接：点击打开链接 题意： 给定n个点m条边的无向图 求从1点开始经过每条边至少一次最后回到1点的最小路程 显然就是找一条路径可重复的欧拉回路 思路： 首先对于欧拉回路的结论是：所有点的度数都为偶数 因为所有边至少经过一次，那么可以把题意转换成加最少多少条边使得图满足以上结论 而加的边目的是为了把奇度数转成偶度数，先floyd一下得到任意点间加边的最小花费 dp[i]表示状态
UVALive 4887 Soccer 状压+模拟
九野的博客
09-05 1213
题目链接：点击打开链接 题意：n个球队，m场比赛 下面n行表示n个球队的名字。 下面m场比赛表示该场比赛的2个队得分。 -1表示我们可以任意填。 这种任意填的比赛场数不超过12场。 求： 胜一场球队得2分，平得1分，败得0分。 求每个球队最好名次与最差名字。 每场只有3个状态，最多只有12场，所以状压一下，3^12个状态。 #include #include
http://localhost/api/xsy-auth-security/loginUser/login 403 (Forbidden)
最新发布
07-28
- *2* [c/c++使用libcurl库做http客户端及封装(HTTP_GET和HTTP_POST)](https://blog.csdn.net/xsy29000/article/details/103181267)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值