【学习笔记】[AGC030F] Permutation and Minimum

套路题。然而思维太混乱了没有做出来 我是丝薄

显然有 ( − 1 , x i ) , ( − 1 , − 1 ) (-1,x_i),(-1,-1) (1,xi),(1,1)两种情况。那么定义 v x i = 1 v_{x_i}=1 vxi=1,限制等价于对于 v i = 1 v_i=1 vi=1的两个元素不能在同一组中。

先不考虑算重。显然需要记录 ( − 1 , x i ) , ( − 1 , y i ) (-1,x_i),(-1,y_i) (1,xi),(1,yi)的数目,其中 x i x_i xi是固定的, y i y_i yi是自己选的,每次加入一个元素,可以放到原有的组中或者新开一个组,这样我们可以算出集合 { B i } \{B_i\} {Bi}的数目。

如果从前往后做的话会比较棘手,这也是这题最坑的地方。那么考虑在 d p dp dp的过程中计算贡献,可以从后往前做,对于 ( − 1 , x i ) (-1,x_i) (1,xi)这样的组我们关心每一组 B i B_i Bi的取值,而对于 ( − 1 , y i ) (-1,y_i) (1,yi)这样的组显然我们只关心值组成的集合,最后乘 m ! m! m!即可。

转移过程不再赘述。复杂度 O ( n 3 ) O(n^3) O(n3)

#include<bits/stdc++.h> #define fi first #define se second #define ll long long #define pb push_back #define inf 0x3f3f3f3f using namespace std; const int mod=1e9+7; int n,m,M,a[605],v[605],b[605],f[305][305],g[305][305]; void add(int &x,int y){ if((x+=y)>=mod)x-=mod; } signed main(){ cin>>n,n<<=1;for(int i=1;i<=n;i++)cin>>a[i]; for(int i=1;i<=n/2;i++){ if(a[2*i-1]!=-1&&a[2*i]!=-1)v[a[2*i-1]]=v[a[2*i]]=1; else if(a[2*i-1]!=-1)v[a[2*i-1]]=2; else if(a[2*i]!=-1)v[a[2*i]]=2; else M++; } for(int i=1;i<=n;i++){ if(v[i]!=1)b[++m]=i,v[m]=(v[i]==2); } reverse(b+1,b+1+m),reverse(v+1,v+1+m); f[0][0]=1; for(int i=1;i<=m;i++){ memset(g,0,sizeof g); for(int j=0;j<=n/2;j++){ for(int k=0;k<=n/2;k++){ if(f[j][k]){ if(v[i]){ add(g[j+1][k],f[j][k]); if(k)add(g[j][k-1],f[j][k]); } else{ add(g[j][k+1],f[j][k]); if(k)add(g[j][k-1],f[j][k]); if(j)add(g[j-1][k],(ll)f[j][k]*j%mod); } } } }memcpy(f,g,sizeof g); }ll res(f[0][0]); for(int i=1;i<=M;i++)res=res*i%mod; cout<<res; }

__EOF__

本文作者仰望星空的蚂蚁
本文链接https://www.cnblogs.com/cqbzly/p/17530055.html
关于博主:评论和私信会在第一时间回复。或者直接私信我。
版权声明:本博客所有文章除特别声明外,均采用 BY-NC-SA 许可协议。转载请注明出处!
声援博主:如果您觉得文章对您有帮助,可以点击文章右下角推荐一下。您的鼓励是博主的最大动力!
posted @   仰望星空的蚂蚁  阅读(7)  评论(0编辑  收藏  举报  
相关博文:
阅读排行:
· 全程不用写代码,我用AI程序员写了一个飞机大战
· MongoDB 8.0这个新功能碉堡了,比商业数据库还牛
· 记一次.NET内存居高不下排查解决与启示
· 白话解读 Dapr 1.15:你的「微服务管家」又秀新绝活了
· DeepSeek 开源周回顾「GitHub 热点速览」
点击右上角即可分享
微信分享提示