BZOJ 2614(梦的困境-桶排数组Dp)

2614: 梦的困境

Time Limit: 3 Sec   Memory Limit: 128 MB
Submit: 30   Solved: 16
[ Submit][ Status]

Description

     “每一个OIer上辈子都是折翼的天使……”。

    OIer 喜欢做梦，尤其是关于OI题目的梦。

    TT 同学便是折翼天使的典范。在NOI2010 day1的前夜，TT梦见自己被困在一个房间里，房间里摆放着n个箱子。他还看到墙上贴着一张告示。作为折翼天使的典范，他就像读题一样开始浏览告示。只见告示上写道：

     天使大牛，在你面前摆放着n个魔法箱，每个魔法箱里各有m个碟子，每个碟子上写着一个正整数。然后我们还会给你两个整数k和l。

     现在要求你从这n个魔法箱中各取出1个或2个碟子(取碟子是有顺序的)，但是请注意：由于某种原因，你从第1个箱子和第n个箱子中都只能各取出一个碟子。

     你取出的所有碟子上的数的乘积必须能被k整除。

     令T:=满足上述条件的取碟子方案数 mod l。(mod为取模符号)

     为了能够解开这个房间大门的密码锁，而这个密码便是T。

    TT 急了，这明摆着不就是一道OI题么。哥明天还得考NOI呢，今晚就要刷题啦？醒醒啊，醒醒啊……

     可就是醒不来，他走到门边，发现上面又有一张告示：如果密码输入错误，你明天比赛必然挂掉……

     请你写一个程序帮助TT求出密码T。

 

Input

。第一行包含两个整数n、m。第二行包含两个整数k、l。接下来给出n×m的矩阵。(x，y)元表示第x个箱子中的第y个碟子上的整数。(数与数之间以空格隔开。)

 

Output

一个整数T。

 

Sample Input

3 3
12 100
5 2 1
2 1 2
3 7 4

Sample Output

12

HINT

Source

O(nm)的转移太坑了。。。

由于题目要求被k整除的数，所以只有k的因子部分是有用的。

这样状态就很小。

转移O(m)继续坑-注意到只有gcd(a[i][j],k)决定的转移的方向，可把它们相同的压为一个点。

这样转移就很小。

&取2个数的转移也能轻松解决。。。

找最大因子数dabiao程序：

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<functional>
#include<iostream>
#include<cmath>
#include<cctype>
#include<ctime>
using namespace std;
#define For(i,n) for(int i=1;i<=n;i++)
#define Fork(i,k,n) for(int i=k;i<=n;i++)
#define Rep(i,n) for(int i=0;i<n;i++)
#define ForD(i,n) for(int i=n;i;i--)
#define RepD(i,n) for(int i=n;i>=0;i--)
#define Forp(x) for(int p=pre[x];p;p=next[p])
#define Lson (x<<1)
#define Rson ((x<<1)+1)
#define MEM(a) memset(a,0,sizeof(a));
#define MEMI(a) memset(a,127,sizeof(a));
#define MEMi(a) memset(a,128,sizeof(a));
#define INF (2139062143)
#define F (100000007)
#define MAXN (200+10)
#define MAXM (10000+10)
#define MAXK (200000+10)
#define MAXL (30000+10)
long long mul(long long a,long long b){return (a*b)%F;}
long long add(long long a,long long b){return (a+b)%F;}
long long sub(long long a,long long b){return (a-b+(a-b)/F*F+F)%F;}
typedef long long ll;
int n;
int main()
{
//	freopen("bzoj2614_dabiao.out","w",stdout);
	int maxtot=0;
	For(n,200000)
	{
		int tot=0;
		For(i,sqrt(n))
		{
			if (n%i==0) tot++;
			if (n%i==0&&i*i!=n) tot++;
		}
		maxtot=max(tot,maxtot);
	}
	cout<<maxtot<<endl;
	
	return 0;
}

实际代码：

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<functional>
#include<iostream>
#include<cmath>
#include<cctype>
#include<ctime>
#include<map>
using namespace std;
#define For(i,n) for(int i=1;i<=n;i++)
#define Fork(i,k,n) for(int i=k;i<=n;i++)
#define Rep(i,n) for(int i=0;i<n;i++)
#define ForD(i,n) for(int i=n;i;i--)
#define RepD(i,n) for(int i=n;i>=0;i--)
#define Forp(x) for(int p=pre[x];p;p=next[p])
#define Lson (x<<1)
#define Rson ((x<<1)+1)
#define MEM(a) memset(a,0,sizeof(a));
#define MEMI(a) memset(a,127,sizeof(a));
#define MEMi(a) memset(a,128,sizeof(a));
#define INF (2139062143)
#define MAXN (200+10)
#define MAXM (10000+10)
#define MAXK (200000+10)
#define MAXL (30000+10)
#define MAXK_factor (160+10)
typedef long long ll;
map<ll,int> h;
int n,m,k,F;
long long mul(long long a,long long b){return (a*b)%F;}
long long add(long long a,long long b){return (a+b)%F;}
long long sub(long long a,long long b){return (a-b+(a-b)/F*F+F)%F;}
ll gcd(ll a,ll b){if (b==0) return a;return gcd(b,a%b);}
ll f[MAXN][MAXK_factor]={0},a[MAXN][MAXK_factor]={0},c[MAXK_factor],tot=0;
ll a2[MAXN][MAXK_factor]={0};
int main()
{
//	freopen("bzoj2614.in","r",stdin);
	scanf("%d%d%d%d",&n,&m,&k,&F);
	
	For(i,sqrt(k)) 
		if (k%i==0) 
		{
			c[++tot]=i;
			if (i*i<k) c[++tot]=k/i;
		}
	sort(c+1,c+1+tot);
	//For(i,tot) cout<<c[i]<<' ';cout<<endl;
	For(i,tot) h[c[i]]=i;
	
	For(i,n)
	{
		For(l,m)
		{
			int j;scanf("%d",&j);j=h[gcd(j,k)];
			a[i][j]++;
		}
	}
	//Fork(i,2,n-1) a[i][1]++;
	memcpy(a2,a,sizeof(a2));
	Fork(i,2,n-1)
	{
		For(j,tot)
			Fork(l,j,tot)
			{
				if (!a[i][j]||!a[i][l]) continue;
				if (j==l&&a[i][j]<2) continue;
				int newj=h[gcd(c[j]*c[l],k)];
				ll v;
				if (j^l) v=a[i][j]*a[i][l]*2;else v=a[i][j]*(a[i][j]-1);
				a2[i][newj]+=v;
			}
	}
	//For(j,tot) a2[1][j]=a[1][j],a2[n][j]=a[n][j];
	
	/*
	For(i,n)
	{
		For(j,tot) cout<<a2[i][j]<<' ';cout<<endl;
	}puts("");
	*/
	
	
	For(i,tot) f[1][i]=a[1][i]%F;
	For(i,n-1) //i->i+1
	{
		For(j,tot)
		{
			For(l,tot)
			{
				if (!f[i][j]||!a2[i+1][l]) continue;
				int newj=h[gcd(c[j]*c[l],k)];
				f[i+1][newj]=add(f[i+1][newj],mul(f[i][j],a2[i+1][l]));
			}
		}
	}
	
	/*
	For(i,n)
	{
		For(j,tot) cout<<f[i][j]<<' ';cout<<endl;
	}*/
	
	cout<<f[n][tot]<<endl;
	
	
	return 0;
}