HJ28 素数伴侣

看到讨论说这是匈牙利算法，多目标跟踪。

求解最大匹配问题

 带你入门多目标跟踪（三）匈牙利算法&KM算法 - 知乎 (zhihu.com)

'''
这题配对很难理解，理解了一天，看程序理解很乱，还是得看输出，逐个跟着。
寻找最大配对数方法如下：
在保持现有配对数不变的情况下，继续寻找可行配对。虽然配对数目不变，但是配对的具体元素改变。
以奇数列表为外循环。
为保证递归不进入死循环，需要两列表对偶数做记号
visited=[]对已经遍历过的偶数下标做标记，
match=[]记录与该下标偶数相加为素数的奇数。
模拟过程如下，
在for循环里为每个元素提供一个初始化列表；
元素visit后下标记录为1，
    若该下标没被奇数占领，记录该新元素奇数；若改下标被占领记录，把占领元素向后寻位置，
    如果有位置则替换位置；若无位置，该奇数元素向后寻位置。
    （如此遍历，在保留素数对数不减的情况下，寻找更多的素数对数。）
    这样，每个奇数与最后可匹配的偶数匹配
方法很巧妙。
'''
def cprime(item):
    if item<=3:
        return True
    for i in range(2,int(item**0.5)+1):
        if item%i==0:
            return False
    else:
        return True
print(cprime(5))
def find_even(odd):
    for k,even in enumerate(evens):        
        if cprime(odd+even) and visited[k]==0:  
            visited[k]=1
            if match[k]==0 or find_even(match[k]):
                match[k]=odd
                return True
    return False
no=int(input())
l=list(map(int,input().split()))
odds,evens=[],[]
for i in l:
    if i%2==0:
        evens.append(i)
    else:
        odds.append(i)
#odds=[1,2,3,7,9,11]
#evens=[4,6,8,10]
count=0
match=[0]*len(evens)
for odd in odds:
    visited=[0]*len(evens)
    if find_even(odd):
        count+=1
print(count)