自己是有多么sb。
密钥
大家都说这是一道普及-的题,一年前我做不起,我可以说我太弱啦,我就普及组水平,今年我还是做不起……
看大佬题解都是:开个桶就好啦!
我:你在说什么……
首先把环拉成链,倍长。
如果确定$i$这个位置是起始位置,那么特征值就是$\sum\limits_{j=1}^{n-1} (p_j!=0 , sum(A_{i+1}...A_{i+j})>0) $。
那么我们先记录一个前缀和,后面所提到的$A$都是前缀和。$\sum\limits_{j=1}^{n-1} (p_j!=0 , A_{i+j} > A_{i})$。
我们可以这样来考虑,我们把$A$看作高度,从左到右连起来看成折线,我每次在$A_{i}$高度处拉根线,比线高的,并且向上走的(就是p!=0),的折线的个数就是特征值
因为$A_i-A_{i-1}$要么是1要么是-1,所以我们每次减掉或者加上某一个取值的个数就可以了,所以这就是开个桶的意义。
然后如果$p!=0$是B的位置的时候,如果图不变,相当于比线低的,往下走的,的折线的个数,就是特征值,
然后因为折线起始点和终止点的高度都是拉的这根线的高度,
所以比线低的,往下走的,和比线低或者相等的,往上走的,是可以一一对应的,数目是相等的,
而比线低或者相等的,往上走的个数,就是k-比线高的,并且向上走的个数
我在说啥……
//Serene
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cstdio>
#include<cmath>
using namespace std;
#define ll long long
#define db double
#define For(i,a,b) for(int i=(a);i<=(b);++i)
#define Rep(i,a,b) for(int i=(a);i>=(b);--i)
const int maxn=4e7+7;
int n,k,seed,now,S,p[maxn],cnt[maxn],A[maxn],ans1,ans2,ans3;
char cc; ll ff;
template<typename T>void read(T& aa) {
aa=0;cc=getchar();ff=1;
while((cc<'0'||cc>'9')&&cc!='-') cc=getchar();
if(cc=='-') ff=-1,cc=getchar();
while(cc>='0'&&cc<='9') aa=aa*10+cc-'0',cc=getchar();
aa*=ff;
}
int getrand(){
seed = ((seed * 12321) ^ 9999) % 32768;
return seed;
}
void generateData(){
scanf("%d%d%d",&k,&seed,&S);
int t = 0;
n = k * 2 + 1;
memset(p, 0, sizeof(p));
for (int i = 1; i <= n; i++)
{
p[i] = (getrand() / 128) % 2;
t += p[i];
}
int i = 1;
while (t > k)
{
while (p[i] == 0)
i++;
p[i] = 0;
t--;
}
while (t < k)
{
while (p[i] == 1)
i++;
p[i] = 1;
t++;
}
}
int main() {
generateData();
For(i,1,n) p[i+n]=p[i];
A[0]=n+1; For(i,1,n<<1) A[i]=A[i-1]+(p[i]? 1:-1);
For(i,1,n) if(p[i]) ++cnt[A[i]],now+=A[i]>A[0];
For(i,1,n) {
if(p[i]) {
now-=cnt[A[i]];//from A[i]-1 to A[i]
--cnt[A[i]];
++cnt[A[i+n]];
now+=A[i+n]>A[i];
}
else {
now+=cnt[A[i]+1];//from A[i]+1 to A[i]
if(now==0) ans1=i;
if(now==S) ans2=i;
if(now==k-S) ans3=i;
}
}
printf("%d\n%d\n%d\n",ans1,ans2,ans3);
return 0;
}
吉夫特
想了很久,还是只会n^2暴力,翻开洛谷题解:
what? 众所周知,只有我不知道?
知道这玩意就可以直接每次枚举子集啦,洛谷上是可以过的,但是bz上就要T掉啦。
//Serene
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cstdio>
#include<cmath>
using namespace std;
#define ll long long
#define db double
#define For(i,a,b) for(int i=(a);i<=(b);++i)
#define Rep(i,a,b) for(int i=(a);i>=(b);--i)
const int maxn=3e5+7;
ll mod=1e9+7;
ll n,a[maxn],p[maxn],dp[maxn];
char cc; ll ff;
template<typename T>void read(T& aa) {
aa=0;cc=getchar();ff=1;
while((cc<'0'||cc>'9')&&cc!='-') cc=getchar();
if(cc=='-') ff=-1,cc=getchar();
while(cc>='0'&&cc<='9') aa=aa*10+cc-'0',cc=getchar();
aa*=ff;
}
int main() {
read(n);
For(i,1,n) read(a[i]),p[a[i]]=i;
For(i,1,n) {
dp[0]+=dp[i]%=mod;
for(int j=a[i];j;j=(j-1)&a[i])
if(p[j]>i) dp[p[j]]+=dp[i]+1;
}
printf("%lld\n",dp[0]%mod);
return 0;
}
所以应该怎么办呢
我们如果枚举子集复杂度大约是3^18的,但是如果我们把18位分成两半,各9位,然后一个是求答案时枚举补集,一个是更新dp数组,枚举子集,就优秀多了
具体意思见代码
//Serene
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cstdio>
#include<cmath>
using namespace std;
#define ll long long
#define db double
#define For(i,a,b) for(int i=(a);i<=(b);++i)
#define Rep(i,a,b) for(int i=(a);i>=(b);--i)
const int maxn=3e5+7,maxt=(1<<9)+7,U=(1<<9)-1;
ll mod=1e9+7;
ll n,o,x,y,dp[maxt][maxt],ans,now;
char cc; ll ff;
template<typename T>void read(T& aa) {
aa=0;cc=getchar();ff=1;
while((cc<'0'||cc>'9')&&cc!='-') cc=getchar();
if(cc=='-') ff=-1,cc=getchar();
while(cc>='0'&&cc<='9') aa=aa*10+cc-'0',cc=getchar();
aa*=ff;
}
void mo(ll& x) {if(x>=mod) x-=mod;}
int main() {
read(n);
For(i,1,n) {
read(o); x=o>>9; y=o&U; now=0;
for(int j=x;j<=U;j=(j+1)|x) now+=dp[j][y];
now%=mod; ans+=now; now++; mo(now);
for(int j=y;;j=(j-1)&y) {
dp[x][j]+=now,mo(dp[x][j]);
if(j==0) break;
}
}
printf("%lld\n",ans%mod);
return 0;
}
来源:oschina
链接:https://my.oschina.net/u/4384347/blog/3988806