有点难度啊。。。
A Minimum Cost Perfect Matching
https://www.nowcoder.com/acm/contest/145/A
队友搞的 应该不难 。。。。把二进制写出来看看
#include<bits/stdc++.h>
using namespace std;
#define LL long long int
#define lson rt<<1,l,m
#define rson rt<<1|1,m+1,r
int res[500500];
int s[30];
int main()
{
int n;
scanf("%d",&n);
s[0]=1;
for(int i=1;i<=30;i++)
s[i]=s[i-1]<<1;
for(int i=n-1;i>0;i--)
{
if(!res[i])
{
int x=upper_bound(s,s+30,i)-s;
x=s[x];
x=x-i-1;
if(!res[x]) res[i]=x,res[x]=i;
}
}
for(int i=0;i<n;i++)
{
if(i) printf(" ");
printf("%d",res[i]);
}
printf("\n");
return 0;
}
J Sudoku Subrectangles
https://www.nowcoder.com/acm/contest/145/J
先处理每个元素向右向下扩张的最大范围,因为最多只有52所以只有1000*1000*104的复杂度。
然后就把每一个元素作为一个子矩阵的左上角,计算一下符合的子矩阵有多少即可。。
注意ans需要是long long 。。。。
#include<bits/stdc++.h>
using namespace std;
char s[1111][1111];
int heng[1111][1111];
int shu[1111][1111];
int vis[222];
int zs[1111][1111];
int main(){
int n,m;
scanf("%d%d",&n,&m);
for(int i=1;i<=n;i++)
{
scanf("%s",s[i]+1);
}
memset(vis,0,sizeof(vis));
int pc=1;
for(int i=1;i<=n;i++)
{
for(int j=1;j<=m;j++)
{
pc++;
int k;
for(k=j;k<=min(m,j+51);k++)
{
if(vis[s[i][k]]==pc)
{
break;
}
vis[s[i][k]]=pc;
}
heng[i][j]=k-j;
pc++;
for(k=i;k<=min(n,i+51);k++)
{
if(vis[s[k][j]]==pc)
{
break;
}
vis[s[k][j]]=pc;
}
shu[i][j]=k-i;
}
}
long long ans=0;
for(int i=1;i<=n;i++)
{
for(int j=1;j<=m;j++)
{
ans+=heng[i][j];
zs[i][j]=shu[i][j];
for(int k=j+1;k<=j+heng[i][j]-1;k++)
{
zs[i][k]=min(zs[i][k-1],shu[i][k]);
}
int len=heng[i][j];
for(int k=i+1;k<=i+shu[i][j]-1;k++)
{
len=min(len,heng[k][j]);
while(zs[i][j+len-1]<k-i+1&&len>0)
{
len--;
}
if(len<=0)
break;
ans+=len;
}
}
}
cout<<ans<<endl;
return 0;
}
C Bit Compression
https://www.nowcoder.com/acm/contest/145/C
dfs 如果遇到全变为0的情况,在这三种操作下是不会再出现1了,所以遇到这种情况就不用再继续dfs了。加上这个剪枝就好了。
#include<bits/stdc++.h>
using namespace std;
char s[(1<<18)];
int da[(1<<19)];
int ans=0;
void dfs(int x)
{
if(x<0)
{
if(da[1]==1)
ans++;
return;
}
int len=(1<<x);
int js=0;
for(int i=0;i<len;i++)
{
int wz=i+len;
da[wz]=da[wz*2]&da[wz*2+1];
if(da[wz]==0)
js++;
}
if(js!=len)
dfs(x-1);
js=0;
for(int i=0;i<len;i++)
{
int wz=i+len;
da[wz]=da[wz*2]|da[wz*2+1];
if(da[wz]==0)
js++;
}
if(js!=len)
dfs(x-1);
js=0;
for(int i=0;i<len;i++)
{
int wz=i+len;
da[wz]=da[wz*2]^da[wz*2+1];
if(da[wz]==0)
js++;
}
if(js!=len)
dfs(x-1);
}
int main(){
int n;
cin>>n;
cin>>s;
int len=(1<<n);
for(int i=0;i<len;i++)
da[i+len]=s[i]-'0';
dfs(n-1);
cout<<ans<<endl;
return 0;
}
E Counting 4 Cliques
构造。
转自:https://blog.csdn.net/hhhhhhxh/article/details/81543926
首先考虑构造一个完全图,那么t个点的完全图一共能有C(t, 4)个大小为4的团。但是在C(t, 4)和C(t+1, 4)之间会有空缺,因此在完全图外放若干个点,每个点与这个完全图中的若干个点连边,最后会形成类似于C(t, 4)+C(x1, 3)+C(x2, 3)+…这样的式子,并且要能补满中间的空隙。
打表发现5个点可以满足。(不知道怎么严格证明)
#include<bits/stdc++.h>
using namespace std;
#define ll long long
int C[2222222];
ll C4(ll x)
{
ll ans=x*(x-1)*(x-2)*(x-3)/(1*2*3*4);
return ans;
}
ll C3(ll x)
{
ll ans=x*(x-1)*(x-2)/(1*2*3);
return ans;
}
int main(){
memset(C,-1,sizeof(C));
for(ll i=1;i<=100;i++)
{
ll ans=i*(i-1)*(i-2)/(1*2*3);
C[ans]=i;
}
ll n;
cin>>n;
ll lim=70;
ll t=4;
while((t+1)<=lim&&C4(t+1)<=n)
{
t++;
}
lim=min(lim,t);
ll a,b,c,d;
ll e=-1;
for(a=2;a<=lim;a++)
{
for(b=a;b<=lim;b++)
{
for(c=b;c<=lim;c++)
{
for(d=c;d<=lim;d++)
{
ll ans=C3(a)+C3(b)+C3(c)+C3(d);
if(ans+C4(t)<=n&&C[n-ans-C4(t)]>=0&&C[n-ans-C4(t)]<=lim)
{
e=C[n-ans-C4(t)];
break;
}
}
if(e!=-1)break;
}
if(e!=-1)break;
}
if(e!=-1)break;
}
long long lss=t*(t-1)/2+a+b+c+d+e;
cout<<t+5<<" "<<lss<<endl;
for(ll i=1;i<=t;i++)
{
for(ll j=i+1;j<=t;j++)
{
cout<<i<<" "<<j<<endl;
}
}
for(ll i=1;i<=a;i++)
{
cout<<t+1<<" "<<i<<endl;
}
for(ll i=1;i<=b;i++)
{
cout<<t+2<<" "<<i<<endl;
}
for(ll i=1;i<=c;i++)
{
cout<<t+3<<" "<<i<<endl;
}
for(ll i=1;i<=d;i++)
{
cout<<t+4<<" "<<i<<endl;
}
for(ll i=1;i<=e;i++)
{
cout<<t+5<<" "<<i<<endl;
}
return 0;
}