Submission #367786

#	Time	Username	Problem	Language	Result	Execution time	Memory
367786		arnold518	Semafor (COI20_semafor)	C++14	60 / 100	976 ms	33388 KiB

semafor

#include <bits/stdc++.h>
using namespace std;

typedef long long ll;
typedef pair<int, int> pii;
typedef pair<ll, ll> pll;
typedef vector<vector<ll>> Mat;

const int MAXN = 1500;
const ll MOD = 1e9+7;

const int MM[10]={12, 8, 5, 25, 10, 19, 20, 9, 23, 27};

int M, X;
ll N, K;

Mat operator + (const Mat &p, const Mat &q)
{
	Mat ret;
	int sz=p.size();
	ret=vector<vector<ll>>(sz, vector<ll>(sz));
	for(int i=0; i<sz; i++)
	{
		for(int j=0; j<sz; j++)
		{
			ret[i][j]=0;
			for(int k=0; k<sz; k++)
			{
				ret[i][j]+=p[i][k]*q[k][j]%MOD;
				ret[i][j]%=MOD;
				//if(ret[i][j]>=MOD) ret[i][j]-=MOD;
			}
		}
	}
	return ret;
}

Mat mypow(Mat x, ll y)
{
	if(y==0)
	{
		int sz=x.size();
		Mat ret=vector<vector<ll>>(sz, vector<ll>(sz, 0));
		for(int i=0; i<sz; i++) ret[i][i]=1;
		return ret;
	}
	if(y%2) return mypow(x, y-1)+x;
	Mat ret=mypow(x, y/2);
	return ret+ret;
}

Mat A, B, C;
Mat P[MAXN+10], Q, R, S, PP;

ll comb[MAXN+10][MAXN+10];
ll ans[100];

int main()
{
	scanf("%d%lld%lld%d", &M, &N, &K, &X);

	A=vector<vector<ll>>(32, vector<ll>(32));
	for(int i=0; i<32; i++)
	{
		for(int j=0; j<32; j++)
		{
			if(__builtin_popcount(i^j)==1) A[i][j]=1;
		}
	}

	if(M==2)
	{
		comb[0][0]=1;
		for(int i=1; i<=K; i++)
		{
			for(int j=0; j<=i; j++)
			{
				if(j==0 || j==i) comb[i][j]=1;
				else comb[i][j]=(comb[i-1][j]+comb[i-1][j-1])%MOD;
			}
		}
		
		P[0]=mypow(A, 0);
		for(int i=1; i<=K; i++)
		{	
			P[i]=P[i-1]+A;
		}

		Q=vector<vector<ll>>(100, vector<ll>(100));
		for(int i=0; i<=K; i++)
		{
			for(int p=0; p<100; p++)
			{
				for(int q=0; q<100; q++)	
				{
					Q[p][q]+=P[i][MM[p/10]][MM[q/10]]*P[K-i][MM[p%10]][MM[q%10]]%MOD*comb[K][i]%MOD;
					Q[p][q]%=MOD;
				}
			}
		}

		R=vector<vector<ll>>(100, vector<ll>(100));
		for(int i=0; i<=N%K; i++)
		{
			for(int p=0; p<100; p++)
			{
				for(int q=0; q<100; q++)	
				{
					R[p][q]+=P[i][MM[p/10]][MM[q/10]]*P[N%K-i][MM[p%10]][MM[q%10]]%MOD*comb[N%K][i]%MOD;
					R[p][q]%=MOD;
				}
			}
		}

		S=mypow(Q, N/K)+R;

		for(int i=0; i<100; i++)
		{
			printf("%lld\n", S[X][i]);
		}
		return 0;
	}
	else
	{
		B=mypow(A, K);
		C=mypow(A, N%K);
	 
		PP=vector<vector<ll>>(10, vector<ll>(10));
		Q=vector<vector<ll>>(10, vector<ll>(10));
	 
		for(int i=0; i<10; i++)
		{
			for(int j=0; j<10; j++)
			{
				PP[i][j]=B[MM[i]][MM[j]];
				Q[i][j]=C[MM[i]][MM[j]];
			}
		}
	 
	 
		R=mypow(PP, N/K)+Q;
	 
		for(int i=0; i<10; i++)
		{
			printf("%lld\n", R[X][i]);
		}
	}
}

Compilation message (stderr)

semafor.cpp: In function 'int main()':
semafor.cpp:60:7: warning: ignoring return value of 'int scanf(const char*, ...)', declared with attribute warn_unused_result [-Wunused-result]
   60 |  scanf("%d%lld%lld%d", &M, &N, &K, &X);
      |  ~~~~~^~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

• Subtask 1: 6 / 6
• Subtask 2: 15 / 15
• Subtask 3: 12 / 12
• Subtask 4: 12 / 12
• Subtask 5: 15 / 15
• Subtask 6: 0 / 40