第1题:
perica ##——1794
思路:裸裸的组合公式送分题,然而又是这个mod——wa了。错的原因也很明显,不能先对余再除,而且上次讲过的公式也忘了…但今天的重点不是这,而是逆元!
逆元:简单来说,就是将除法变乘法,减法变加法。如-x=+(-x),1/x=1*x^(-1)。 这里就是要用到有关这个的费马小定理——a^(p-2)=a(mod p)。
#include<bits/stdc++.h>
#define REP(i,f,t) for(int i=(f),i##_end_=(t);i<=i##_end_;i++)
#define LL long long
#define N 100005
#define M 55
using namespace std;
const int P=1e9+7;
int n,k;
int A[N];
LL dp[N][M];
LL ans;
struct P1{
void solve(){
REP(i,0,n)dp[i][0]=1;
REP(i,1,n)
REP(j,1,k)
dp[i][j]=(dp[i-1][j]+dp[i-1][j-1])%P;
REP(i,k,n)ans+=1LL*A[i]*dp[i-1][k-1],ans%=P;
cout<<ans<<endl;
}
}p1;
struct P2{
int Pow(int x,int y) {
int res=1;
while(y) {
if(y&1)
res=(1LL*res*x)%P;
x=(1LL*x*x)%P;
y>>=1;
}
return res;
}
void work(int a,int b){
if(a==b){dp[a][b]=1;return;}
int res=Pow(a-b,P-2);
dp[a][b]=(1LL*dp[a-1][b]*a)%P;
dp[a][b]=(1LL*dp[a][b]*res)%P;
}
void solve(){
REP(i,k,n){
work(i-1,k-1);
ans+=1LL*dp[i-1][k-1]*A[i],ans%=P;
}
cout<<ans<<endl;
}
}p2;
int main(){
srand(time(NULL));
scanf("%d%d",&n,&k);
REP(i,1,n)scanf("%d",&A[i]);
sort(A+1,A+1+n);
int f=rand()%2;
if(f)p1.solve();
else p2.solve();
return 0;
}
第2题:
delete_num ##——2988
思路:又是求方案的玩意,这一题显然是dp,但是这个序列的元素过int,所以不能单纯的存一维,所以马上就想到状压dp。(其实此题也是可以滚动,也很明显,只存在两层之间的关系)。
#include<bits/stdc++.h>
#define REP(i,a,b) for(int i=(a),i_##END_=(b);i<=i_##END_;++i)
#define DREP(i,a,b) for(int i=(a),i_##BEGIN_=(b);i>=i_##BEGIN_;--i)
#define N 105
#define M 100005
#define LL long long
#define INF 0x3f3f3f3f
#define P 1000000007
using namespace std;
LL A;
int n;
int len[N],B[M>>1];
int dp[N][M];
LL num[N],Num[N][M];
LL ans;
LL Change(LL x,int y) {
int L[15],mark[15];
memset(mark,0,sizeof(mark));
int c=0;
while(x)L[++c]=x%10,x/=10;
int i=1,a=1;
while(i<=y){
if(y&i)mark[a]=1;
i<<=1;
a++;
}
LL tmp=0;
DREP(i,c,1){
if(mark[i])continue;
tmp=1LL*(tmp*10+L[i]);
}
return tmp;
}
void solve(){
REP(i,1,n+1)REP(j,0,(1<<len[i])-2)Num[i][j]=Change(num[i],j);
REP(i,1,n+1)sort(Num[i],Num[i]+(1<<len[i])-1);
REP(i,0,(1<<len[1])-2)dp[1][i]=1;
int a=0;
memset(B,0,sizeof(B));
REP(i,0,(1<<len[1])-2){
while(Num[1][i]>Num[2][a]&&a<=(1<<len[2])-1)a++;
if(a==(1<<len[2])-1)break;
B[a]+=dp[1][i];
}
REP(i,0,(1<<len[2])-2)if(i>0)B[i]+=B[i-1],B[i]%=P;
REP(i,2,n+1){
REP(j,0,(1<<len[i])-2)dp[i][j]+=B[j],dp[i][j]%=P;
int b=0;
if(i<n+1){
memset(B,0,sizeof(B));
REP(j,0,(1<<len[i])-2){
while(Num[i][j]>Num[i+1][b]&&b<=(1<<len[i+1])-1)b++;
if(b==(1<<len[i+1])-1)break;
B[b]+=dp[i][j],B[b]%=P;
}
REP(j,0,(1<<len[i+1])-2)if(j>0)B[j]+=B[j-1],B[j]%=P;
}
}
REP(i,0,(1<<len[n+1])-2)ans+=dp[n+1][i],ans%=P;
cout<<ans<<endl;
}
int main() {
cin>>A>>n;
REP(i,1,n+1){
num[i]=A+i-1;
LL k=num[i];
while(k)len[i]++,k/=10;
}
solve();
return 0;
}
第3题:
four ##——2989
思路:枚举中间两个点,再从两点找Mx,其实考试是有想到,但自己的bellman写复杂了,当场否定了这个想法…(同时,今天给小c gay了,bellman其实是卡过去的,它的复杂度是不好证明的(然而我并不造…),遇到链的话就会退化了…)
#include<bits/stdc++.h>
#define REP(i,f,t) for(int i=(f),i##_end_=(t);i<=i##_end_;i++)
#define DREP(i,f,t) for(int i=(f),i##_end_=(t);i>=i##_end_;i--)
#define N 2005
#define M 5005
#define INF 0x3f3f3f3f
#define LL long long
// LL mod 内存 调试 文件名
using namespace std;
int n,m,ans;
void chkmax(int &x,int y){if(x<y)x=y;}
struct node{
int from,to,cost;
}E[M];
struct Node{
int Mx[2],cnt[2],Id[2];
}A[N];
int dis[N][N];
void Bellman(int s){
dis[s][s]=0;
while(1){
bool f=0;
REP(i,1,m)
if(dis[s][E[i].from]!=INF && dis[s][E[i].to]>dis[s][E[i].from]+E[i].cost)dis[s][E[i].to]=dis[s][E[i].from]+E[i].cost,f=1;
if(!f)break;
}
}
int main(){
scanf("%d%d",&n,&m);
REP(i,1,m)scanf("%d%d",&E[i].from,&E[i].to),E[i].cost=1;
memset(dis,INF,sizeof(dis));
REP(i,1,n)Bellman(i);
REP(i,1,n){
REP(j,1,n){
if(i==j || dis[i][j]==INF)continue;
if(dis[i][j]>A[i].Mx[0])A[i].Id[0]=j,A[i].cnt[0]=1,A[i].Mx[0]=dis[i][j];
else if(dis[i][j]==A[i].Mx[0])A[i].cnt[0]++;
if(dis[i][j]>A[j].Mx[1])A[j].Id[1]=i,A[j].cnt[1]=1,A[j].Mx[1]=dis[i][j];
else if(dis[i][j]==A[j].Mx[1])A[j].cnt[1]++;
}
}
REP(i,1,n)REP(j,1,n){
if(dis[i][j]==INF || i==j || !A[i].cnt[1] || !A[j].cnt[0] || (A[i].Id[0]==A[j].Id[1] && A[i].cnt[0]==1 && A[j].cnt[1]==1))continue;
chkmax(ans,A[i].Mx[1]+dis[i][j]+A[j].Mx[0]);
}
cout<<ans<<endl;
return 0;
}
小结 ##:
这次考得不好,还是由于第1题wa了,虽然也是出错在mod,但其实公式跟mod半毛钱关系都没有,而且公式还是讲过的…