BZOJ3746 : [POI2015]Czarnoksiężnicy okrągłego stołu
NOIP前做了几道POI,现在终于能在BZOJ上提交了…
交上去最后几个点WA,看了数据发现p=0的特判错了…
p=0,1时特判
p=2时构造两种情况判断
p=3时不考虑1的座位进行DP
可以发现对于i+1的位置安排,我们只关心i-2,i-1,i的相对顺序以及它们的相邻、边界情况
所以设f[i][j][S1][S2]表示已经安排了前i个人的座位,i-2,i-1,i的顺序为j,是否有人在两端点为S1,是否有人相邻为S2的方案数
答案最后再除以n
这样复杂度有点飞…
#include<cstdio>
typedef long long ll;
const int N=1000010,P=1000000007;
int n,m,p,i,j,S1,S2,ans;bool ban[N][7];//0:-3 3:0 6:3
inline int abs(int x){return x>0?x:-x;}
inline void read(int&a){char c;while(!(((c=getchar())>='0')&&(c<='9')));a=c-'0';while(((c=getchar())>='0')&&(c<='9'))(a*=10)+=c-'0';}
namespace Subtask2{
int a[N];
bool flag;
inline void Main(){
if(n&1){
flag=1;
for(a[j=1]=1,i=2;i<n;i+=2)a[++j]=i;
for(i=n;i>2;i-=2)a[++j]=i;
for(i=1;i<n;i++)if(ban[a[i]][a[i+1]-a[i]+3])flag=0;
if(ban[a[n]][a[1]-a[n]+3])flag=0;
if(flag)ans++;
flag=1;
for(a[j=0],i=1;i<=n;i+=2)a[++j]=i;
for(i=n-1;i;i-=2)a[++j]=i;
for(i=1;i<n;i++)if(ban[a[i]][a[i+1]-a[i]+3])flag=0;
if(ban[a[n]][a[1]-a[n]+3])flag=0;
if(flag)ans++;
}else{
flag=1;
for(a[j=1]=1,i=2;i<=n;i+=2)a[++j]=i;
for(i=n-1;i>2;i-=2)a[++j]=i;
for(i=1;i<n;i++)if(ban[a[i]][a[i+1]-a[i]+3])flag=0;
if(ban[a[n]][a[1]-a[n]+3])flag=0;
if(flag)ans++;
flag=1;
for(a[j=0],i=1;i<n;i+=2)a[++j]=i;
for(i=n;i;i-=2)a[++j]=i;
for(i=1;i<n;i++)if(ban[a[i]][a[i+1]-a[i]+3])flag=0;
if(ban[a[n]][a[1]-a[n]+3])flag=0;
if(flag)ans++;
}
}
}
namespace Subtask3{
struct E{int nxt,j,S1,S2,l,r;/*1:i-2 2:i-1 3:i*/}e[200];
int g[6][4][4],ed,k;
int loc[6][3]={{1,2,3},{1,3,2},{2,1,3},{2,3,1},{3,1,2},{3,2,1}};
int hash[500],tmp[10],r,pos[3],pl[5];
int f[2][6][4][4];
bool isl,isr,cl1,cl2;
inline ll pow(ll a,int b){ll t=1;for(;b;b>>=1,a=a*a%P)if(b&1)t=t*a%P;return t;}
inline void add(int l,int r){
int nj=0,i,t=0,p;
for(i=l;i<=r;i++)if(tmp[i]>=2)pl[++t]=i,(nj*=10)+=tmp[i];
e[++ed].j=hash[nj];
e[ed].S1=((pl[1]==l)<<1)|(pl[3]==r);
e[ed].S2=((pl[1]+2>=pl[2])<<1)|(pl[2]+2>=pl[3]);
for(t=0,i=l;i<=r;i++)if(tmp[i]){
pl[++t]=i;
if(tmp[i]==4)p=t;
}
if(p>1&&pl[p-1]+2>=pl[p])e[ed].l=tmp[pl[p-1]];
if(p<4&&pl[p]+2>=pl[p+1])e[ed].r=tmp[pl[p+1]];
e[ed].nxt=g[j][S1][S2];
g[j][S1][S2]=ed;
}
inline void init(){
hash[243]=1;hash[324]=2;hash[342]=3;hash[423]=4;hash[432]=5;
for(j=0;j<6;j++)for(S1=0;S1<4;S1++)for(S2=0;S2<4;S2++){
isl=S1>>1&1,isr=S1&1,cl1=S2>>1&1,cl2=S2&1;
tmp[r=0]=0;
if(!isl)tmp[++r]=0;
tmp[pos[1]=++r]=loc[j][0];
tmp[++r]=0;if(!cl1)tmp[++r]=0;
tmp[pos[2]=++r]=loc[j][1];
tmp[++r]=0;if(!cl2)tmp[++r]=0;
tmp[++r]=loc[j][2];
if(!isr)tmp[++r]=0;
if(isl&&isr)tmp[0]=4,add(0,r),tmp[0]=0,tmp[r+1]=4,add(1,r+1),tmp[r+1]=0;
if(cl1)tmp[pos[1]+1]=4,add(1,r),tmp[pos[1]+1]=0;
if(cl2)tmp[pos[2]+1]=4,add(1,r),tmp[pos[2]+1]=0;
}
}
inline void Main(){
init();
for(j=0;j<6;j++)f[1][j][3][3]=1;
for(i=3;i<n;i++){
for(j=0;j<6;j++)for(S1=0;S1<4;S1++)for(S2=0;S2<4;S2++)f[i&1^1][j][S1][S2]=0;
for(j=0;j<6;j++)for(S1=0;S1<4;S1++)for(S2=0;S2<4;S2++)if(f[i&1][j][S1][S2])for(k=g[j][S1][S2];k;k=e[k].nxt){
if(e[k].l==1)if(ban[i-2][6])continue;
if(e[k].r==1)if(ban[i+1][0])continue;
if(S2>>1&1){
if(loc[j][0]==1&&ban[i-2][loc[j][1]+2]&&e[k].l!=1)continue;
if(loc[j][1]==1&&ban[i+loc[j][0]-3][4-loc[j][0]]&&e[k].r!=1)continue;
}
if(S2&1){
if(loc[j][1]==1&&ban[i-2][loc[j][2]+2]&&e[k].l!=1)continue;
if(loc[j][2]==1&&ban[i+loc[j][1]-3][4-loc[j][1]]&&e[k].r!=1)continue;
}
if(S1==3){
if(loc[j][0]==1&&ban[i+loc[j][2]-3][4-loc[j][2]]&&e[k].l&&e[k].r)continue;
if(loc[j][2]==1&&ban[i-2][loc[j][0]+2]&&e[k].l&&e[k].r)continue;
}
if(S1&1){
if(loc[j][2]==1&&ban[i-2][6]&&(!e[k].l||!e[k].r))continue;
}
if(S1>>1&1){
if(loc[j][0]==1&&ban[i+1][0]&&(!e[k].l||!e[k].r))continue;
}
(f[i&1^1][e[k].j][e[k].S1][e[k].S2]+=f[i&1][j][S1][S2])%=P;
}
}
for(j=0;j<6;j++)for(S1=0;S1<4;S1++)for(S2=0;S2<4;S2++)if(f[n&1][j][S1][S2]){
if(S2>>1&1)if(ban[n+loc[j][0]-3][loc[j][1]-loc[j][0]+3])continue;
if(S2&1)if(ban[n+loc[j][1]-3][loc[j][2]-loc[j][1]+3])continue;
if(S1==3)if(ban[n+loc[j][2]-3][loc[j][0]-loc[j][2]+3])continue;
(ans+=f[n&1][j][S1][S2])%=P;
}
ans=(ll)ans*pow(n,P-2)%P;
}
}
int main(){
read(n),read(m),read(p);
if(n==1)return puts("1"),0;
if(p==0)return puts("0"),0;
if(n==2)return puts(m?"0":"1"),0;
if(p==1)return puts("0"),0;
while(m--){
read(i),read(j);
if(abs(i-j)<=p)ban[i][j-i+3]=1;
}
if(p==2)Subtask2::Main();else Subtask3::Main();
return printf("%d",ans),0;
}
BZOJ3746 : [POI2015]Czarnoksiężnicy okrągłego stołu的更多相关文章
- BZOJ #3746: [POI2015]Czarnoksiężnicy okrągłego stołu 动态规划
转载请注明出处:http://www.cnblogs.com/TSHugh/p/8823423.html 读完题就会发现p=0.1的情况以及n=1.2的情况都可以直接判掉,而p=2的时候也可以直接构造 ...
- BZOJ 4276: [ONTAK2015]Bajtman i Okrągły Robin [线段树优化建边]
4276: [ONTAK2015]Bajtman i Okrągły Robin 题意:\(n \le 5000\)个区间\(l,r\le 5000\),每个区间可以选一个点得到val[i]的价值,每 ...
- [ONTAK2015]Bajtman i Okrągły Robin
bzoj 4276: [ONTAK2015]Bajtman i Okrągły Robin Time Limit: 40 Sec Memory Limit: 256 MB Description 有 ...
- Bajtman i Okrągły Robin
Bajtman i Okrągły Robin 题目描述 你是一个保安,你发现有n个强盗,其中第i个强盗会在[a[i],a[i]+1],[a[i]+1,a[i]+2],...,[b[i]-1,b[i] ...
- 4276: [ONTAK2015]Bajtman i Okrągły Robin
4276: [ONTAK2015]Bajtman i Okrągły Robin Time Limit: 40 Sec Memory Limit: 256 MBSubmit: 345 Solved ...
- 【BZOJ4276】[ONTAK2015]Bajtman i Okrągły Robin 线段树优化建图+费用流
[BZOJ4276][ONTAK2015]Bajtman i Okrągły Robin Description 有n个强盗,其中第i个强盗会在[a[i],a[i]+1],[a[i]+1,a[i]+2 ...
- BZOJ_4276_[ONTAK2015]Bajtman i Okrągły Robin_线段树优化建图+最大费用最大流
BZOJ_4276_[ONTAK2015]Bajtman i Okrągły Robin_线段树优化建图+最大费用最大流 Description 有n个强盗,其中第i个强盗会在[a[i],a[i]+1 ...
- BZOJ4276 : [ONTAK2015]Bajtman i Okrągły Robin
建立线段树, S向每个叶子连边,容量1,费用0. 孩子向父亲连边,容量inf,费用0. 每个强盗向T连边,容量1,费用为c[i]. 对应区间内的点向每个强盗,容量1,费用0. 求最大费用流即可. #i ...
- BZOJ 4276: [ONTAK2015]Bajtman i Okrągły Robin
最大权值匹配,贪心匈牙利即可. 检查一些人是否能被全部抓住可以采用左端点排序,右端点优先队列处理. By:大奕哥 #include<bits/stdc++.h> using namespa ...
随机推荐
- unity StreamingAssets路径
原地址:http://blog.csdn.net/nateyang/article/details/8493791 我们在读写例如XML和TXT文件的时候,在电脑上和手机上路径不一致,造成了很多麻烦, ...
- Centos镜像使用帮助
https://lug.ustc.edu.cn/wiki/mirrors/help/centos
- MAC的OS X10.10更新以后进入用户界面就死机
我用的是搜狗输入法,所以把搜狗卸载就好了.(注意是卸载,不是单纯的从输入源里移除) 下载一个搜狗输入法的DMG,打开后选择卸载搜狗输入法. 转自: http://zhidao.baidu.com/qu ...
- ENGINE=InnoDB
最开始用MySQL Administrator建数据库的时候,表缺省是InnoDB类型,也就没有在意.后来用Access2MySQL导数据的时候发现只能导成 MyISAM类型的表 区别如下原来是MyI ...
- wsp反编译
最后出于好奇,我把wsp文件解压缩,看看里面是什么(如果您的机器上的压缩软件不能直接解压,可尝试修改后缀名为cab.).我看到的首先是一个清单文件(manifest.xml),一个DLL文件(Shar ...
- Telegram
官网:https://www.telegram.org/ GitHub:https://github.com/telegramdesktop/tdesktop
- 【转】Eclipse中查看jar包中的源码
(简单的方式:通过jd-gui来进行反编译,最简单!,参考我的另一篇博文, 地址:http://www.cnblogs.com/gmq-sh/p/4277991.html) Java Decompil ...
- 读取STL模型
读取二进制格式的STL模型文件 std::ifstream fin;fin.open(stlFilePath, std::ios::in | std::ios::binary);bool isBina ...
- ajax(ajax开发与入门)
AJAX即“Asynchronous Javascript And XML”(异步JavaScript和XML),是指一种创建交互式网页应用的网页开发技术.AJAX是一种进行页面局部异步刷新的技术,局 ...
- Android之自定义ViewGroup
概述 在写代码之前,我必须得问几个问题: 1.ViewGroup的职责是啥? ViewGroup相当于一个放置View的容器,并且我们在写布局xml的时候,会告诉容器(凡是以layout为开头的属性, ...